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Abstract 


This  study  examines  the  flow  of  groundwater  through  a 
constructed  treatment  wetland.  The  wetland  was  built  to 
explore  the  viability  of  constructed  wetlands  as  a 
treatment  technology  for  groundwater  contaminated  with 
perchloroethylene,  and  it  employs  an  upward  vertical  flow 
design.  A  major  goal  of  the  study  is  to  determine  whether 
the  system  design  facilitates  uniform  vertical  flows 
through  the  subsurface  soil  sediments  or  if  preferential 
flows  occur.  Conceptually,  uniform  flows  will  achieve  the 
most  efficient  degree  of  contaminant  removal  possible  by 
evenly  dispersing  the  groundwater  contaminants  throughout 
the  full  volume  of  the  subsurface  media. 

A  three-dimensional  grid  of  piezometers  is  used  to 
generate  potentiometric  contour  maps,  and  in-situ  tests  of 
hydraulic  conductivity  facilitate  construction  of  a 
numerical  computer  model.  The  contours  of  hydraulic  head, 
measured  conductivities,  and  numerical  model  simulations 
imply  preferential  flows  and  a  wetland  operating  at  less 
than  an  optimal  level  of  efficiency.  Hydraulic  residence 
times  for  representative  water  particles  released  in  the 
model  range  from  16.5  hours  to  15  days  with  a  mean 
residence  time  of  three  days.  The  divergence  from  the 
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uniform  flow  ideal  suggests  an  alternative  construction 
approach  may  be  appropriate. 
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GROUNDWATER  FLOW  THROUGH  A  CONSTRUCTED  TREATMENT  WETLAND 


I .  Introduction 

The  purpose  of  this  research  is  to  characterize  the 
flow  of  groundwater  through  a  constructed  treatment  wetland. 
The  location,  purpose  and  description  of  the  wetland  are 
covered  below.  Hydraulic  head  measurements  from  a  three- 
dimensional  grid  of  piezometers  will  be  used  to  construct 
potentiometric  contour  maps,  conduct  simulation  modeling, 
and  calculate  hydraulic  residence  times.  An  understanding 
of  the  flow  through  the  wetland  will:  (1)  provide  insight 
into  where  important  wetland  and  contaminant  interactions 
may  be  occurring;  (2)  indicate  the  amount  of  time  that  the 
interactions  have  to  occur;  and  (3)  substantiate  or  dispel 
concerns  over  the  possibility  of  water  preferentially 
flowing  through  the  wetland  media. 

Two  wetland  cells  were  constructed  at  Wright  Patterson 
Air  Force  Base  in  Dayton,  Ohio,  during  the  summer  of  2000. 

The  cells  were  constructed  by  the  Air  Force  Institute  of 
Technology  (AFIT)  to  study  the  ability  of  constructed 
wetlands  to  naturally  degrade  perchloroethylene  (PCE) .  The 
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study  represents  a  joint  effort  between  students  and 
faculty  of  AFIT,  the  Air  Force  Research  Laboratory,  and 
Wright  State  University. 

Both  wetland  cells  are  situated  over  an  aquifer  that 
is  contaminated  with  a  plume  of  PCE .  The  site  has  been 
identified  and  documented  with  the  Ohio  Environmental 
Protection  Agency  with  no  current  requirement  for 
remediation.  The  origin  of  the  plume  is  unknown,  but  dry 
cleaning  operations  are  known  to  have  existed  in  the  area 
(Shelley,  2001) .  The  concentration  of  PCE  in  the  aquifer 
in  the  aqueous  phase  is  approximately  .05  mg/L. 

PCE  -  or  tetrachloroethylene  as  it  is  also  known  -  is 
a  volatile  organic  compound  (VOC)  and  a  prevalent 
groundwater  contaminant;  it  is  among  the  three  most 
frequently  detected  groundwater  contaminants  nationwide 
(National  Research  Council,  1994) .  PCE  is  used  as  a 
solvent  in  applications  such  as  metal  degreasing  and 
commercial  dry  cleaning.  PCE  is  a  potential  carcinogen  and 
regulated  drinking  water  contaminant;  it  has  an  assigned 
maximum  contaminant  level  of  .005  mg/L  (Masters,  1998) . 

A  common  approach  to  treating  groundwater  contaminated 
with  VOCs  is  to  employ  energy  intensive  technologies  such 
as  pump  and  treat  systems.  This  is  because  VOCs  -  by  their 
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nature  -  will  volatilize  at  room  temperature.  Although 
costly,  technologies  that  aerate  contaminated  groundwater 
to  encourage  vaporization  allow  for  a  high  degree  of 
control  during  the  remediation  process  (Masters,  1998). 

Another  remediation  strategy  lies  at  the  other  end  of 
the  energy-use  and  technology  spectrum.  It  is  called 
natural  attenuation,  and  it  is  essentially  letting  nature 
run  its  course.  Natural  attenuation  relies  on  natural 
processes  and  energies  to  degrade  contaminants  into  less 
harmful  substances.  Research  indicates  that,  under  the 
right  conditions,  PCE  will  degrade  as  a  result  of  microbial 
activity  into  less  chlorinated  degradation  products  and 
possibly  even  all  the  way  down  into  innocuous  substances 
(Lee  et  al,  1998) .  Natural  attenuation  may  require  some 
degree  of  human  intervention,  however,  to  create  the  right 
conditions  for  the  contaminant  degradation  to  occur.  For 
example,  a  contaminated  aquifer  could  be  seeded  with 
microorganisms  that  are  known  to  decontaminate  PCE  as  well 
as  a  substrate  to  serve  as  a  source  of  carbon  and  energy 
for  the  microbes  (Fogel  et  al,  1995;  Wu  et  al,  1995) .  One 
disadvantage  with  such  a  method  is  that  it  may  be  difficult 
to  control  the  attenuation  process  and  achieve  the  desired 
remediation  outcome  in  the  natural  setting. 
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The  wetland  cells  at  Wright  Patterson  Air  Force  Base 


represent  a  low  cost,  low  energy,  pump  and  treat  system 
that  relies  on  natural  processes  to  attenuate  PCE .  The 
concept  is  supported  by  findings  at  a  U.S.  Department  of 
Defense  site  in  Aberdeen,  Maryland,  where  a  natural, 
groundwater-fed  wetland  appears  to  demonstrate  the  complete 
destruction  of  trichloroethylene  -  a  related  VOC  and  also  a 
degradation  product  of  PCE  (Lorah  and  Olsen,  1999) . 

Because  the  wetland  cells  are  constructed  above  ground, 
their  design  accommodates  a  greater  degree  of  process 
control  than  that  available  with  natural  attenuation. 

The  design  of  the  wetland  cells  is  unique  to  wetlands 
constructed  for  the  purpose  of  treating  contaminants. 

Both  wetland  cells  are  approximately  120'  long  and  60'  wide, 
and  they  function  in  the  following  manner.  Contaminated 
water  from  the  aquifer  is  pumped  into  the  bottom  of  each 
wetland  cell  and  distributed  through  three  parallel  3" 
diameter  PVC  pipes.  The  pipes  are  perforated  along  the 
sides  and  run  lengthwise  through  a  9"  gravel  layer  along 
the  bottom  of  the  wetland  cells.  (An  impermeable  liner 
isolates  the  water  in  the  cells  from  the  surrounding 
environment.)  The  water  forced  through  these  pipes 
gradually  works  its  way  up  to  the  surface  of  the  wetland 
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cells  through  approximately  54"  of  wetland  sediment.  The 
water  then  flows  across  the  wetland  surface  to  exit  weirs 
located  at  one  end  of  each  wetland  cell.  Plots  of 
different  emergent  wetland  plant  species  grow  in  both 
wetland  cells.  Figure  1-1  is  a  sketch  of  the  wetland  cell 
construction . 


Note:  drawing  not  to  scale. 


Figure  1-1.  Concept  Design  of  the  Wright-Patterson  Air  Force  Base 
Wetland  Treatment  Cells 

The  soil  composition  of  the  subsurface  media  differs 
slightly  between  cells.  In  one  cell,  18"  of  historically 
saturated  or  hydric  soil  (i.e.,  soil  characterized  by 
anaerobic  or  reducing  conditions  when  saturated)  with  10% 
woodchip  amendment  sits  over  the  bottom  gravel  layer 
followed  by  two  18"  layers  of  hydric  soil  without  any  added 
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woodchips .  The  woodchips  in  the  bottom  layer  were  added  to 
provide  a  more  concentrated  source  of  carbon  for  the 
anaerobic  microbes.  In  the  second  cell,  an  18"  layer  of 
locally  obtained,  iron-rich  fill  is  sandwiched  between  two 
18"  layers  of  hydric  soil  (without  any  woodchips  added) . 

The  iron-rich  layer  of  indigenous  soil  was  added  to  observe 
the  effects  of  Fe^^  on  the  reduction  of  vinyl  chloride  -  a 
degradation  product  of  PCE  and  potent  human  carcinogen. 
Vinyl  chloride  is  known  to  readily  degrade  under  iron 
reducing  conditions  (Bradley  and  Chapelle,  1997) .  The 
differences  in  cell  design  will  allow  researchers  to 
compare  the  efficiency  of  the  different  types  of  media  and 
recommend  an  optimal  design  approach.  Both  cells  are 
planted  with  various  plots  of  emergent  wetland  vegetation. 
Table  1-1  provides  a  summary  of  the  subsurface  sediment 
composition . 


Table  1-1.  Composition  of  the  Soil  Layers 


Laver 

Cell  1 

Cell  2 

Top 

Hydric  Soil  (likely  root  zone) 

Hydric  Soil  (likely  root  zone) 

Middle 

Hydric  Soil 

Local,  Iron-rich  Fill 

Bottom 

Hydric  Soil  (organic  matter  added) 

Hydric  Soil 

Since  their  construction,  researchers  active  in  the 
study  have  voiced  concerns  that  water  may  be  preferentially 
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flowing  through  certain  portions  of  the  wetland  media.  Of 
equal  concern  is  whether  this  particular  design  approach 
(i.e.,  water  distributed  through  pipes  in  a  gravel  layer 
with  uniform  sediment  layers  above  and  an  outfall  at  one 
end  of  the  cell)  facilitates  uniform  vertical  flows  and 
whether  soil  heterogeneities  -  or  pockets  of  low  hydraulic 
conductivity  -  disrupt  that  flow.  This  research  seeks  to 
characterize  the  behavior  of  water  within  the  wetland  cells 
and  address  such  concerns.  An  understanding  of  how  water 
moves  through  the  wetland  media,  and  how  long  it  resides  in 
the  wetland  media,  is  critical  to  understanding  what 
processes  may  be  at  work  to  degrade  the  PCE .  The  results 
of  this  study  will  add  to  efforts  that  seek  to  explore  the 
potential  of  wetlands  to  de-chlorinate  PCE  and  other 
related  substances.  The  ultimate  goal  of  this  and  related 
research  is  to  develop  design  guidelines  should  such  a 
remediation  approach  prove  relevant  and  successful. 

Problem  Statement 

Early  indications  are  that  the  Wright-Patterson  Air 
Eorce  Base  wetland  cells  are  effectively  de-chlorinating 
PCE  from  the  contaminated  groundwater.  Little  is  known, 
however,  about  the  specific  processes  that  are  affecting 
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the  contaminant's  removal,  where  such  processes  are 
occurring,  and  at  what  rate  they  occur.  A  better 
understanding  of  how  water  flows  through  the  wetland  media 
will  facilitate  efforts  that  seek  to  answer  such  questions. 

Research  Questions 

(1)  Is  the  flow  of  water  uniform  through  the  wetland 
sediments  or  has  the  water  found  preferential  flow  paths? 

(2)  Does  the  behavior  of  groundwater  flow  change  with 
varying  loading  rates  or  environmental  conditions? 

(3)  What  is  the  approximate  residence  time  for  groundwater 
particles  moving  through  the  subsurface  media? 

Scope  and  Limitations 

This  research  will  characterize  the  flow  of 
groundwater  through  the  wetland  media  and  give  subsequent 
researchers  an  idea  of  where  to  focus  their  investigative 
efforts.  Observations  of  hydraulic  head  will  provide  an 
indication  of  flow  behavior  in  the  wetland.  Analysis  of 
various  wetland  soil  parameters  will  allow  simulation 
modeling  and  provide  a  visual  representation  of  the  water 
dynamics  within  the  media.  Subsequent  runs  of  the  model 
will  aim  to  fit  the  observed  hydraulic  head  observations 
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with  the  numerical  model's  calculations.  The  fit  data  will 


enable  calculations  of  hydraulic  residence  times  and  may 
contribute  to  the  development  of  design  parameters  that 
apply  to  like  systems. 
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II .  Literature  Review 


Wetlands 

"Wetlands  are  defined  as  land  where  the  water  table  is 
at  (or  above)  the  ground  surface  long  enough  each  year  to 
maintain  saturated  soil  conditions  and  the  growth  of 
related  vegetation"  (Reed  et  al . ,  1995).  Wetlands  can 
exist  as  transitional  areas  between  terrestrial  and  aquatic 
ecosystems,  and  they  can  also  occur  as  isolated  ecosystems 
wherever  ground  waters  intersect  the  earth's  surface. 
Wetlands  can  derive  their  water  from  surface  sources  -  such 
as  precipitation,  runoff  from  rainfall  events,  and  flooding 
or  overflows  from  adjacent  bodies  of  water  -  or  they  can  be 
groundwater  fed. 

Wetlands  are  complex  ecosystems  that  perform  a  variety 
of  beneficial  functions  in  nature.  As  depressions  in  the 
ground,  wetlands  can  buffer  downstream  locations  from  the 
effects  of  heavy  rainfall  events  and  reduce  the  potential 
for  flooding  (Wetlands  Research  Program,  1993)  .  When 
situated  next  to  other  bodies  of  water,  they  provide 
erosion  control  from  rainfall  and  tidal  flooding.  The 
damping  motion  caused  by  wetland  vegetation  also  reduces 
the  velocity  of  waters  passing  through  the  wetland;  this 
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discourages  the  re-suspension  of  sediments,  which  in  turn 
improves  water  clarity  and  quality  (Wetlands  Research 
Program,  1993).  As  highly  productive  ecosystems,  wetlands 
can  support  life  at  several  ecological  levels  and  provide  a 
suitable  habitat  for  migratory  water  foul  and  many  other 
kinds  of  species.  Their  service  as  nature's  "kidneys  of 
the  landscape"  (Mitsch  and  Gooselink,  1993)  makes  them 
particularly  beneficial  in  areas  where  they  recharge 
surface  and  ground  waters.  An  appreciation  of  this 
capability  has  made  wetlands  a  popular  alternative  to 
remedy  some  of  the  more  concentrated  streams  of  pollution 
generated  by  man. 

There  are  essentially  two  broad  classifications  of 
wetlands:  natural  and  manmade.  While  the  different  sub¬ 
classifications  of  natural  wetlands  may  be  of  interest  in 
ecological  studies,  a  discussion  on  natural  wetlands  here 
is  relevant  only  where  the  mechanisms  of  pollutant  removal 
apply  to  manmade  wetlands  as  well.  Because  natural 
wetlands  are  considered  waters  of  the  United  States,  they 
are  subject  to  regulatory  control  under  the  Clean  Water  Act. 
The  Clean  Water  Act  requires  a  permit  for  the  addition  or 
discharge  of  any  pollutant  from  a  point  source  into  waters 
of  the  United  States,  and  such  permits  normally  require 
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pre-treatment  to  specified  effluent  standards  (Gallagher, 
1999) .  This  regulatory  status  generally  limits  their  use 
in  pollutant  remediation.  Furthermore,  because  uniform 
hydraulic  flows  seldom  occur  in  natural  wetlands,  only  a 
small  portion  of  the  wetland  may  actually  come  into  contact 
with  the  waterborne  pollutants  that  enter  the  system.  It 
is  often  not  possible  to  correct  such  treatment 
inefficiencies  barring  a  significant  engineering  effort, 
and  this  -  along  with  the  addition  of  any  contaminants  - 
would  undoubtedly  alter  the  ecology  of  these  naturally 
developed  systems  (Reed  et  al . ,  1995). 

Manmade  wetlands  constructed  to  treat  pollutants,  on 
the  other  hand,  are  largely  free  from  the  regulatory 
oversight  and  ecological  concerns  inherent  with  natural 
aquatic  ecosystems  (Reed  et  al . ,  1995).  Furthermore,  the 
construction  of  wetlands  to  treat  pollutants  allows  for  an 
optimum  design  consistent  with  known  parameters.  While 
there  are  always  difficulties  in  engineering  efforts, 
constructed  wetlands  can  often  be  designed  to  accommodate 
flows  that  achieve  a  desired  level  of  contaminant  treatment. 

There  are  two  main  types  of  manmade  wetlands:  free 
water  surface  and  subsurface  flow  wetlands.  The  free  water 
surface  variety  is  characterized  by  a  water  surface  that  is 
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exposed  to  the  atmosphere.  Free  water  surface  wetlands  may 


contain  floating,  submerged,  and  emergent  vegetation.  The 
concept  behind  free  water  surface  wetlands  is  that 
pollutants  in  the  water  flow  across  the  surface  of  the 
wetland  and  degrade  through  various  natural  processes 
active  in  the  largely  aerobic  (i.e.  plant)  portion  of  the 
wetland.  The  subsurface  flow  wetland,  in  contrast,  is 
characterized  by  a  water  table  that  is  below  the  surface  of 
the  wetland.  Water  flows  through  a  subsurface  flow  wetland 
in  a  similar  fashion  as  groundwater  flows  through  an 
aquifer.  The  concept  is  that  the  pollutants  in  the  water 
flow  uniformly  through  the  wetland  media  (rather  than  over 
the  surface  of  it)  and  degrade  through  processes  that  occur 
in  this  largely  anaerobic  region.  A  sufficiently  porous 
media  is  a  pre-requisite  for  this  type  of  wetland  to 
function  properly.  Wetland  plants,  if  present,  may  provide 
for  a  slightly  aerobic  region  in  the  root  zone  strata  of 
subsurface  flow  wetlands. 

The  subsurface  flow  wetland  is  generally  considered 
the  more  efficient  of  the  two  types  of  constructed  wetlands 
(Reed  et  al,  1995) .  This  may  be  somewhat  counterintuitive 
to  persons  familiar  with  wastewater  treatment  processes  as 
rates  of  aerobic  degradation  generally  far  exceed  those  of 
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anaerobic  processes.  In  manmade  wetlands,  however,  the 
degradation  processes  believed  to  be  at  work  occur  largely 
as  the  result  of  attached  growth  organisms,  and  the  larger 
reactor  volume  of  subsurface  flow  wetlands  make  up  for  the 
(generally  slower)  anaerobic  processes  at  work  in  them. 

The  wetland  cells  at  Wright-Patterson  Air  Force  Base 
represent  a  hybrid  design  where  water  moves  vertically  up 
though  the  wetland  media  until  it  pools  and  drains  off  the 
surface  of  the  wetland.  While  the  water  is  moving  up 
through  the  wetland  media,  the  cells  function  as  subsurface 
flow  wetlands.  As  the  water  moves  across  the  surface  of 
the  wetland,  the  cells  function  as  free  water  surface 
wetlands . 

Contaminant  Fate  in  Wetlands 

Microorganisms,  plants  and  wetland  sediments  may  all 
play  an  active  role  in  the  pollutant  degradation  processes 
at  work  in  both  natural  and  manmade  wetlands.  The 
temperature,  pH,  dissolved  oxygen  concentration,  and  water 
depth  of  the  wetland  influence  these  components  and  affect 
the  overall  system  performance  (Reed  et  al,  1995)  . 

Microorganisms .  Bacteria  are  the  unchallenged 
champions  of  pollutant  degradation  in  wetland  ecosystems  as 


2-5 


they  are  in  most  wastewater  treatment  systems.  Bacteria 
manufacture  enzymes  that  enable  them  to  consume  or  modify 
pollutants  in  the  environment.  In  most  cases,  the  presence 
of  the  right  kind  of  bacteria  will  result  in  some  degree  of 
pollutant  degradation,  while  their  absence  may  mean  that 
the  wetland  is  only  capable  of  facilitating  transport  to 
destinations  unknown.  Conceptually,  aerobic  bacteria  are 
the  predominant  consumers  of  nutrients  in  free  water 
surface  wetlands,  while  anaerobic  bacteria  inhabit  the 
sediments  of  subsurface  flow  wetlands. 

Bacteria  are  not  the  only  microorganisms  that  take  up 
their  residence  in  wetlands.  The  presence  of  wastes, 
bacteria  and  decaying  matter  will  attract  other  consumers 
(primary,  secondary  and  tertiary)  of  the  animal  world.  For 
a  contaminant  such  as  PCE  however,  the  right  kind  of 
bacteria  are  believed  to  be  the  sole  consumers  or  modifiers 
of  the  chlorinated  solvent.  The  presence  of  higher  animals 
-  sometimes  waterfowl,  fish  and  even  small  mammals  - 
signals  a  vibrant  ecosystem  with  several  trophic  levels 
present  above  that  occupied  by  the  active  microorganisms. 

Plants.  Plants  play  a  beneficial  role  in  wetland 
ecosystems,  although  there  is  disagreement  over  the  extent 
to  which  they  contribute  to  pollutant  removal.  Some 


2-6 


possible  contributions  may  include  plant  uptake,  support  of 
microorganisms,  and  modification  of  soil  properties. 

During  transpiration,  plants  draw  in  water  through  the 
roots  to  support  various  plant  processes.  Plants 
internally  transport  this  water  to  the  leaves  to  support 
photosynthesis,  and  they  transpire  the  excess  water 
(actually  most  of  what  they  take  in  through  the  roots)  to 
the  atmosphere  through  openings  in  the  leaves  called  stoma. 
Research  indicates  that  plants  may  also  be  able  to  expel 
volatile  waterborne  contaminants  along  with  the  water  they 
transpire  (Newman  et  al . ,  1997;  Nietch  et  al .  ,  1999). 
Another  possible  fate  for  waterborne  contaminants  is 
storage  in  the  internal  plant  tissues  by  diffusion. 
Contaminants  warehoused  in  this  manner  may  or  may  not 
exceed  the  plant's  tolerance  and  unleash  a  toxic  effect  on 
the  plant.  Additionally,  stored  contaminants  are  often 
released  back  into  the  wetland  when  the  plants  die  and 
decompose . 

During  transpiration,  plants  also  take  in  carbon 
dioxide  and  oxygen  through  their  stoma  in  the  leaves  to 
support  photosynthesis  and  respiration.  Plants  transport 
some  of  the  oxygen  generated  from  these  processes  to  the 
roots  along  with  the  products  of  photosynthesis  to  maintain 
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cell  functions.  Various  studies  indicate  that  some  of  the 


oxygen  transported  to  the  roots  diffuses  into  the 
surrounding  soil  creating  aerobic  microzones  in  the  root 
zone  area  of  the  wetland  sediments  (Brix,  1994  and 
referenced  sources) .  While  the  amount  of  oxygen  delivered 
in  this  manner  may  be  minimal,  it  could  provide  attached 
aerobic  microorganisms  with  enough  oxygen  to  transform 
contaminants  into  substances  more  readily  taken  up  by  the 
plant,  assimilated  into  its  tissues,  or  volatilized  into 
the  air  from  stoma  in  the  leaf.  The  submerged  portions  of 
wetland  plants  also  provide  a  surface  for  microorganisms  in 
free  water  surface  wetlands  to  attach  themselves  to  (Reed 
et  al . ,  1995) .  In  addition  to  these  potential 
contributions,  plants  may  also  alter  the  physical 
properties  of  the  wetland  soil  as  will  be  discussed  later 
in  the  chapter. 

Soil.  In  their  saturated  condition,  wetland  soils  are 
a  mixture  of  mineral  sediments,  organic  matter  and  water. 
The  water  occupies  -  and  in  the  case  of  subsurface  flow 
wetlands  moves  through  -  the  voids  or  pore  spaces  that 
exist  between  the  individual  soil  particles.  Various 
physical  properties  influence  the  manner  and  rate  at  which 
water  moves  through  these  voids,  and  these  will  be 
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introduced  later  in  the  chapter. 


For  now  it  is  sufficient 


to  note  that  water  will  move  through  different  soils  at 
different  rates,  and,  generally,  the  longer  the  water 
remains  in  the  soil,  the  longer  contaminants  have  to  react 
in  the  wetland  sediment. 

Soils  also  differ  in  chemical  composition,  and  this 
can  affect  the  fate  of  some  contaminants  in  the  water.  For 
example,  soils  that  are  high  in  organic  content  (such  as 
many  clays  and  silts)  can  react  readily  with  various 
waterborne  pollutants,  while  soils  that  are  low  in  organic 
content  (such  as  gravels  and  sands)  are  often  relatively 
inert.  Metals  and  many  organic  contaminants  often  exhibit 
a  strong  affinity  for  soils  that  are  high  in  organic 
content  and  may  actually  sorb  -  or  attach  -  to  the 
individual  soil  particles.  Sorption  is  not  necessarily  a 
permanent  phenomenon,  and  often  contaminants  that  have 
sorbed  to  soil  particles  will  de-sorb  back  into  the  water 
due  to  a  change  in  environmental  conditions.  A  shift  in 
the  water's  pH,  for  example,  can  generate  such  a  change  in 
affinity . 

Like  the  surfaces  and  roots  of  wetland  plants,  the 
soil  particles  in  a  wetland  also  provide  a  surface  for  the 
microorganisms  active  in  pollutant  degradation  to  attach 
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themselves  to.  Finally,  wetland  soils  filter  suspended 
solids  out  of  passing  waters,  and  they  contain  and  store 
the  nutrients  that  are  generated  from  decomposing  plant  and 
animal  matter  in  the  wetland;  such  nutrients  are  essential 
to  support  new  plant  and  animal  growth. 

Organic  Contaminants 

While  there  are  a  number  of  processes  that  can  affect 
contaminant  fate  in  the  environment,  organic  contaminants 
are  particularly  susceptible  to  the  transformation, 
volatilization,  and  sorption  pathways  (Sawyer  et  al .  ,  1994). 

Transformation  reactions  can  occur  through  biotic  and 
abiotic  mechanisms.  Bacteria  can  transform  contaminants 
through  microbial  processes,  and  this  is  frequently  the 
more  rapid  of  the  transformation  mechanisms. 

Transformation  reactions  relevant  to  organic  contaminants 
in  the  environment  include  oxidation,  reduction,  hydrolysis 
and  photolysis.  Of  these,  oxidation,  reduction,  hydrolysis 
reactions  may  occur  through  biotic  and  abiotic  mechanisms 
(Sawyer  et  al . ,  1994). 

Volatilization  is  the  process  that  occurs  when  a 
contaminant  vaporizes  upon  being  exposed  to  the  atmosphere. 
PCE,  as  well  as  all  of  its  chlorinated  degradation  products. 
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is  highly  volatile,  although  like  other  organic  compounds 
it  is  generally  less  volatile  at  lower  temperatures  (Sawyer 
et  al . ,  1994). 

Sorption  was  discussed  earlier  in  the  chapter.  It  is 
a  process  by  which  a  contaminant  sorbs  or  attaches  to  a 
solid  particle's  surface.  Sorption  is  a  general  term  that 
is  often  used  to  describe  natural  processes  (rather  than 
engineered  processes  where  the  terms  adsorption  and 
absorption  carry  more  precise  meanings) .  Hydrophobic 
organics  like  PCE  can  be  particularly  susceptible  to  such  a 
fate  in  wetlands  due  to  the  high  organic  content  of  many 
hydric  soils  (Sawyer  et  al . ,  1994;  Charbeneau,  2000).  As 
indicated  earlier,  sorption  is  not  necessarily  a  permanent 
phenomenon:  contaminants  can  sorb  to  soil  particles  and  de¬ 
sorb  back  into  the  water  as  a  result  of  changing 
environmental  conditions.  A  contaminant  sorbed  to  a  soil 
particle's  surface,  however,  may  be  more  available  to  the 
microorganisms  that  are  also  resident  there. 

PCE  Degradation 

The  transformation  of  PCE  to  un-chlorinated  substances 
is  believed  to  occur  through  different  biologically 
mediated  mechanisms,  each  requiring  the  right  environmental 
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conditions  at  the  right  stage  in  the  process.  A  basic 
understanding  of  how  PCE  is  believed  to  degrade  in  the 
environment  will  set  the  stage  for  why  it  is  important  to 
understand  the  flow  of  water  through  a  wetland  designed  for 
its  remediation. 

PCE  is  a  chlorinated  ethene  that  is  represented  as 
follows : 

Cl  Cl 

\  / 

c=  c 

/  \ 

Cl  Cl 

As  seen  in  this  representation,  PCE  is  a  fully  chlorinated 
ethene.  If  a  molecule  of  PCE  looses  a  chlorine  molecule 
and  gains  a  hydrogen  molecule,  it  becomes  trichloroethylene 
(TCE) .  If  PCE  exchanges  two  chlorine  molecules  for  two 
hydrogen  molecules,  it  becomes  dichloroethylene  (DCE) .  The 
chlorinated  ethene  with  only  one  chlorine  molecule  and 
three  hydrogen  molecules  is  known  as  chloroethylene  or 
vinyl  chloride.  The  de-chlorination  of  PCE  is  known  to 
occur  sequentially,  and  a  conceptual  pathway  for  PCE  de¬ 
chlorination  is  demonstrated  in  figure  2-1. 

Although  similar  in  their  basic  structure,  these 
related  chlorinated  ethenes  degrade  through  different 
processes  and  under  different  environmental  conditions. 
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Cl  H  Cl  Cl  H  Cl 
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/  \  /  \  /  \ 

Cl  H  H  H  Cl  H 


1,1 -DCE 


trans-1 ,2-DCE 

/■ 
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Vinyl  chloride  (VC) 


Bhene 


Figure  2-1.  Proposed  PCE  Dechlorination  Pathway  by  PCE/TCE  De- 
chlorinating  Granules.  (Source:  Fogel  et  al.  “PCE  Treatment  in  Saturated 
Soil  Columns  with  Methanogens”  in  Bioremediation  of  Chiorinated 
Soivents.  Eds.  R.E.  Hinchee,  A.  Leeson,  and  L.  Semprini.  Columbus,  OH: 
Battelle  Press,  1995.) 
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Hoefer  (1999)  summarizes  the  likely  oxidation  and  reduction 
reactions  that  occur  as  the  result  of  microbial  activity. 
Such  reactions  may  produce  energy  for  the  microbes  as  the 
intended  result  of  the  their  enzymatic  activity,  or  they 
may  be  co-metabolic  and  occur  as  the  result  of  enzymatic 
activity  intended  for  other  purposes. 

Of  the  chlorinated  ethenes,  vinyl  chloride  can  readily 
degrade  under  aerobic  conditions  as  an  energy-yielding 
oxidation  reaction.  TCE,  cis-DCE  and  vinyl  chloride  can 
degrade  as  a  result  of  co-metabolic  oxidations,  while  all 
of  the  chlorinated  ethene  variants  can  degrade  under  the 
reducing  conditions  present  in  anaerobic  environments.  The 
initial  dechlorination  of  PCE  is  strictly  limited  to 
anaerobic  conditions,  however,  and  can  only  occur  as  a  co- 
metabolic  reductive  reaction  (Lee  et  al . ,  1998). 

In  addition  to  these  biotic  transformation  pathways, 
the  chlorinated  ethenes  may  also  be  able  to  exit  the 
wetland  system  as  a  result  of  volatilization.  This  may  be 
true  especially  for  vinyl  chloride  which  has  a  particularly 
high  Henry's  constant  value  relative  to  the  other 
chlorinated  ethenes.  (The  Henry's  constant  reflects  the 
degree  to  which  a  substance  will  partition  between  an 
aqueous  and  gaseous  state  (Sawyer  et  al . ,  1994);  the  higher 
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the  Henry's  constant,  the  greater  the  tendency  to 
volatilize  or  partition  to  a  gaseous  state.)  Also,  the 
tendency  for  these  organic  contaminants  to  sorb  to  soil 
particles  is  high  given  their  hydrophobic  nature.  While, 
as  indicated,  sorption  is  not  necessarily  a  permanent 
condition,  contaminants  sorbed  to  soil  particles  in  the 
wetland  sediments  may  be  more  available  to  microorganisms 
that  are  resident  and  active  there. 

Exactly  what  processes  are  taking  place  in  the  Wright- 
Patterson  Air  Force  Base  wetland  treatment  cells,  as  well 
as  where  they  are  taking  place,  is  a  matter  of  ongoing 
research.  The  volatilization,  sorption  and  transformation 
reaction  pathways  may  all  be  relevant.  Transformation 
processes  do  require  time,  however,  and  the  amount  of  time 
that  water  remains  in  or  transits  through  the  wetland 
system  could  determine  the  effectiveness  of  the  processes 
responsible  for  PCE  degradation.  It  is  reasonable  to 
assume  that  if  the  flow  of  groundwater  through  the  wetland 
cells  occurs  too  fast,  the  PCE  will  have  less  of  an 
opportunity  to  react  with  the  anaerobic  bacteria 
responsible  for  cleaving  off  the  first  chlorine  molecule. 
Likewise,  if  the  flow  of  groundwater  occurs  predominately 
through  preferential  flow  paths,  only  a  small  portion  of 
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the  wetland  may  be  active  in  de-chlorinating  the  PCE  and 
related  degradation  products.  Any  non-uniform  flow  of 
water  in  the  treatment  cells  will  reduce  their  efficiency 
as  will  any  short-circuiting  of  the  wetland  media. 

Groundwater  Flow 

Various  physical  properties  influence  the  flow  of 
water  through  soil.  In  wetlands,  the  type  of  soil,  the 
strength  of  hydraulic  gradients  and  pressures,  and  other 
factors  such  as  the  roots  of  plants  will  impact  the  manner 
and  rate  at  which  water  moves  through  the  wetland  sediments. 

Soil  Composition  and  Porosity .  Soils  are  a  porous 
media.  The  observation  of  ocean  waters  seeping  into  sand 
on  a  beach  illustrates  this  fact.  A  representative  volume 
of  soil  consists  of  both  solid  particles  and  void  spaces. 

Both  liquids  (normally  water)  and  gasses  (normally  air)  can 
occupy  the  spaces  that  are  empty  or  void  of  the  solid  soil 
particles.  If  the  voids  in  a  soil  are  100%  full  of  water, 
the  soil  is  called  saturated .  Soils  that  are  not  fully 
saturated  are  labeled  unsaturated .  While  this  distinction 
may  be  obvious,  it  distinguishes  between  applicable 
groundwater  equations. 
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The  property  of  soil  that  characterizes  the  amount  of 
void  space  relative  to  total  soil  volume  is  porosity . 
Porosity,  n,  equals 

n  =  —  (1) 

Vt 

where  Vy  is  volume  of  the  void  spaces  and  Vt  is  the  total 
soil  volume.  The  intrinsic  properties  of  a  soil  define  its 
primary  porosity ,  while  various  other  influences,  such  as 
the  roots  of  plants  or  fractures  in  rocks,  contribute  to 
its  secondary  porosity .  Collectively  these  factors  define 
the  percentage  of  water  a  fully  saturated  volume  of  soil 
can  hold. 

Another  aspect  of  porosity  is  also  important  to  note. 
When  water  moves  through  the  ground,  it  moves  through  some 
but  not  all  of  the  soil  pores.  Some  of  the  pores  may  be 
isolated  and  not  facilitate  water  movement.  The  effective 
porosity ,  rig,  of  a  soil  represents  the  volume  of 
interconnected  voids  in  a  given  volume  of  soil  through 
which  the  water  is  actually  able  to  flow.  A  soil's 
effective  porosity  will  be  less  than  the  overall  porosity 
for  the  soil,  especially  for  fine-grained  sediments  (Kresic, 
1997)  . 
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Darcy's  Law. 


Darcy' s  law  can  be  used  to  analyze  the 


flow  of  water  through  wetland  sediments.  Darcy's  Law 
determines  the  amount  of  water  that  moves  through  a 
representative  volume  of  soil  under  saturated  conditions. 
Darcy' s  law  is  defined  as 

Q^-K—  A  (2) 

dl 

where  Q  is  the  flow  rate  or  discharge  (in  units  of  cubic 
length  per  unit  time  [L^/T] ) ,  K  is  hydraulic  conductivity 
[L/T] ,  dh/dl  is  hydraulic  gradient  [dimensionless],  and  A 
is  area  [L^] .  (The  negative  convention  recognizes  that 
water  flows  from  a  position  of  higher  to  lower  hydraulic 
head.)  Another  usage  of  Darcy's  law  defines  the  specific 
discharge  or  Darcy  velocity  as  q  =  Q/A  in  units  of  [L/T] . 
Whereas  Q  reflects  the  amount  of  water  that  can  move 
through  a  representative  volume  of  soil  per  unit  of  time,  q 
reflects  the  rate  at  which  the  water  moves  across  a  given 
cross-sectional  area. 

The  Darcy  velocity  is  a  measure  of  velocity  on  a  macro 
or  superficial  scale.  In  reality,  groundwater  doesn't  move 
through  the  entire  cross  section  of  a  representative  soil 
volume.  The  Darcy  velocity  neglects  the  existence  of  the 
microscopic  pores  within  a  representative  volume  of  soil 
through  which  water  actually  flows.  Because  the  pores 
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represent  only  a  fraction  of  the  total  soil  volume,  water 
will  move  through  the  pores  at  a  rate  faster  than  that 
reflected  by  the  Darcy  velocity.  Freeze  and  Cherry  (1979) 
define  this  rate  as 

Vl--  (3) 
n 

where  Vj,  is  average  linear  velocity  [L/T]  ,  q  is  Darcy 
velocity  [L/T],  and  n  is  porosity  [unitless] .  Kresic 
(1997)  presents  a  similar  variable  but  substitutes 
effective  porosity  (rig)  for  overall  porosity. 

Although  it  can  be  derived  analytically,  groundwater 
flow  calculations  based  on  Darcy' s  law  are  generally  only 
accurate  to  within  an  order  of  magnitude  precision.  This 
is  due  to  the  inherent  variability  that  exists  in  nature. 
For  example,  even  very  similar  soils  can  be  observed  to 
move  water  at  different  rates.  For  groundwater  flow,  such 
variability  manifests  itself  in  the  hydraulic  conductivity 
parameter  in  Darcy's  law. 

Hydraulic  conductivity .  Hydraulic  conductivity 
represents  the  ease  at  which  water  can  move  through  a 
representative  volume  of  soil  (U.S.  EPA,  1992) .  Hydraulic 
conductivity  depends  on  both  the  properties  of  the  liquid 
moving  through  the  soil  pores  and  on  the  properties  of  the 
soil  that  yield  a  unique  configuration  of  pore  spaces.  A 
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fluid  such  as  motor  oil,  for  example,  will  move  through  a 
porous  media  at  a  slower  rate  than  a  less  viscous  fluid 
such  as  water.  In  a  similar  manner,  a  soil  comprised  of 
rough,  angular  particles  will  exhibit  more  drag  on  passing 
water  molecules  than  a  soil  comprised  primarily  of  smooth, 
spherical  particles.  The  contribution  of  the  liquid  and 
soil  particles  can  be  seen  in 


where 

K  =  hydraulic  conductivity  [L/T] 
p  =  fluid  density  [M/L3] 
p  =  viscosity  [M/(L*T)] 
g  =  gravity  [UT2] 
k  =  intrinsic  permeability  [L2] 

In  equation  (4),  the  grouping  of  variables  inside  the 
parentheses  (pg/p)  represents  the  properties  of  the  fluid 
that  contribute  to  hydraulic  conductivity,  while  a  soil's 
intrinsic  permeability,  k,  reflects  the  influence  of  the 
soil  particles  on  the  ease  at  which  a  fluid  moves  through  a 
porous  media. 

Hydraulic  conductivity  can  vary  with  respect  to 
location  and  direction.  If  a  soil's  hydraulic  conductivity 


2-20 


is  constant  with  respect  to  location,  the  soil  is  said  to 
be  homogeneous  (otherwise  it  is  heterogeneous) .  If  the 
horizontal  and  vertical  components  of  a  soil's  hydraulic 
conductivity  are  constant  with  respect  to  direction,  a  soil 
is  labeled  isotropic  (otherwise  it  is  anisotropic) . 

Homogenous  and  isotropic  soil  conditions  permit 
relatively  easy  calculations  of  groundwater  flow.  Soils 
that  diverge  from  the  homogenous,  isotropic  ideal  are 
commonplace  however.  When  layers  of  different  homogenous 
and  isotropic  soils  overlay  each  other  in  the  same  area  of 
interest,  the  condition  is  termed  layered  heterogeneity . 

To  account  for  the  differing  hydraulic  conductivity  values 
between  soil  layers,  the  formation  can  be  treated  as  a 
single  composite  layer  with  notional  horizontal  and 
vertical  components  of  hydraulic  conductivity.  (Figure  2-2 
illustrates  this  concept.)  The  horizontal  and  vertical 
components  of  hydraulic  conductivity  for  this  composite 
layer  are  determined  by 

(5a) 

i=l 

Kv  -  (5b) 

'tbi/Ki 

i=1 

where  n  is  the  number  of  layers. 
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Figure  2-2.  Layered  Heterogeneity.  (Adapted  from  Freeze  and  Cherry,  1979) 


Anisotropic  soil  conditions  also  complicate  analysis 
of  groundwater  flow.  For  anisotropic  conditions,  the  one¬ 
dimensional  form  of  Darcy's  law  (equation  2)  may  be  written 
in  three  dimensions  as 

ur  dh 

=  -Kx  —  (6a) 

cx 

dh 

Vy^-Ky^  (6b) 

8h 

i/z  =  -Kz  —  (6c) 

oz 

where  Vy  and  Wz  represent  directional  components  of  an 

overall  velocity  vector  Vg,*  K^,  Ky  and  Kz  represent  the 
hydraulic  conductivity  values  in  these  directions;  and 
dh/dx,  dh/dy,  dh/dz  represent  the  partial  derivatives  of 
the  hydraulic  head  h  in  the  x,  y  and  z  directions.  This 
simplification  for  anisotropic  conditions  applies  as  long 
as  the  X,  y  and  z  coordinate  axes  coincide  with  the 
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principal  directions  of  hydraulic  conductivity  (Freeze  and 


Cherry,  1979) . 

Various  factors  can  influence  the  conductivity  of 
wetland  sediments.  The  intrinsic  properties  of  the  soil 
have  a  major  influence;  the  presence  of  clay  particles,  for 
instance,  will  generally  reduce  hydraulic  conductivity. 

The  temperature  of  the  environment  will  affect  the  water' s 
viscosity,  with  higher  temperatures  resulting  in  lower 
viscosities  and  a  faster  moving  liquid.  Compaction  of  the 
wetland  sediments  could  degrade  hydraulic  conductivity  by 
reducing  the  pore  volume  available  for  fluid  flow.  The 
presence  of  waterborne  sediments  and  decomposing  plant  and 
animal  matter  will  add  fines  to  the  wetland  media  and 
reduce  the  hydraulic  conductivity  of  the  soil  through 
sedimentation.  The  roots  of  plants  appear  to  degrade 
hydraulic  conductivity  by  forming  a  very  dense  mass  around 
which  the  water  prefers  to  flow  (Bowmer,  1987;  Fisher, 

1990;  McIntyre  and  Riha,  1991;  Hilton,  1993;  Brix,  1994; 
Kadlec  and  Knight,  1996;  and  Marsteiner,  1997) .  Most  of 
these  studies  concern  the  effects  of  roots  in  sand  or 
gravel  based  constructed  treatment  wetlands,  however,  and 
Kadlec  and  Knight  (1996)  suggest  that  reductions  in  fine¬ 
grained  sediments  may  not  be  significant  due  to  the  already 
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low  hydraulic  conductivities  of  these  soils  compared  to  the 
more  permeable  medias .  The  main  concern  over  potential 
soil  clogging  from  sedimentation  and  root  zone  development 
is  that,  eventually,  hydraulic  conductivity  could  degrade 
to  such  an  extent  where  the  water  in  the  wetland  will  seek 
the  path  of  least  resistance  and  short-circuit  portions  of 
the  subsurface  media. 

Hydraulic  conductivity  can  be  measured  in  the  field 
and  in  the  laboratory.  Fetter  (1994)  discusses  relevant 
tests  but  cautions  that  the  results  of  laboratory  tests  may 
differ  from  those  obtained  under  field  conditions.  Ward 
and  Dorsey  (1995)  indicate  that  the  results  of  laboratory 
tests  yield  lower  hydraulic  conductivity  estimates  than  in- 
situ  tests,  while  Herzog  and  Morse  (1990)  report  that 
laboratory  and  field  tests  on  fine  grained  sediments  can 
vary  by  more  than  two  orders  of  magnitude.  In  any  case, 
the  results  obtained  from  any  hydraulic  conductivity  test 
are  merely  order  of  magnitude  estimates  of  a  soil  property 
that  is  highly  variable  in  nature.  Methods  used  to  measure 
hydraulic  conductivity  in  this  study  are  presented  in  the 
following  chapter. 

Hydraulic  Head  and  Gradient.  A  fundamental  concept  of 
groundwater  flow  is  that  water  moves  from  a  position  of 


2-24 


higher  hydraulic  head  (or  energy)  to  a  position  of  lower 
hydraulic  head.  This  principle  reveals  itself  by  the 
hydraulic  gradient  variable  (dh/dl)  in  equation  (2) . 

Hydraulic  gradient  is  defined  as  the  change  in  the 
potentiometric  surface  divided  by  the  horizontal  distance 
over  which  that  change  is  observed.  The  potentiometric 
surface  is  the  elevation  to  which  water  would  rise  in  an 
observation  well  or  piezometer  that  penetrates  an  aquifer. 
The  hydraulic  gradient,  therefore,  is  merely  the  slope  of 
hydraulic  head  plotted  against  horizontal  distance.  Figure 
2-3  illustrates  this  concept. 
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Figure  2-3.  Hydraulic  Gradient. 
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Hydraulic  head  consists  of  three  components:  pressure 
head,  elevation  head,  and  velocity  head.  Pressure  head  is 
the  energy  available  to  a  particle  of  water  due  to  external 
forces  being  exerted  upon  it.  Pressure  head  is  positive  at 
depths  lower  than  the  water  table  and  negative  above  it. 
(The  water  table  is  assigned  an  atmospheric  pressure  value 
of  0.)  Elevation  head  is  the  potential  energy  available  to 
a  particle  of  water  due  to  its  height  above  some  arbitrary 
datum  point.  Velocity  head  is  the  energy  inherent  in  a 
water  particle  due  to  the  particle's  velocity.  Total 
hydraulic  head  is  represented  as 

hj  ^h^+hp  +  hy  (7) 

where  hr  is  total  hydraulic  head,  hz  is  elevation  head,  hp 
is  pressure  head,  and  hy  is  velocity  head.  Additionally, 

h^=z  (8) 

(9) 


where  z  is  the  elevation  above  the  datum  point  [L]  ,  \|/  is 

the  distance  of  fluid  surface  over  the  point  of  measurement 
[L] ,  V  is  velocity  [L/T] ,  and  g  is  gravity  [M/T^]  .  Due  to 
the  characteristically  slow  velocity  of  groundwater  flows, 
the  velocity  head  element  of  the  total  hydraulic  head 
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equation  is  negligible  relative  to  the  other  two  components. 
Figure  2-4  depicts  the  contribution  of  the  pressure  and 
potential  head  components  to  total  hydraulic  head. 

ground  surface 
water  level 

point  of  measurement 
arbitrary  datum  (z=  0) 

Figure  2-4.  Hydraulic  Head,  Pressure  Head,  and  Elevation  Head  for  a  Field 
Piezometer.  (Source:  Freeze  and  Cherry,  1979).  Note:  the  contribution  of 
velocity  head  (v^/2g)  to  total  hydraulic  head  (h)  is  negligible  compared  to  the 
pressure  head  ( \|/ )  and  elevation  head  (z)  and,  therefore,  is  not  depicted. 


Limits  of  Darcy's  Law.  Darcy's  law  does  not  apply  for 
extremely  slow  flows  -  or  "flows  through  low-permeability 
sediments  under  very  low  gradients"  -  and  it  does  not  apply 
to  flows  that  approach  turbulent  conditions  (Freeze  and 
Cherry,  1979) .  The  low-end  limit  is  difficult  to  quantify, 
but  the  high-end  limit  begins  with  Reynolds  numbers  of  10 
and  higher.  For  the  purposes  of  this  discussion,  the 
Reynolds  number.  Re,  for  groundwater  flow  is  defined  as 
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where  p  is  the  fluid  density  [M/L^]  ,  p  is  the  fluid 
viscosity  [M/ (L*T) ] ,  vis  the  specific  discharge  [L/T] ,  and 
d  is  the  mean  grain  diameter  [L] . 

Groundwater  Flow  Determination 

Hydraulic  head  measurements  can  be  used  to  determine 
the  flow  patterns  of  water  in  the  ground.  Because  water 
flows  from  a  position  of  higher  to  lower  hydraulic  head, 
pieziometric  measurements  can  indicate  the  direction  of 
groundwater  flow  and  permit  the  construction  of  flow  nets. 

Flow  Direction.  Most  groundwater  literature  assumes 
horizontal  flow  (i.e.  the  horizontal  component  of  flow 
greatly  exceeds  that  in  the  vertical  direction  as  depicted 
in  Figure  2-3.)  The  assumption  of  horizontal  flow  is 
generally  true  in  aquifers  -  or  highly  permeable  sediments 
-  relative  to  confining  layers.  Highly  permeable  layers 
act  as  conduits  and  facilitate  horizontal  flows  while 
layers  of  low  permeability  impede  flow  and  encourage  water 
to  move  vertically  along  the  most  direct  route  out  of  the 
media.  Areas  of  recharge  -  where  waters  recharge  the 
ground  below  -  or  discharge  -  where  water  discharges  into 
the  area  -  can  also  have  substantial  vertical  gradients. 
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Natural  springs  are  another  phenomena  where  the  vertical 
component  of  flow  can  be  significant. 

Observation  wells  and  piezometers  sited  at  the  same 
elevation  will  exhibit  head  loss  in  the  direction  of  flow 
provided  the  movement  of  water  is  predominately  horizontal. 

If  a  significant  vertical  component  to  flow  is  suspected, 
the  use  of  nested  piezometers  can  reveal  such  a  phenomenon. 

A  nest  of  piezometers  is  a  collection  of  piezometers  that 
penetrate  a  groundwater  formation  in  very  close  proximity 
to  each  other  but  at  different  depths.  The  existence  of  an 
upward  vertical  gradient  will  manifest  itself  in  hydraulic 
head  readings  that  increase  with  piezometer  depth  (Sprecher, 
2000) .  In  other  words,  the  static  water  level  in  deep 
piezometers  will  exceed  that  in  shallow  piezometers  of  the 
same  nest.  (The  opposite  situation  would  characterize  a 
recharge  zone.)  For  a  wetland  with  a  water  table  above 
ground,  the  existence  of  an  upward  vertical  gradient  will 
manifest  itself  in  piezometer  water  levels  that  are  above 
the  surface  of  the  water  (Haynos,  1991)  . 

Flow  Nets.  A  flow  net  is  a  graphical  representation 
of  water  particle  trajectories.  A  flow  net  consists  of 
representative  particle  flow  paths  (or  flow  lines)  that 
traverse  equipotential  lines  in  the  direction  of  decreasing 
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hydraulic  head.  Equipotential  lines  represent  the  locus  of 
all  points  of  equal  hydraulic  head.  The  combination  of 
equipotential  and  particle  flow  lines  results  in  a  flow  net. 
Figure  2-5  is  an  example  of  a  simple  flow  net. 


Figure  2-5.  Simple  Flow  Net 


Various  useful  references  present  techniques  for  flow 
net  construction  (Freeze  and  Cherry,  1979;  Cedergren,  1989; 
Fetter,  1994;  Kresic,  1997) .  All  of  these  authors  remark 
that  the  construction  of  flow  nets  is  somewhat  of  an  art. 
For  homogeneous  and  isotropic  conditions,  the  procedures 
for  flow  net  construction  are  as  follows: 

(1)  Flow  lines  intersect  equipotential  lines  at  right 
angles  to  form  a  mesh  of  curvilinear  squares. 

(2)  Equal  gradients  of  hydraulic  head  separate 
consecutive  equipotential  lines  as  indicated  in  Figure  2-5. 
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(3)  Equipotential  lines  form  right  angles  with 
impermeable  boundaries  while  flow  lines  run  parallel  to 
them  (there  can  be  no  flow  across  an  impermeable  boundary) . 

(4)  Equipotential  lines  run  parallel  to  constant  head 
boundaries  while  flow  lines  form  right  angles  with  them. 

(5)  Elow  lines  form  oblique  angles  with  water  table 
boundaries  where  discharge  or  recharge  conditions  exist  but 
remain  parallel  to  it  in  the  absence  of  discharge/recharge. 

In  heterogeneous  and  isotropic  systems,  flow  lines 
intersect  equipotential  lines  at  right  angles  within 
formations  of  constant  conductivities,  but  they  refract  at 
boundaries  of  dissimilar  conductivities.  When  crossing  a 
conductivity  boundary,  groundwater  obeys  the  tangent  law 
where 

Ki  _  tan( 01 )  ^^2) 

K2  tan(92) 

as  demonstrated  in  Figure  2-6.  In  Figure  2-6,  the  water  is 
moving  from  a  media  of  lower  to  higher  hydraulic 
conductivity . 

Because  of  the  refraction  of  flow  lines,  curvilinear 
squares  cannot  exist  throughout  the  entire  region  of 
interest  in  heterogeneous  systems.  If  a  flow  net  is  drawn 
with  curvilinear  squares  in  the  region  of  Ki,  for  instance. 
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Figure  2-6.  Refraction  of  a  Flow  Line  Crossing  a  Hydraulic 
Conductivity  Boundary  (Adapted  from  Fetter,  1994). 


rectangles  will  result  in  areas  where  K  #  Ki  (Freeze  and 
Cherry,  1979) . 

Anisotropic  conditions  further  complicate  flow  net 
construction.  In  anisotropic  conditions,  flow  lines  do  not 
traverse  equipotential  lines  at  right  angles.  In 
anisotropic  systems,  the  angle  of  intersection  can  be 
determined  graphically  by  using  the  inverse  hydraulic 
conductivity  ellipse  (Freeze  and  Cherry,  1979) .  With 
reference  to  Figure  2-7,  the  technique  presented  by  Fetter 
(1994)  is  as  follows: 

(1)  Draw  an  inverse  hydraulic  conductivity  ellipse 
with  semi-axes  oriented  along  the  major  axes  of  anisotropy. 

(2)  Draw  an  equipotential  line  (dashed  line)  through 
the  origin  of  the  ellipse. 
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Figure  2-7.  Method  to  Determine  Direction  of  Flow  in  Anisotropic 
Media  (Source:  Fetter,  1994). 


(3)  Draw  an  arrow  (thin  arrow)  in  the  direction  of 
the  hydraulic  gradient  that  emanates  from  the  origin  of  the 
ellipse  and  is  perpendicular  to  the  equipotential  line. 

(4)  Draw  a  tangent  through  the  point  on  the  ellipse 
where  the  hydraulic  gradient  arrow  intersects  the  ellipse. 

(5)  Draw  an  arrow  (bold  arrow)  from  the  origin  of  the 
ellipse  that  intersects  the  tangent  line  at  a  right  angle. 

When  properly  constructed,  flow  nets  can  facilitate 
calculation  of  flow  rates  and  hydraulic  residence  times 
provided  the  hydraulic  conductivity  of  the  media  is  known. 
This  is  possible  by  applying  Darcy' s  law  for  flow  between 
flow  lines  realizing  that  the  quantity  of  water  moving 
between  flow  lines  remains  constant  throughout  the  system. 
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The  sum  of  the  individual  flow  segments,  then,  equals  the 
flow  through  the  entire  system.  This  graphical  method  of 
flow  calculation  is  possible  for  relatively  simple  systems 
that  diverge  little  from  the  homogeneous  and  isotropic 
ideal.  For  highly  complex  systems  where  heterogeneous  and 
anisotropic  conditions  prevail,  such  a  method  is  seldom 
practical  however.  Where  quantitative  analysis  is  not 
possible,  flow  nets  can  still  provide  valuable  qualitative 
insight  into  groundwater  flow  and  behavior.  For  example, 
equipotential  lines  closely  spaced  can  indicate  regions  of 
low  conductivity  where  the  aquifer  thickness  remains 
constant  (Kresic,  1997).  Likewise,  converging  flow  lines 
can  indicate  discharge  areas  while  flow  lines  that  diverge 
can  signal  recharging  conditions  (Fetter,  1994). 

Vertical  gradients  complicate  the  interpretation  of 
potentiometric  contour  maps.  If  a  significant  vertical 
element  to  groundwater  flow  is  suspected,  vertical  contour 
profiles  are  necessary  to  characterize  flow  throughout  the 
three-dimensional  region. 

Various  computer  programs  facilitate  potentiometric 
contouring  of  hydraulic  head.  Surfer®  by  Golden  Software 
is  a  popular  contouring  program  that  permits  quick  display 
of  groundwater  contours.  Surfer  is  capable  of  numerous 
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interpolation  options  (e.g.  triangulation  with  linear 
interpolation  and  geostatistical  methods  such  as  kriging) 
to  fill  in  missing  data  and  produce  contours  of  equal 
hydraulic  head.  The  potentiometric  contour  maps  included 
in  this  study  are  generated  using  Surfer. 

Groundwater  Flow  Analysis 

A  development  of  Darcy' s  law  in  three  dimensions 
provides  the  basic  equations  for  analytical  solution  and 
numerical  modeling  of  groundwater  flow  systems. 

Development  of  Darcy's  Law  in  Three  Dimensions .  The 
one-dimensional  form  of  Darcy' s  law  can  be  applied  to 
three-dimensions  by  considering  the  flow  through  a 
representative  control  volume  (Freeze  and  Cherry,  1979). 

By  applying  a  conservation  of  mass  approach,  the  flow 
through  a  representative  control  volume  is  illustrated  in 
Figure  2-8  where  p  is  the  fluid  density;  v^,  Vy  and  Vz  are 
the  velocity  vectors  of  the  fluid  in  the  x,  y,  and  z 
directions,  respectively;  and  the  operators  djdx.  ,  djdy  ,  and 
djdz  represent  the  change  in  the  product  of  fluid  density 
and  velocity  in  the  x,  y,  and  z  directions,  respectively. 
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pv^  -  a/0z(pv^) 


Figure  2-8.  Three-dimensional  Flow  Through  a  Representative  Control  Volume. 


A  mass  balance  for  steady-state  flow  (i.e.  where  hydraulic 
head  does  not  change  with  time)  yields 

By  assuming  an  incompressible  fluid  -  Freeze  and  Cherry 
(1979)  demonstrate  this  is  a  valid  assumption  for 
groundwater  flow  -  and  by  substituting  the  directional 
equivalent  of  Darcy' s  law  for  Vx,  Wy  and  Vz  vectors 
(equations  6a,  6b,  and  6c,  respectively).  Equation  (13) 
becomes 


A 

dx 


( ^  dh) 

1  ^ 

dh] 

1  ^ 

dll'] 

Kx  — 

1  5xj 

By 

1  ^  5yJ 

dz 

1  dz) 

=  0 


for  an  anisotropic  media.  For  homogenous 
conditions,  =  Ky  =  Kz  and  Equation  (14) 


(14) 

and  isotropic 
reduces  to 
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(15) 


h  h  8^ h 
dx^  dy^  dz^ 

Equation  (15)  is  a  well-known  partial  differential  equation 
known  as  the  Laplace  Equation. 

Freeze  and  Cherry  (1979)  present  a  similar,  but 
slightly  more  rigorous,  derivation  for  transient  flow  in 
anisotropic  conditions  to  obtain 


d(^  dh^  8  8h]  8f^  8h']  S  8h 

8x[^  8x )  8y[  ^  8y  J  8z\  8z)  b  8t 


(16) 


where  S  is  the  storage  coefficient ,  b  is  the  average 
aquifer  thickness  for  confined  systems  or  average  saturated 
thickness  for  unconfined  systems,  and  dh/dt  represents  the 
change  in  hydraulic  head  per  change  in  time.  The  storage 
coefficient  is  a  dimensionless  measure  of  the  volume  of 


water  per  unit  surface  area  stored  or  expelled  from  storage 
for  a  unit  change  in  head  (Fetter,  1994) .  In  confined 
aquifers,  water  is  stored  or  expelled  from  storage  due  to 
compression  or  expansion  of  the  mineral  skeleton  and  liquid 
and 


S  =  bSs  (17) 

where  Ss  is  specific  storage  [L“^].  Specific  storage, 
therefore,  is  a  measure  of  the  amount  of  water  stored  or 
expelled  from  storage  due  to  compression  or  expansion  of 
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the  mineral  skeleton  per  volume  of  media  for  a  unit  change 
in  head  (Fetter,  1994).  Fetter  (1994)  demonstrates  that 
specific  storage  is  represented  by 

Ss^Pw  9((^  +  n P)  (18) 

where 

Pw  =  density  of  water  [M/L^] 
g  =  gravity  [L/T^] 

a  =  compressibility  of  the  aquifer  skeleton  [1/(M/LT^)] 
n  =  porosity  [L^/  L^] 

P  =  compressibility  of  water  [1/M/LT^]] 

Values  of  specific  storage  are  generally  very  small  -  on 
the  order  of  10”^  ft”^  or  less.  For  homogenous  and  isotropic 
conditions.  Equation  (16)  reduces  to 

d^h  d^h  d^h  Ss  ( dll'] 

8x^  8y^  8z^  \dt  J 

for  transient  flow  in  confined  systems. 

In  unconfined  aquifers,  water  is  produced  from  gravity 

drainage  of  the  pore  spaces  as  the  water  table  declines  in 

addition  to  the  water  being  expelled  by  the  compressing 

mineral  skeleton  and  expanding  fluid.  For  unconfined 

systems , 

S^bSs+Sy  (20) 
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where  Sy  is  specific  yield  [dimensionless] .  Specific  yield 
-  or  drainable  porosity  -  is  a  ratio  of  the  volume  of  water 
that  will  drain  under  the  force  of  gravity  from  a  saturated 
volume  of  soil.  Because  the  specific  yield  is  generally 
orders  of  magnitude  larger  than  the  product  of  bSsr  equation 
(20)  is  often  approximated  as 


S  =  Sy  (21) 

for  unconfined  aquifers.  (Fetter  (1994)  cautions,  however, 
that  values  for  specific  yield  may  approach  the  same  order 
of  magnitude  as  the  product  of  bSs  in  fine-grained 
sediments.)  Although  the  saturated  thickness  of  unconfined 
systems  can  change  with  time,  if  the  change  is  small 
compared  to  the  average  aquifer  thickness  then  the  general 
flow  equation  for  unconfined,  transient  conditions  can  be 
approximated  by 


A 

dx 


sh'] 

1  ^ 

( ur  dh] 

1  ^ 

db\ 

Ky  - 

Kv  — 

Ky  - 

1  dx) 

8y 

1  ^  dy] 

8z 

1  dz) 

=  (22: 

b  dt 


in  anisotropic  systems  and 


d^h  d^h  d^h  Sy  f Ob'] 
dx^  8y^  8z^  bKy8t  J 

in  homogeneous  and  isotropic  systems  (Fetter,  1994)  . 

Homogenous  and  isotropic  soil  conditions  permit 
relatively  easy  calculations  of  groundwater  flow.  Under 
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such  conditions,  analytical  solutions  -  which  combine  the 
foregoing  differential  equations  with  relevant  boundary 
conditions  -  are  fairly  straightforward  and  normally  the 
method  of  choice  for  quantitative  analysis.  Soils  that 
diverge  from  the  homogenous,  isotropic  ideal  are 
commonplace  however,  and  aquifers  characterized  by  complex 
heterogeneities  and  anisotropic  conditions  require  the  use 
of  numerical  methods  for  solution. 

Numerical  Models.  Numerical  models  recast  relevant 
differential  equations  (i.e.  Equations  14,  16  or  22)  in 
algebraic  form  and  apply  numerical  approximations  to 
achieve  approximate  solutions  (Fetter,  1994).  A  numerical 
model  solves  one  approximation  for  each  cell  in  a 
subdivided  area  of  interest,  essentially  solving  a  series 
of  X  equations  and  x  unknowns  (where  x  is  the  number  of 
cells  in  the  model) .  The  inherent  complexity  of  numerical 
models  requires  the  use  of  computers  for  solution.  The  two 
most  common  types  of  numerical  modeling  approaches  are  the 
finite  difference  method  and  the  finite  element  method. 

Both  methods  partition  the  area  of  interest  into  a  grid 
that  attempts  to  characterize  the  actual  formation  with 
relevant  hydrological  parameters  and  boundary  conditions. 
The  finite  difference  method  is  simpler  to  conceptualize 
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and  use;  it  divides  the  aquifer  into  a  grid  of  rectangular 
cells  and  is  best  suited  to  simple  aquifer  geometries.  The 
finite  element  method  solves  over  a  grid  of  irregularly 
shaped  triangles  or  quadrilaterals;  although  it  is  more 
difficult  in  practice,  it  can  accommodate  complex 
hydrological  settings  such  as  multiple  principle  directions 
of  anisotropy  (Freeze  and  Cherry,  1979;  U.S.  EPA,  1992) . 

To  build  a  numerical  model,  the  modeler  assigns  known 
parameter  values  (such  as  hydraulic  conductivity  and 
porosity)  and  boundaries  (such  as  no  flow  or  constant  head) 
at  relevant  points  or  cell  nodes  in  the  model.  Once  built, 
the  modeler  attempts  to  simulate  the  hydraulic  head 
distribution  of  an  observed  reference  condition  through 
various  runs  of  the  model;  the  modeler  calibrates  the  model 
by  adjusting  the  input  parameter  values  or  boundary 
conditions  to  correct  deviations  from  the  observed 
conditions.  Once  calibrated,  the  modeler  must  then  verify 
the  model  through  simulations  under  a  different  set  of 
observed  conditions  (such  as  a  flood  event  or  drought  for 
instance) .  A  verified  model  is  one  that  can  accurately 
mimic  many  observed  conditions.  (The  root  mean  square 
error  between  the  calculated  and  observed  values  of 
hydraulic  head  provides  a  quantitative  assessment  of  model 
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accuracy.)  The  strength  of  a  verified  model  is  that  it  can 
be  used  to  explore  the  effect  of  inputs  or  events  yet  to  be 
observed . 

MODFLOW.  MODFLOW  is  a  finite  difference  model  that 
was  developed  by  McDonald  and  Harbaugh  (1988)  for  the  U.S 
Geological  Survey.  MODFLOW  is  one  of  the  more  popular 
numerical  modeling  computer  models,  and  its  accuracy  has 
been  verified  against  numerous  applications  (Pohll,  1993; 
Fetter,  1994;  Kresic,  1997) .  MODFLOW  is  a  highly  flexible 
numerical  modeling  application:  it  contains  a  couple  of 
basic  modules  to  run  the  numerical  simulations  and  a  number 
of  independent  packages  or  sub-routines  to  handle  features 
such  as  constant  sources  of  head,  evapotranspiration  and 
recharge,  wells,  and  drains.  Kresic  (1997)  provides  a 
brief  overview  with  application  and  the  authors  (McDonald 
and  Harbaugh,  1988)  provide  a  more  comprehensive  review. 
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III.  Methodology 


Overview 

A  three-dimensional  grid  of  piezometers  was  installed 
in  wetland  cell  number  one  to  facilitate  the  objectives  of 
this  study  and  the  research  of  students  investigating  PCE 
degradation  in  the  wetland.  Measurements  of  hydraulic  head 
taken  from  water  levels  in  the  piezometer  grid  were  used  to 
determine  the  direction  and  gradients  of  groundwater  flow. 
Contours  of  hydraulic  head  were  produced  from  this  data 
using  an  automated  computer  software  application  to  portray 
the  predominate  behavior  (uniform  or  otherwise)  of  flow. 
Slug  and  pump  test  methodologies  were  employed  on  numerous 
piezometers  in  the  grid  to  obtain  in-situ  estimates  of 
hydraulic  conductivity.  Laboratory  tests  of  soil  samples 
collected  from  various  locations  in  the  wetland  were  used 
to  estimate  porosity  and  effective  porosity.  The  measured 
hydrological  parameters  were  incorporated  into  a  numerical 
computer  model  of  the  wetland  to  further  explore  the 
behavior  and  estimate  the  rate  of  groundwater  flow.  The 
observed  values  of  hydraulic  head  on  a  representative  day 
facilitated  model  calibration,  and  an  alternate  set  of 
boundary  conditions  was  used  to  validate  the  model. 


3-1 


Sampling  Grid 

A  three-dimensional  grid  of  piezometers  was  installed 
in  wetland  cell  number  one  during  the  late  summer  and  early 
fall  months  of  2001.  The  grid  was  designed  to  accommodate 
the  aims  of  this  study  and  the  efforts  of  researchers 
investigating  PCE  degradation  in  the  wetland.  A  systematic 
sampling  grid  was  selected  to  adequately  characterize  the 
flow  of  groundwater  and  facilitate  numerical  modeling; 
although  a  systematic  grid  can  introduce  bias,  it  is  often 
preferred  in  environmental  studies  for  estimating  patterns 
of  contamination  (Gilbert,  1987)  or  in  this  case  trends  in 
groundwater  flow.  Additionally,  Millard  and  Neerchal 
(2001)  point  out  that  the  magnitude  of  such  bias  varies 
inversely  with  grid  resolution:  the  finer  the  grid,  the 
less  the  bias.  The  fine  resolution  of  the  sampling  grid  in 
this  wetland  (relative  to  many  groundwater  investigations) 
alleviates  some  of  the  concern  over  introduced  bias.  Nests 
of  piezometers  were  employed  to  verify  the  existence  and 
magnitude  of  vertical  flow  gradients.  The  nests  were 
spaced  to  ensure  at  least  three  nests  per  plot  of  wetland 
plant  species  to  allow  for  statistical  analysis  of 
contaminant  degradation.  A  total  of  66  nests,  with  three 
piezometers  per  nest,  were  installed  as  shown  in  Figure  3-1. 
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Wetland  Cell  #1 


Figure  3-1.  Piezometer  Placement 
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Piezometer  Installation 


A  total  of  198  Solinst®  Model  615S  Shielded  Drive-Point 
Piezometers  were  assembled,  measured,  and  driven  into  place 
to  construct  the  sampling  grid.  The  stainless  steel 
piezometers  were  chosen  to  minimize  the  potential  of 
reactions  between  the  organic  contaminants  under  study  and 
the  sampling  equipment  (Neislen  and  Schalla,  1991);  an 
interior  W  diameter  Teflon-lined  tube  also  facilitates  a 
high  degree  of  chemical  sampling  integrity.  Most  of  the 
piezometer  tips  were  affixed  to  ’H"  stainless  steel  riser 
pipes;  galvanized  steel  risers,  which  may  be  more 
susceptible  to  corrosion,  were  used  on  13  piezometers  in 
the  top  layer  near  the  weir  end  of  the  wetland  when  the 
supply  of  the  more  expensive  stainless  steel  pipes  ran 
short.  A  grid  of  surveyed-in  stakes  and  string  lines 
ensured  the  proper  placement  of  the  piezometer  nests.  The 
measured  piezometers  were  driven  down  to  a  predetermined 
depth  below  the  surface  of  the  wetland  and  then  retracted 
six  inches  to  separate  the  shield  from  the  screened  portion 
of  the  piezometer  tip.  (The  shield  feature  of  this  model 
piezometer  reduced  the  potential  for  sediment  smearing  and 
clogging  of  the  screens  during  installation.)  Figure  3-2 
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is  a  profile  representation  of  piezometer  placement  in  the 


wetland . 


1 


6"  Sand 

I  //////////////////////////////////////////////////////////// 


Figure  3-2.  Profile  View  of  Piezometer  Placement. 
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Bentonite  pellets  were  placed  around  the  steel  riser 
pipes  to  form  an  effective  seal  between  the  pipes  and  the 
wetland  sediments  following  piezometer  installation.  This 
is  recommended  by  various  sources  to  inhibit  the  vertical 
migration  of  water  and  contaminants  along  the  surface  of 
the  pipes  (Sprecher,  2000;  Nielsen  and  Schalla,  1991), 
although  the  use  of  bentonite  is  not  without  issue.  As  a 
clay  material,  bentonite  has  a  high  cation  exchange 
capability,  and  it  generally  sets  up  with  a  high  pH; 

Nielsen  and  Schalla  (1991)  caution  that  these  two  factors 
could  interfere  with  chemical  monitoring  investigations. 
Additionally,  the  presence  of  a  high  chloride  concentration, 
as  well  as  some  organic  solvents,  may  "interfere  with  the 
ability  of  the  bentonite  to  form  an  effective  seal" 

(Nielsen  and  Schalla,  1991) .  Despite  these  concerns,  the 
use  of  a  bentonite  seal  was  warranted  to  provide  stability 
and  to  keep  many  of  the  piezometers  from  sinking  down  into 
the  wetland.  There  was  an  initial  concern  over  an  apparent 
migration  of  bentonite  from  around  the  riser  pipes  shortly 
after  flooding  of  the  wetland.  Approximately  1^"  of  quick 
setting  cement  was,  subsequently,  poured  around  the  risers 
over  the  bentonite  in  an  effort  to  impede  the  suspected 
migration.  In  hindsight,  a  cap  of  indigenous  wetland  soil 
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as  recommended  by  Sprecher  (2000)  may  have  worked  better 
than  the  Quickcrete®,  as  many  of  the  cement  caps  cracked 
shortly  after  hydration.  (Soil  caps  would  not  crack  upon 
hydration  and,  conceptually,  they  could  contain  the 
bentonite  in  the  annular  space  if  thick  enough.) 
Fortunately,  a  subsequent  eruption  of  bentonite  from  under 
the  cracked  caps  did  not  ensue,  and  additional  migration 
has  not  been  an  issue. 

The  piezometers  were  developed  approximately  one  month 
after  installation  of  the  sampling  grid.  Well  development 
is  a  common  practice  used  to  maximize  the  efficiency  and 
yield  from  production  wells  in  aquifers;  it  often  involves 
the  extraction  or  surging  of  large  amounts  of  water  to 
clear  away  fine  sediments  from  the  well  intake  area.  Its 
practice  in  wetland  studies  is  not  as  well  documented, 
although  Sprecher  (2000)  comments  on  its  use  if  clogging 
has  occurred.  The  requirement  to  develop  the  piezometers 
was  initially  ignored  as  the  researchers  involved  with 
constructing  the  grid  assumed  the  piezometer  shields  would 
prevent  clogging  during  installation. 

The  decision  to  develop  some  of  the  piezometers  was 
based  on  the  requirement  of  the  researchers  investigating 
PCE  degradation  to  be  able  to  purge  and  extract  enough 
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water  for  chemical  sampling;  many  of  the  piezometers  in  the 
top  and  middle  layer  recovered  too  slowly  and  were  not 
yielding  enough  water  to  obtain  representative  samples  for 
analysis.  Twenty-six  piezometers  in  the  top  layer  and  10 
piezometers  in  the  middle  layer  were  initially  developed 
with  this  aim  in  mind.  The  piezometers  were  developed  by 
using  a  Solinst  Model  410  Peristaltic  Pump  in  the  following 
manner:  the  piezometers  were  pumped  dry  and  then  water  from 
the  wetland  was  pumped  back  in,  often  until  bubbles  were 
observed  coming  up  through  the  wetland  sediments;  this 
process  was  repeated  once  again.  A  decision  to  develop  the 
remaining  piezometers  was  based  on  the  results  of  slug 
tests  conducted  before  and  after  the  initial  round  of 
development  that  differed  by  more  than  an  order  of 
magnitude.  The  goal  of  the  second  round  of  development  was 
to  homogenize  the  sampling  grid,  and  the  remaining 
piezometers  in  the  top  and  middle  layers  were  subsequently 
developed . 

A  complete  survey  of  the  sampling  grid  was  performed 
approximately  two  months  after  piezometer  emplacement  using 
Sokkia  Set  2100  total  station  survey  equipment.  The  survey 
was  conducted  by  civil  engineering  personnel  of  the  88^*^ 
Airbase  Wing,  Wright-Patterson  Air  Force  Base  and  is  based 
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on  a  second  order  survey  of  the  base.  The  longitude  and 
latitude  coordinates  were  converted  to  local  coordinates 
for  use  in  the  numerical  computer  model. 

Piezometric  Measurements 

Water  level  measurements  were  taken  on  eight  separate 
occasions  to  observe  trends  in  the  contours  of  hydraulic 
head.  The  measurements  were  taken  using  a  Solinst  Model 
lOlM  Water  Level  Meter  and  sometimes  a  ruler.  Comparison 
of  measurements  taken  with  the  Solinst  water  level  meter 
and  ruler  provided  instrument  familiarization  and  an 
accurate  interpretation  of  the  sound  of  the  meter. 

Each  round  of  water  level  measurements  (i.e.,  198 
piezometers)  took  about  three  hours.  Although  not  a  true 
snapshot  in  time,  the  duration  of  each  sweep  was  assumed 
short  enough  to  ignore  slight  changes  in  barometric 
pressures,  at  least  during  periods  of  clear,  steady  weather. 
(Water  levels  in  small  diameter  piezometers  will  adjust 
rapidly  to  changing  barometric  pressures  and  at  a  much 
faster  rate  than  that  which  the  pressures  are  transmitted 
through  the  wetland  media  (Bouwer,  1978.)) 

The  water  levels  in  the  piezometers  were  recorded  by 
measuring  down  from  the  top  of  the  Teflon-lined  tube.  The 
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distance  from  the  top  of  the  tube  to  the  top  of  the 
surveyed,  above  ground  riser  coupling  permitted  calculation 
of  water  levels  in  the  piezometers  above  the  reference 
datum  (in  this  case  the  approximate  elevation  of  the 
impervious  liner  underlying  the  wetland) .  The  water  levels 
in  the  piezometers  relative  to  this  datum  provide  the 
magnitude  of  hydraulic  head  at  the  midpoint  elevations  of 
the  piezometer  screens  in  the  wetland. 

Displays  of  potentiometric  contour  maps  for  the  three 
soil  layers  were  generated  using  the  Surfer®  groundwater 
contouring  program  by  Golden  Software.  In  order  to 
compensate  for  slight  vertical  deviations  in  piezometer 
screen  intake  elevations,  the  measured  heads  were  adjusted 
prior  to  contouring  to  generate  maps  of  hydraulic  head  at 
the  mean  midpoint  elevation  of  piezometer  screens  for  each 
layer  of  the  wetland.  Least  squares  regression  was  used  to 
establish  a  linear  relationship  of  hydraulic  head  versus 
elevation  for  each  piezometer  nest,  and  the  values  of 
hydraulic  head  were  adjusted  based  on  the  deviation  of  the 
recorded  measurements  above  or  below  the  regression  line. 
The  adjusted  values  were  contoured  using  the  default 
settings  of  the  kriging  interpolation  option  available  in 
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Surfer®  (ordinary  kriging,  linear  variogram,  nugget  effect 
equal  to  zero)  to  interpret  the  data  and  generate  the  maps. 

Parameter  Estimation 

A  substantial  amount  of  time  and  effort  was  spent 
conducting  hydraulic  conductivity  tests  on  the  installed 
grid  of  piezometers.  Sufficient  characterization  of  the 
wetland  was  necessary  to  facilitate  construction  of  the 
numerical  model.  Analysis  using  slug  test  methodology 
(description  to  follow)  was  generally  possible  for  the  top 
and  middle  layer  piezometers  due  to  the  relatively  low 
hydraulic  conductivities  in  this  portion  of  the  wetland. 

Slug  tests  were  not  a  practical  means  to  measure  the 
conductivity  of  the  bottom  layer  sediments,  however,  due  to 
the  nearly  instantaneous  rate  of  water  recovery  in  these 
piezometers.  As  a  result,  pump  tests  were  attempted  on  a 
few  piezometers  to  characterize  this  portion  of  the  wetland. 
Tests  for  porosity  were  also  conducted  and  are  discussed 
herein . 

Slug  Tests.  The  literature  is  open  to  interpretation 
when  selecting  a  relevant  slug  test  for  a  hydrological 
setting  similar  to  that  of  the  Wright-Patterson  Air  Force 
Base  wetland  treatment  cells.  Dawson  and  Istok  (1991) 
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instruct  that  the  Hvorslev  method  applies  only  to  confined 
aquifers  while  the  Bouwer  and  Rice  method  is  appropriate 
for  unconfined  or  leaky  confined  settings.  Although  the 
Wright-Patterson  wetland  cells,  conceptually,  represent  an 
unconfined  hydrological  setting,  the  presence  of  a 
piezometric  surface  above  the  surface  of  the  water  table 
mimics  the  behavior  of  a  confined  aquifer  system. 

Similarly,  while  the  Bouwer  and  Rice  method  measures  the 
rate  of  recovery  in  wells  relative  to  the  initial  water 
table  elevation,  the  Hvorslev  test  measures  the  rate  of 
recovery  relative  to  the  initial  potentiometric  surface. 
Finally,  the  Bouwer  and  Rice  test  was  intended  to  test 
hydraulic  conductivity  from  wells,  while  the  Hvorslev  test 
is  normally  associated  with  hydraulic  conductivity  tests 
using  piezometers.  In  spite  of  these  concerns,  the  Bouwer 
and  Rice  method  was  selected  for  use  in  this  study  because 
the  geometry  of  the  piezometer  intakes  conceptually 
mimicked  that  of  the  Bouwer  and  Rice  test,  and  an 
appropriate  "shape  factor"  for  use  with  the  Hvorslev  method 
(1951)  was  not  found  in  the  literature  until  late  in  study. 
The  observed  rate  of  recovery  in  this  study,  therefore, 
measures  the  hydraulic  head  differential  from  the  initial 
piezometric  surface  rather  than  the  initial  water  table 
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elevation.  Thompson  (1993)  also  used  the  Bouwer  and  Rice 
test  with  piezometers  to  evaluate  hydraulic  conductivity  in 
a  natural  wetland. 

The  Bouwer  and  Rice  slug  test  -  like  the  Hvorslev  test 
-  involves  the  displacement  of  a  known  volume  of  water  in 
an  observation  well  or  piezometer  to  create  a  sudden  change 
in  hydraulic  head;  the  observed  rate  of  water-level 
recovery  following  displacement  is  related  to  Darcy' s  law 
and  can  be  used  to  determine  the  hydraulic  conductivity  of 
the  media  immediately  surrounding  the  intake  area  of  the 
monitoring  device.  Like  the  Hvolslev  test,  the  Bouwer  and 
Rice  test  can  be  used  on  partially  penetrating  and 
partially  screened  wells  and  piezometers  (i.e.,  wells  or 
piezometers  that  do  not  penetrate  and  are  not  perforated 
along  the  entire  saturated  thickness  of  the  aquifer) .  The 
water  for  the  test  can  be  physically  added  to  or  extracted 
from  the  piezometer;  alternatively,  the  water  can  also  be 
displaced  by  introducing  or  extracting  a  solid  metal  object 
to  produce  the  necessary  instantaneous  change  in  head. 

The  well  geometry  for  the  Bouwer  and  Rice  test  is 
illustrated  in  Figure  3-3 
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Figure  3-3.  Geometry  for  Bouwer  and  Rice 
Slug  Test  (Source:  Charbeneau,  2000). 


where 

Lg  =  effective  height  of  the  screened  area  or  perforated  zone 

r„  =  radius  of  the  piezometer  point  or  deveioped  zone 

/■g  =  radius  of  the  standpipe  through  which  the  water  ievei  rises/faiis 

D  =  distance  between  the  equiiibrium  water  ievei  and  bottom  of  the 
piezometer  screen 

H  =  height  of  equiiibrium  water  ievei  above  the  impermeabie  strata 
y  =  water  ievei  foiiowing  dispiacement  from  the  equiiibrium  position 

The  geometry  for  the  piezometers  used  in  this  study  is 
depicted  in  Figure  3-4,  and  the  same  notation  applies. 
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The  basic  equation  to  determine  hydraulic  conductivity 
using  the  Bouwer  and  Rice  test  is 


rc^  ln{Re/r^  )  1 


-  In  ^ 


where  Re  is  the  effective  radial  distance  over  which  the 
head  change  is  dissipated  in  the  media  (Bouwer,  1978) . 
Bouwer  and  Rice  (1976)  present  the  following  equations  to 
compute  the  value  of  In  (Re/r„)  : 
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(25) 


(26) 


Bouwer  and  Rice  (1976)  developed  a  graph  that  relates  the 
constants  A,  B,  and  C  to  the  known  ratio  of  Le/r„.  This 
graph  is  widely  available  in  many  groundwater  texts  that 
discuss  the  test  (e.g.  Bouwer,  1978;  Fetter,  1994;  Kresic, 
1997) .  Bouwer  and  Rice  (1976)  indicate  that  the  variable 
In ( (H-D) /r„)  in  Equation  (25)  has  an  upper  bound  of  6  based 
on  the  results  of  analog  simulations,  so  a  value  of  6 
should  be  used  for  In  ( (H-D) /r„)  if  the  calculated  value 
exceeds  6.  Because  water  flows  into  or  out  of  the 


piezometer  intake  in  a  predominately  radial  -  rather  than 
vertical  -  manner,  the  test  primarily  yields  values  of  Kh 
(Bouwer,  1978) . 

Ideally,  plots  of  the  change  in  hydraulic  head  versus 
time  should  appear  as  indicated  in  Figure  3-5.  Although 
the  variable  yo  is  hard  to  measure  in  a  practical  sense 
without  the  aid  of  automated  measuring  equipment,  its  value 
can  be  closely  computed  using  the  method  of  least  squares 


3-16 


Plot  of  Observations 

a,  14  - 

>. 

DO 

• 

"O 

to 

^  10 

• 

• 

X 

• 

I 

^  6- 

c  ° 
to 

4 

•• 

• 

o  ^ 

••• 

( 

D  5 

Tim 

0  1 

e,  t 

Plot  for  Bouwer  and  Rice 


100 


O)  to 

g  bi 

O  o 

■D  ^ 
a; 

0) 

X 


10 


50 

Time,  t 


100 


Figure  3-5.  Arithmetic  and  Semi-log  Plots  of  Change  in  Head  Versus 
Time  for  the  Bouwer  and  Rice  Slug  Test. 

regression.  Alternatively,  Fetter  (1994)  points  out  that 
any  two  points  along  the  straight  line  portion  of  the  semi¬ 
log  plot  can  be  used  in  Equation  (24)  provided  the  variable 
t  reflects  the  time  interval  between  the  two  points. 
Sometimes,  the  semi-log  plot  of  change  in  head  versus  time 
will  yield  a  curve  with  two  straight-line  portions,  the 
first  having  a  steeper  slope  than  the  second.  Kresic 
(1997)  indicates  this  may  be  due  to  the  more  permeable  sand 
or  gravel  packs  that  surround  many  well  screens  relative  to 
the  undisturbed  surrounding  aquifer  soil.  Values  of  y  and 
t  from  the  second  straight-line  segment  should  be  used  in 
Equation  (24)  to  calculate  hydraulic  conductivity  if  this 
situation  occurs.  Also,  Bouwer  (1989)  reports  that  the 
data  points  will  characteristically  taper  off  and  away  from 
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the  straight-line  portion  of  the  semi-log  plot  during  the 
later  portion  of  the  test. 

Most  texts  recommend  taking  measurements  until  the 
water  has  recovered  to  90%  of  the  pre-test,  equilibrium 
levels  to  obtain  a  representative  amount  of  data  points. 
Kraemer  et  al .  (1990)  suggest  that  a  goal  of  50%  recovery 

may  suffice  in  low  conductivity  sediments  whose  slow 
recovery  rates  do  not  make  a  recovery  of  90%  practical, 
although  they  do  not  provide  a  supportive  basis  for  this. 

Slug  Test  Methodology .  The  tests  were  conducted  by 
adding  or  extracting  water  through  the  interior  W  Teflon- 
lined  tube.  To  measure  the  rate  of  water  level  recovery, 
masking  tape  was  affixed  to  the  outside  of  the  plastic 
piezometer  tube,  and  the  rate  of  recovery  was  marked  at 
regular  intervals  on  the  tape;  the  ability  to  see  the  water 
level  through  the  transparent  tube  greatly  enhanced  the 
degree  of  measuring  precision.  Hydraulic  conductivity 
values  were  computed  using  Equations  (24)  and  (25) .  Values 
for  Yo,  Yt  and  t  were  extracted  from  the  straight-line  fit 
of  data  points  on  the  semi-log  plot  of  y  versus  t.  Where 
the  data  deviated  significantly  from  the  theoretical 
straight  line  fit  (generally  for  an  of  less  than  98%) , 
points  forming  a  straight  line  during  the  middle  portion  of 


3-18 


the  test  were  used  to  compute  hydraulic  conductivity.  A 
total  of  64  and  61  slug  tests  were  conducted  on  the  top  and 
middle  layer  piezometers,  respectively;  some  of  these  were 
repeat  tests  to  assess  the  impact  of  piezometer  development. 

Pump  Tests.  Pump  tests  were  performed  on  select 
piezometers  in  the  bottom  layer  of  the  wetland  to  obtain 
estimates  of  hydraulic  conductivity.  (As  already  noted, 
the  rapid  recovery  rate  of  water  in  the  bottom  layer 
piezometers  precluded  the  use  of  slug  tests  to  characterize 
hydraulic  conductivity  in  this  portion  of  the  wetland.) 

Where  feasible,  pump  tests  are  often  preferred  over  slug 
tests  to  obtain  in-situ  estimates  of  hydraulic  conductivity 
because  they  evaluate  a  much  larger  region  of  soil. 
Additionally,  other  hydrologic  parameters  -  such  as 
specific  storage  and  specific  yield  -  can  also  be  estimated 
depending  on  the  particular  method  used  to  evaluate  data 
from  the  test. 

The  majority  of  the  pump  tests  presented  in  the 
literature  involve  observations  of  water  level  (or 
potentiometric  surface)  decline  in  wells  or  piezometers  as 
another  well  is  pumped,  normally  at  a  constant  rate.  When 
water  is  pumped  out  of  a  well,  a  cone  of  depression  forms 
in  the  water  table  (or  potentiometric  surface)  and  expands 
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outward  from  the  well  until  a  recharge  boundary  is  reached. 
Analytical  equations  that  compute  hydraulic  conductivity 
from  the  drawdown  once  this  steady-state  condition  is 
reached  are  straightforward  but  require  wells  that  fully 
penetrate  the  aquifer.  An  analysis  of  transient  time- 
drawdown  data  is  necessary  for  wells  that  only  partially 
penetrate  an  aquifer  as  well  as  for  tests  that  do  not  reach 
steady  state  conditions. 

Various  methods  for  analyzing  transient  time-drawdown 
data  from  aquifer  pump  tests  exist  in  the  literature.  Many 
employ  a  curve- fitting  technique  that  compares  the 
transient  time-drawdown  data  obtained  from  the  test  with 
the  plot  of  a  theoretical  type-curve .  In  addition  to  the 
extent  of  well  penetration,  the  type  of  aquifer  (confined, 
leaky,  or  unconfined)  also  influences  the  choice  of  an 
appropriate  test  for  transient  conditions.  Dawson  and 
Istok  (1991)  provide  a  useful  summary  of  aquifer  tests  that 
are  applicable  over  a  wide  range  of  geological  settings. 

The  piezometers  installed  in  wetland  cell  number  one 
only  partially  penetrate  the  soil  layer  of  interest,  and  by 
design  the  hydrogeologic  setting  is  conceptually  that  of  an 
unconfined  aquifer.  The  appearance  of  transient  time- 
drawdown  data  collected  during  the  pump  tests,  however. 
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mimicked  that  of  a  confined,  leaky  aquifer  with  influence 
from  a  source  bed  (i.e.,  an  aquifer  with  an  overlying 
confining  layer  that  is  permeable  enough  to  transmit  water 
down  into  the  aquifer  from  a  constant  source  of  recharge) 
as  illustrated  by  Dawson  and  Istok  (1991)  .  Fetter  (1994) 
cautions  that  the  effects  of  partial  penetration  can 
produce  a  time-drawdown  curve  that  appears  similar  to  that 
for  a  confined,  leaky  setting,  but  he  also  considers  an 
aquifer  confined  when  it  is  overlain  by  strata  with  a 
conductivity  that  is  two  orders  of  magnitude  less  than  that 
of  the  aquifer  of  interest.  The  evidence  supporting 
confined  and  leaky  conditions  (rather  than  unconfined)  for 
the  bottom  soil  layer  of  the  wetland  was  two-fold.  First, 
the  water  level  recovery  observed  for  slug  tests  attempted 
in  the  bottom  layer  was  considerably  faster  (virtually 
instantaneous)  compared  to  that  for  tests  conducted  in  the 
top  two  layers  of  soil  indicating  a  hydraulic  conductivity 
much  higher  than  that  above.  Second,  the  pump  tests 
measured  the  decline  of  the  potentiometric  surface  -  rather 
than  a  declining  water  table  -  which  again  is  consistent 
with  a  confined  aquifer  setting. 

Hantush  and  Jacob  (1955)  present  a  method  for 
analyzing  the  time-drawdown  data  obtained  from  pump  tests 
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conducted  in  confined,  leaky  aquifers.  Hantush  (1960) 
adapted  this  method  to  settings  where  the  confining 
aquitard  releases  water  from  storage  during  the  test  and  to 
account  for  partial  penetration  of  the  pumping  and/or 
observation  wells  or  piezometers  used  (Hantush,  1961; 
Hantush  1964) .  Boulton  and  Streltsova  (1975)  present  a 
method  that  accommodates  all  of  these  conditions,  but  it 
uses  equations  that  are  considerably  more  complex  than 
those  presented  by  Hantush. 

The  governing  equations  for  drawdown  at  an  observation 
well  that  fully  penetrates  a  confined,  leaky  aquifer  are 

s  =  -^W{u,rlB)  (27) 

4  kT 


where 


W{u,rlB)= 


u  =  - 


S 

4Tt 


dy 


(28) 


(29) 


(30) 


s  =  drawdown  [L] 

W(u,r/B)  =  well  function  for  confined,  leaky  aquifers 
u  =  dimensionless  parameter 

B  =  leakage  from  the  overlying  aquitard  [dimensionless] 
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Q  =  pump  rate  [L^/T] 

T  =  aquifer  transmissivity  [L^/T] 

S  =  aquifer  storativity  [unitiess] 

r  =  radius  from  the  pumping  weii  to  the  observation  weii  [L] 
b’  =  aquitard  thickness  [L] 

t  =  time  eiapsed  from  the  beginning  of  the  aquifer  test  [T] 
y  =  variabie  of  integration  [dimensioniess] 

Transmissivity  is  defined  as 

T^Kb  (31) 

where  K  is  the  hydraulic  conductivity  [L/T]  and  b  is  the 
aquifer  thickness  [L] .  Aquifer  storativity  is  the  same  as 
the  storage  coefficient  defined  by  Equation  (17)  in  chapter 

2. 


The  equation  governing  drawdown  at  an  observation  well 
that  fully  penetrates  a  confined,  leaky  aquifer  when  the 
overlying  aquitard  releases  water  from  storage  is 

s^[QI4  7rT)H{u,l3)  (32) 

where  the  variables  s,  Q  and  T  are  previously  defined  and 
H(u,/3)  is  the  well  function  that  accounts  for  aquitard 
storage  in  leaky  aquifers  [dimensionless] .  The  leaky  well 
function,  H(u,(5),  is  defined  as 


H(»v)=r 

•'t/ 


-erfc 


y 


p4u 
y(y  -  u) 


dy 


(33) 
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where  the  variables  u  and  y  were  previously  defined,  is  a 
dimensionless  parameter  (defined  below),  and  erfc(x)  is  the 
complementary  error  function  (also  defined  below) .  The 
complementary  error  function,  erfc(x),  and  parameter  are 
defined  by 

erfc(x)^^=\^  e~^  dy  (34) 


K'  S 
TSb' 


11/ 2 


(35) 


where 

T  =  transmissivity  [L^/T] 

S  =  aquifer  storativity  [dimensioniess] 

S’  =  storativity  of  the  overtying  aquitard  [dimensioniess] 

K’  =  verticai  hydrauiic  conductivity  of  the  overtying  aquitard  [UT] 
r  =  radius  from  the  pumping  weii  to  the  observation  weii  [L] 
b’  =  aquitard  thickness  [L] 
y  =  variabie  of  integration 

Equations  (33)  through  (35)  only  apply  for  t<  (b' S' )  /  (lOK' ) 
Hantush  (1964).  Hantush  (1960)  presents  other  equations 
for  analyzing  aquifer  tests  that  extend  beyond  this  time 
window . 

The  equations  governing  drawdown  around  partially 
screened,  partially  penetrating  wells  and  piezometers  are 
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^ave 


^_Q_' 

4nKrb 


[W(Ur,  rIBr)  +  F(Ur,  X,  rjBr  ,  l/b ,  d/b  ,  lo/b  ,  da/b)] 


/•OO 

1 

\  I 

—  exp 

JU 

y 

L  yj 

dy 


S 

4Krbt 


K^bb' 

K' 


(36) 

(37) 

(38) 

(39) 


F(Ur,  X,  rjBr ,  l/b ,  d/b ,  lo/b ,  do/b)  -- 


2b^ 


nUl-d)(lo-do) 


where 

Save  =  average  drawdown  over  the  screened  portion  of  the  well  [L] 
W(ur,r/Br)  =  well  function  for  confined,  leaky  aquifers  [dimensionless] 
ur  =  dimensionless  parameter 

Br  =  leakage  from  the  overlying  aquitard  [dimensionless] 

F(ur,x,r/Br, l/b, d/b, lo/b, do/b)  =  dimensionless  parameter  adjusting  for  partial 

penetration 


Q  =  pump  rate  [L3/T] 

S  =  aquifer  storativlty  [unitless] 

Kr  =  aquifer  hydraulic  conductivity  In  the  radial  direction  [L/T] 
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K’  =  vertical  hydraulic  conductivity  of  the  overlying  aquitard  [UT] 
r  =  radius  from  the  pumping  well  to  the  observation  well  [L] 
b  =  aquifer  thickness  [L] 
b’  =  aquitard  thickness  [L] 

I  =  vertical  distance  from  the  bottom  of  the  pumping  well  screen  to  the 
bottom  of  the  confining  layer  [L] 

lo  =  vertical  distance  from  the  bottom  of  observation  well  screen  to  the 
bottom  of  the  confining  layer  [L] 

d  =  vertical  distance  from  the  top  of  the  pumping  well  screen  to  the  bottom 
of  the  confining  layer  [L] 

do  =  vertical  distance  from  the  top  of  the  observation  well  screen  to  the 
bottom  of  the  confining  layer  [L] 

t  =  time  elapsed  from  the  beginning  of  the  aquifer  test  [T] 

y  =  variable  of  integration  [dimensionless] 

n  =  range  variable  for  the  summation  operator  [dimensionless] 

Equations  (36)  through  (40)  assume  that  the  amount  of  water 
released  from  storage  in  the  overlying  aquitard  during  the 
test  is  negligible  (Hantush,  1961).  Hantush  (1964)  points 
out  that  for  r>1 .5  b-^K^  jK^  (where  r,  b,  and  Kr  are  previously 
defined  and  is  the  hydraulic  conductivity  of  the  aquifer 
in  the  vertical  direction) ,  the  effects  of  partial 
penetration  are  negligible  and  Equations  (27)  through  (31) 
may  be  used  in  the  place  of  Equations  (36)  through  (40) . 
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To  determine  hydraulic  conductivity  using  the  time- 
drawdown  data  collected  from  a  pump  test  in  a  confined, 
leaky  aquifer,  the  following  procedure  applies. 

(1)  Prepare  a  logarithmic  graph  of  type-curves  by 
plotting  the  applicable  well  function  on  the  vertical  axis 
and  u  (or  Ur)  ranging  from  10~®  to  10  on  the  horizontal  axis 
for  various  values  of  r/B,  (3  ,  or  r/Br  (depending  on  the 
particular  method  used) . 

(2)  Prepare  a  logarithmic  plot  of  (observed)  drawdown 
on  the  vertical  axis  and  time  on  the  horizontal  axis  using 
same  scale  as  that  used  for  the  type-curves. 

(3)  Overlay  one  of  the  plots  against  the  other  keeping 
the  axes  aligned  until  a  close  match  between  the  observed 
data  points  and  one  of  the  type-curves  is  achieved. 

Annotate  a  single  match  point  on  each  graph  that  reflects 
this  best  fit,  and  record  the  values  of  the  well  function, 

u  (or  Ur),  drawdown,  and  time  for  the  match-point. 

(4)  Use  the  recorded  match  point  values  to  determine 

transmisivity  from  Equation  (27),  (32),  or  (36).  Use  the 

computed  value  of  transmisivity  to  determine  hydraulic 
conductivity  from  Equation  (31),  and  use  Equation  (29)  or 
(38)  to  determine  storativity  if  required.  With  Equations 
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(30)  and  (39) ,  estimates  of  hydraulic  conductivity  for  the 
overlying  aquitard  are  also  possible. 

Pump  Test  Methodology .  The  pump  tests  were  conducted 
by  draining  water  from  one  piezometer  (the  pumping  well) 
while  observing  the  drawdown  of  water  in  another  piezometer 
(the  observation  well) .  The  aboveground  riser  pipes  were 
removed  from  both  piezometers  and  masking  tape  was  affixed 
to  the  clear  plastic  tubing  of  the  observation  piezometer 
to  record  the  position  of  the  water  level  over  time.  To 
initiate  the  test,  the  end  of  the  clear  plastic  tubing  in 
the  pumping  piezometer  was  lowered  to  a  position  below  that 
of  the  hydraulic  head  in  the  piezometer;  this  allowed  water 
to  freely  drain  from  the  pumping  piezometer.  The  water 
draining  from  this  piezometer  was  directed  away  from  the 
test  site  and  discharged  directly  into  the  effluent  weir  of 
the  wetland  by  using  a  garden  hose;  this  prevented  the 
water  from  mounding  over  the  test  site  and  interfering  with 
the  test  results.  The  discharge  rate,  Q,  was  determined  by 
measuring  the  amount  of  water  collected  in  a  plastic 
container  of  known  volume  over  time  from  the  end  of  the 
garden  hose.  The  tests  were  run  until  the  water  levels  in 
the  observation  piezometers  ceased  to  decline.  The  data 
was  then  analyzed  using  was  Equations  (27)  through  (40) . 
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For  the  analysis,  the  top  two  soil  layers  of  the  wetland 
comprised  the  aquitard,  while  the  bottom  soil  layer  and 
underlying  gravel  and  sand  layers  served  as  a  composite 
aquifer  layer.  Estimates  of  hydraulic  conductivity  for  the 
sand  and  gravel  layers  were  assumed  based  on  the  mean  of 
values  presented  in  Charbeneau  (2000)  for  these  materials, 
and  Equation  (5b)  was  then  used  to  compute  the  hydraulic 
conductivity  of  the  bottom  soil  layer.  The  geometric  mean 
of  conductivities  from  the  four  pump  tests  conducted 
provided  a  baseline  conductivity  estimate  for  the  bottom 
soil  layer  of  the  wetland. 

Porosity .  Porosity  tests  were  conducted  on  two 
occasions  for  a  total  of  eight  samples  tested.  Eully 
saturated  soil  samples  were  collected  at  various  locations 
around  the  wetland  from  surface  sediments.  The  samples 
were  collected  in  metal  containers  of  known  weight  and 
volume  by  depressing  the  container  into  the  sediments.  The 
saturated  samples  were  covered  with  plastic  wrap  and  placed 
in  an  airtight  container  to  prevent  evaporation  losses 
until  they  could  be  weighed  in  the  laboratory.  The  samples 
were  weighed  on  a  model  AB204-S  Mettler  Toledo  scale  to  the 
nearest  1/1000^*^  of  a  gram.  The  first  batch  of  samples  was 
dried  in  an  oven  at  105°C  until  a  constant  weight  was 
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achieved  (approximately  24  hours) .  Because  these  samples 
shrunk  considerably  during  the  drying  process,  subsequent 
samples  were  dried  at  60°C  until  a  constant  weight  was 
achieved  (approximately  48  hours)  as  recommended  by  ASTM 
D2216-63T  (1963)  for  soils  high  in  organic  content.  The 
weight  of  the  dried  samples  was  subtracted  from  that  of  the 
saturated  samples  and  the  result  was  divided  by  the  density 
of  water  (at  20°C)  to  obtain  the  volume  of  water  -  and 
consequently  the  volume  of  voids  -  in  the  original  sample. 
Porosity  was  then  computed  using  Equation  (1) . 

An  additional  calculation  was  also  performed  using 

n  =  1--^  (41) 

Pm 

where  is  dry  bulk  density  [M/L^]  and  is  particle 
density  [M/L^]  .  The  equations  for  dry  bulk  and  particle 
density  are 

Pb^'^sNl  (42) 

Pm-^sl^S  (43) 

where  Ws  is  the  mass  of  the  solid  particles  [M]  ,  Vt  is  the 
original  sample  volume  [L^]  ,  and  Vs  is  the  volume  of  the 
solid  particles  [L^] .  The  mass  of  the  solid  particles  was 
determined  by  subtracting  the  weight  of  the  collection  can 
from  the  weight  of  the  dried  sample.  The  volume  of  solid 
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soil  particles  was  determined  by  subtracting  the  computed 
volume  of  water  from  the  total  sample  container  volume. 
According  to  Fetter  (1994),  a  comparison  of  Equation  (41) 
with  Equation  (1)  yields  an  internal  check  of  consistency. 

An  evaluation  of  particle  density  can  give  an  idea  of 
whether  organic  matter  is  present  in  the  sampled  sediments. 
For  example,  while  the  particle  density  of  mineral  soils  is 
often  assumed  as  2.65  g/cm^,  organic  matter  weighs  much  less 
with  particle  densities  ranging  from  1.2  to  1.5  g/cm^ 
(Charbeneau,  2000)  .  Samples  with  particle  densities  less 
than  the  mineral  average  imply  that  organic  material  may  be 
present  in  the  media. 

Specific  Yield.  The  specific  yield,  Sy,  of  the  soil 
samples  was  approximated  by  gravity  draining  the  porosity 
test  samples  prior  to  oven  drying.  The  saturated  samples 
were  drained  in  an  airtight  container  to  prevent 
evaporation  losses,  and  after  72  hours,  the  samples  were 
weighed  and  the  specific  yield,  Sy,  was  estimated  with 


^sample  ^gravity  drained )/ Pwater 
^sample 


Although  soils  have  been  observed  to  drain 
column  experiments  (Prill,  Johnson  and  Morr 


(44) 

for  months 
is,  1965), 


in 

72 
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hours  of  drainage  was  assumed  adequate  for  the  relatively 
small  sample  volumes  used  in  these  porosity  tests. 

Effective  Porosity .  The  effective  porosity,  rig,  of  the 
soil  samples  was  estimated  by  comparing  the  values  obtained 
by  gravity  draining  the  samples  with  advice  contained  in 
the  literature.  Although  a  soil's  specific  yield  -  or 
drainable  porosity  -  is  not  the  same  as  its  effective 
porosity  (it  ignores  the  volume  of  water  retained  in  the 
pores  due  to  capillary  forces) ,  it  provides  a  closer 
approximation  of  effective  porosity  than  total  porosity  for 
fine-grained  sediments.  Effective  porosities  as  low  as  0- 
5%  are  commonplace  for  soils  with  significant  amounts  of 
fine  particles  such  as  clays  and  silts  (Kresic,  1997)  . 

Nvimerical  Model 

A  three-dimensional  finite-difference,  numerical 
computer  model  was  constructed  using  Visual  MODFLOW  by 
Waterloo  Hydrogeologic,  Inc.  The  model  was  partitioned 
into  90  columns,  46  rows,  and  6  layers.  The  individual 
cells  were  evenly  distributed  throughout  the  model  with 
column  and  row  dimensions  of  IW  on  each  side.  Inactive 
cells  were  assigned  to  portions  of  the  model  that  fell 
outside  of  the  wetland  media  to  define  the  physical 
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boundaries  of  the  wetland  as  depicted  in  Figure  3-6.  The 
impermeable  wetland  liner  was  assigned  an  elevation  of  zero 
for  the  model  to  define  the  bottom  layer.  The  soil  layers 
represented  in  the  model  are  the  same  as  those  shown  in 
Figure  3-2  with  the  addition  of  a  notional  3"  soil  layer 
overlying  the  top  sediments  of  the  wetland  to  represent  the 
free  water  surface.  All  layers  were  designated  as  confined 
in  the  model  except  for  the  top,  fictional  layer. 


Figure  3-6.  Overhead  View  of  the  Wetland  Model  in  Visual 
MODFLOW®  (Waterloo  Hydrogeologic,  Inc.).  Dark  parallel 
lines  indicate  location  of  constant  head  pipes  in  the  model. 


The  piezometers  were  installed  as  observation  wells  at 
relevant  points  in  the  model.  The  weir  was  simulated  as  a 
constant  source  (i.e.  sink)  of  head  with  a  value  equal  to 
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the  elevation  of  the  water  flowing  over  it.  The  water 
distribution  pipes  were  initially  simulated  as  constant 
sources  of  head  to  facilitate  model  calibration;  the 
elevation  of  the  aboveground  water  table  surface  was  also 
initially  simulated  as  a  constant  source  of  head  with  this 
aim  in  mind.  The  distribution  pipes  were  later  simulated 
as  a  series  of  injection  wells,  and  the  constant  head 
restriction  for  the  water-surface  was  removed  to  allow  the 
model  to  solve  for  its  elevation. 

Hydraulic  conductivity  values  for  active  cells  were 
kriged  for  each  soil  layer  of  the  wetland  using  the  default 
settings  of  the  (two-dimensional)  kriging  interpolation 
option  in  Surfer®.  The  hydraulic  conductivity  values  for 
the  sand  and  gravel  layers  were  assumed  to  be  the 
arithmetic  means  of  the  range  of  values  for  these  materials 
presented  in  the  literature  (Charbeneau,  2000)  .  Hydraulic 
conductivity  values  for  the  top  (fictional  soil)  layer  were 
obtained  by  trial  and  error  until  the  model  presented  a 
logical  solution.  All  values  of  hydraulic  conductivity 
were  imported  into  the  model  using  a  Fortran  code  written 
by  Huang  (2001) .  Values  for  porosity  and  effective 
porosity  were  entered  directly  into  the  model  using  Visual 
MODFLOW' s  window  prompts.  Estimates  for  evapotranspiration 
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(Renwick,  2002)  were  also  entered  using  the  program's 
default  interface. 

The  model  simulations  were  run  using  a  convergence 
criterion  of  0.01  feet.  The  deviation  between  calculated 
and  observed  hydraulic  heads  provided  the  basis  for  model 
calibration.  Calibration  involved  adjusting  the  baseline 
hydraulic  conductivity  estimates  at  relevant  piezometer 
locations  until  a  good  fit  between  the  calculated  and 
observed  hydraulic  heads  was  achieved.  (Visual  MODFLOW® 
provides  a  comparative  graphing  feature  that  also  computes 
the  root  mean  square  error  between  the  calculated  and 
observed  hydraulic  heads  to  aid  in  model  calibration.)  The 
two  bottom  layers  of  soil  were  calibrated  in  this  manner  to 
simulate  the  hydraulic  head  distribution  observed  on 
November  1,  2001.  Calibration  of  these  layers  took  time 
but  was  relatively  simple  as  the  two  demonstrated  little 
relation  to  each  other  (i.e.  a  change  of  values  in  one 
layer  had  little  effect  on  the  values  in  the  other  layer) . 
Calibration  of  the  top  layer  proved  much  more  difficult  -  a 
change  in  the  top  layer  conductivities  produced  significant 
changes  in  the  those  of  the  middle  layer  -  and  the  measured 
values  from  hydraulic  conductivity  tests  were  kriged  and 
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used  without  further  modification  for  this  layer  of  the 
model . 

The  modeling  effort  produced  essentially  two  models. 
The  first  model  used  the  constant  head  package  to  simulate 
the  distribution  pipes,  effluent  weir,  and  the  water  table 
surface.  (The  values  of  head  in  the  pipes  were  determined 
by  trial  and  error  until  they  produced  a  good  fit  with  the 
baseline  hydraulic  conductivity  estimates  and  the 
piezometric  heads  observed  on  November  1,  2001.  The  values 
of  head  for  the  water  surface  were  determined  by  measuring 
down  from  the  known  elevation  of  the  above  ground  riser 
coupling  to  the  water  surface,  and  the  value  at  the  weir 
was  determined  in  a  similar  manner.)  These  conditions 
produced  a  rapid  solution  and  in  this  respect  facilitated 
model  calibration,  but  they  were  somewhat  restrictive:  they 
essentially  forced  the  water  in  the  model  to  flow  from  a 
position  of  higher  head  in  the  pipes  to  the  lower  head  at 
the  surface  of  the  wetland  and  in  the  weir. 

To  test  the  calibrated  model  under  less  restrictive 
conditions,  the  constant  head  for  the  water  table  surface 
was  removed.  This  allowed  the  model  to  determine  the 
proper  location  of  the  water  table  surface  based  on  the 
water  inputs  and  other  boundary  conditions.  The  pipes  were 
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also  simulated  as  a  series  of  injection  wells  rather  than 
constant  sources  of  head.  Sporadic  readings  from  the 
inflow  meter  display  were  averaged  to  obtain  a  mean  value 
of  flow  into  the  wetland;  the  output  per  well  was 
determined  by  evenly  dividing  this  value  by  the  number  of 
wells  used  to  simulate  the  distribution  pipes.  Only  one 
solver  -  the  Preconditioned  Conjugate-Gradient  Package 
(PCG2)  using  the  Modified  Incomplete  Cholesky  pre¬ 
conditioning  method  -  produced  logical  results  for  these 
conditions.  The  model  required  approximately  five  minutes 
to  solve  using  this  solver  package  for  the  less  restrictive 
boundary  conditions. 

Validation  was  attempted  by  reducing  the  flow  rate 
into  the  wetland  to  test  whether  the  calibrated  hydraulic 
conductivities  produced  results  similar  to  that  for  the 
baseline  model  conditions.  The  piezometric  data  obtained 
on  February  19,  2002  were  used  for  this  purpose. 

Visual  MODFLOW' s  particle  tracking  feature  was  used  to 
estimate  the  hydraulic  residence  time  of  particles  placed 
along  the  simulated  distribution  pipes. 
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IV.  Results 


Sampling  Grid 

The  sampling  grid  provided  the  resolution  necessary  to 
reveal  spatial  trends  in  the  distribution  of  hydraulic  head. 
The  existence  of  significant  soil  heterogeneities  precluded 
a  rigorous  application  of  flow  net  construction;  however, 
on-site  observations  provided  valuable  insight  to  confirm 
what  the  hydraulic  head  contours  imply.  The  piezometers  in 
the  grid  also  facilitated  measurements  of  hydraulic 
conductivity  for  numerical  modeling  and  calculations  of 
groundwater  flow  rates  in  the  wetland. 

Piezometer  Installation 

Although  the  drive  point  method  of  installation  was 
found  to  be  more  challenging  than  expected,  all  of  the 
piezometers  -  with  the  exception  of  number  29C  which  hit  a 
rock  -  were  successfully  hand  driven  and  retracted  to  the 
target  depth  below  the  surface  of  the  wetland.  The  results 
of  the  total  station  survey  (Appendix  A)  reveal  a  standard 
deviation  about  the  mean  piezometer  intake  elevations  of 
0.83",  0.71",  and  0.74"  for  the  top,  middle  and  bottom  soil 
layers,  respectively.  The  general  method  of  emplacement 
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(i.e.  assuming  a  level  wetland  surface,  driving  the 
piezometer  point  down  to  a  set  distance,  and  then  pulling 
back  6"  to  extract  the  shield  from  the  screened  portion  of 
the  piezometer)  worked  well  enough  to  achieve  an  acceptable 
variance  of  sampling  point  elevation  in  the  wetland.  (The 
assumption  of  a  level  surface  was  based  on  the  construction 
history  of  the  wetland  as  well  as  a  hasty  elevation  survey 
conducted  prior  to  installing  the  piezometer  grid.)  The 
small  variance  in  piezometer  intake  elevations  facilitated 
meaningful  interpretation  of  the  potentiometric  contours 
for  each  soil  layer  in  the  wetland. 

A  few  observations  concerning  the  installation  of  the 
piezometers  may  assist  with  similar  groundwater  studies  in 
wetlands.  The  bentonite  seal  greatly  improved  the 
stability  of  the  piezometers  in  the  soft  wetland  sediments, 
but  its  potential  impact  on  the  monitoring  objectives  of 
researchers  studying  contaminant  degredation  is  unknown. 

The  presence  of  a  connector  coupling  located  below  the 
surface  of  the  wetland  on  some  of  the  bottom  layer 
piezometers  also  seems  to  have  markedly  improve  stability: 
the  piezometers  with  these  couplings  were  significantly 
less  prone  to  vertical  movement  and  settling  than  those 
constructed  with  just  a  single  riser  tube.  Last,  as 
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already  acknowledged,  an  annular  cap  of  indigenous  wetland 
soil  may  have  provided  a  more  enduring  seal  than  the 
Quickcrete®  caps  used  given  the  observation  that  many  of 
the  latter  cracked  shortly  after  hydration  and  did  not 
provide  an  optimal  covering. 

Developing  the  piezometers  seems  to  have  altered  the 
hydrogeologic  properties  surrounding  the  intake  portion  of 
the  piezometers.  Prior  to  well  development,  slug  tests 
performed  in  the  top  and  middle  soil  layers  yielded  mean 
coefficients  of  determination  (R^)  of  0.995  and  0.973, 
respectively.  Following  well  development,  the  mean 
values  for  the  top  and  middle  layers  were  0.948  and  0.911, 
respectively.  In  other  words,  the  data  from  the  slug  tests 
performed  prior  to  piezometer  development  closely  resembled 
the  ideal  straight-line  appearance  of  the  Bouwer  and  Rice 
semi-log  plots  while  deviations  from  a  straight-line  fit 
after  piezometer  development  were  commonplace.  Such 
deviations  were  characterized  by  a  rapid  water-level 
recovery,  generally  during  the  first  minute  of  the  test;  in 
some  cases,  a  straight-line  plot  of  data  points  on  the 
semi-logarithmic  plot  did  not  appear  until  towards  the  very 
end  of  the  test.  Figure  4-1  illustrates  these  points  for  a 
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comparison  slug  test  conducted  before  and  after  piezometer 


development . 
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Figure  4-1.  General  Effect  of  Developing  the  Piezometers.  Later 
date  is  test  conducted  on  developed  piezometer;  note  the  deviation 
from  the  straight  line  appearance  and  more  rapid  rate  of  water  level 
recovery  relative  to  the  earlier  date. 


The  rapid  water  level  recovery  was  not  due  to  the 
presence  of  a  higher  permeability  sand  or  gravel  pack 
surrounding  the  piezometer  screen  as  one  does  not  exist; 
the  possibility  that  sediments  may  have  been  clogging  the 
piezometer  screens  prior  to  development  is  also  small  owing 
to  the  use  of  the  piezometer  shields  during  installation. 

A  more  plausible  explanation  is  that  developing  the 
piezometers  altered  the  soil  properties  in  close  proximity 
to  the  piezometer  screen;  alternatively,  development  may 
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also  have  disturbed  the  bentonite  seal  around  the  riser 
pipe . 

The  appearance  of  data  points  on  the  semi-log  plot  was 
not  the  only  change  observed  in  slug  tests  as  a  result  of 
well  development.  More  significantly,  comparisons  of  slug 
tests  conducted  before  and  after  piezometer  development  on 
19  piezometers  demonstrated  order  of  magnitude  increases  in 
computed  conductivities;  seven  of  these  differed  by  more 
than  two  orders  of  magnitude  while  only  three  of  the  19 
increased  by  less  than  one  order  of  magnitude. 

The  possibility  that  developing  the  piezometers  may 
have  altered  the  soil  properties  in  the  area  immediately 
surrounding  the  intake  screens  is  based  on  the  elastic 
properties  of  soils.  Demir  and  Narasimhan  (1994) 
demonstrate  that  data  points  for  the  Hvorslev  test  deviate 
from  the  theoretical  straight  line  for  certain  piezometer 
shape  factors  as  specific  storage,  Ss,  increases,  especially 
where  radial  flow  through  the  piezometer  screen  occurs. 
Although  the  author's  findings  pertain  to  the  Hvorslev  test, 
the  analytical  form  of  the  Bouwer  and  Rice  and  Hvorslev 
equations  are  similar  (Charbeneau,  2001)  and  the  device 
intake  geometries  in  question  both  involve  radial  flow. 

Demir  and  Narasimhan  (1994)  make  it  clear  that  deviations 
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are  commonplace  in  unconsolidated  or  "soft  sediments"  where 
Ss  can  exceed  values  of  0.001  m"^ .  As  presented  in  chapter 
2,  specific  storage,  Ss,  is 

Ss^Pw  9(c(  +  n /3)  (20) 

where 

Pw  =  density  of  water  [M/L^] 
g  =  gravity  [L/T^] 

a  =  compressibility  of  the  aquifer  skeleton  [1/(M/LT^)] 
n  =  porosity  [L^/  L^] 

P  =  compressibility  of  water  [1/M/LT^]] 

While  the  variables  ,  g  and  (3  are  unlikely  candidates  for 
change,  an  increase  in  either  n  or  «  may  have  occurred  as  a 
result  of  surging  too  much  water  into  the  piezometer  during 
development  (assuming  the  media  immediately  surrounding  the 
piezometer  intake  had  an  opportunity  to  expand  and  deform) ; 
the  resultant  increase  in  specific  storage  relative  to  the 
undeveloped  condition  could  account  for  the  characteristic 
deviations  from  straight-line  appearance  observed  on  the 
semi-log  plot  that  were  common  after  developing  the 
piezomters.  Likewise,  although  hydraulic  conductivity  does 
not  depend  on  either  n  or  « ,  both  parameters  are 
represented  in  the  intrinsic  properties  and  permeability,  k, 
of  the  soil,  so  an  increase  in  either  n  or  a  could  account 
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for  the  observed  order  of  magnitude  increases  in  hydraulic 
conductivity  after  piezometer  development. 

It  is  unlikely  that  the  bentonite  seal  was  completely 
compromised  as  the  result  of  developing  the  piezometers:  an 
ineffective  annular  seal  would  manifest  itself  in  static 
water  levels  equal  to  the  water  table  elevation  in  the 
presence  of  vertical  groundwater  flow  gradients  (Sigel  and 
Glaser,  1987) .  Developing  may  have  compromised  a  portion 
of  the  seal  close  to  the  piezometer  screen  however. 
Expansive  clays  such  as  bentonite  can  have  porosities  that 
range  as  high  as  80%  (Kresic,  1997),  and  as  indicated  by 
Equation  (20) ,  a  high  porosity  can  result  in  a  high  value 
of  specific  storage  (Ss)  .  The  data  point  deviation  from  a 
straight-line  appearance  could  indicate  that  some  of  the 
water  introduced  or  extracted  during  the  slug  tests  was 
moving  into  or  out  of  the  highly  plastic  bentonite  material 
but  was  not  fully  penetrating  the  annular  seal. 

Piezometric  Surface 

The  results  of  the  piezometric  measurements  indicate 
areas  of  preferential  flow  in  the  wetland.  Potentiometric 
surfaces  for  November  1,  2001  and  Eebruary  19,  2002  are 
given  in  appendix  B  and  C,  respectively.  These  contour 
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plots  indicate  a  fairly  uniform  distribution  of  hydraulic 


head  in  the  bottom  soil  layer  relative  to  the  two  upper 
layers.  The  potentiometric  surfaces  for  the  two  upper 
layers  depict  regions  of  higher  hydraulic  head  and 
preferential  flows  towards  the  north  side  and  weir  end  of 
the  wetland.  (Figure  3-1  indicates  which  direction  is 
north  relative  to  the  wetland.)  The  behavior  observed  in 
Appendix  B  was  typical  of  that  observed  on  seven  occasions 
over  a  three-week  period  around  the  date  represented.  The 
magnitude  of  variation  observed  in  Appendix  C  appears  less 
than  that  observed  in  Appendix  B  due  to  a  reduced  flow  rate 
into  the  wetland.  (Appendix  D  contains  the  water  level 
data  collected  on  all  days  measured.) 

The  variability  depicted  in  Appendices  B  and  C  is 
quantified  in  table  4-1.  Although  the  magnitude  of 
variation  for  all  layers  might  not  be  significant  if  it 
occurred  gradually  over  the  full  length  of  the  wetland. 


Table  4-1.  Comparison  of  Hydraulic  Head  Distribution 
for  November  1,  2001  and  February  19,  2002. 


11/1/2001 

2/19/2002 

Top 

Middie 

Bottom 

Top 

Middie 

Bottom 

Mean  (ft) 

6.2 

6.61 

8.04 

5.983 

6.249 

7.616 

Standard  Deviation  (ft) 

0.335 

0.469 

0.093 

0.18 

0.351 

0.068 

Maximum  Difference  in  Head  (ft) 

1.657 

2.015 

0.581 

0.812 

1.47 

0.319 

Maximum  Hydrauiic  Gradient 

0.108 

0.123 

0.055 

0.065 

0.138 

0.019 
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differences  in  head  as  large  as  12  to  18  inches  are  common 
in  piezometers  located  close  to  each  other  in  the  top  two 
soil  layers.  Additionally,  the  large  differences  are  not 
isolated  to  single  piezometers  (which  could  then  be 
dismissed  as  mere  anomalies  or  outliers) ;  rather  they  occur 
between  adjacent  regions  or  groups  of  piezometers. 

The  case  for  preferential  flow  is  further  supported  by 
observations  on  the  surface  of  the  wetland  near  the 
locations  of  higher  hydraulic  head.  The  sediments  around 
the  perimeter  of  the  wetland,  particularly  along  the  north, 
are  highly  plastic  and  generally  do  not  support  the  weight 
of  foot  traffic  well.  The  sediments  down  the  center  of  the 
wetland,  from  piezometers  3  and  4  to  33  and  34,  are  also 
very  soft.  (This  contrasts  with  other  areas  of  the  wetland 
where  the  sediments  generally  support  foot  traffic  well.) 
Additionally,  moisture  above  the  water  table  and  close  to 
the  liner  is  evident  between  piezometers  54  and  60.  Figure 
4-2  displays  these  areas. 

The  most  convincing  observation  of  preferential  flow 
is  not  entirely  supported  by  the  piezometric  data  however. 
The  presence  of  boiling  or  heave  is  evident  approximately 
two  feet  to  the  north  of  nest  29.  Generally  attributed  to 
cohesionless  soils  like  sand  or  silt,  boiling  or  heave  can 


4-9 


Figure  4-2.  Regions  of  Particularly  Poor  Trafficability  in  the  Wetland. 
Shaded  regions  represent  areas  of  suspected  preferential  flows. 

Arrow  pointing  away  from  nest  29  indicates  observable  direction  of 
flow  off  of  the  mound. 

occur  in  the  presence  of  vertical  gradients  that  approach  a 
value  of  unity  (Charbeneau,  2000) .  The  generic  name 
applied  to  this  condition  is  quicksand .  The  phenomenon 
occurs  when  the  seepage  forces  acting  up  on  a  soil  mass 
balance  the  weight  of  the  media  and  water  acting  down  upon 
it.  Although  the  hydraulic  gradient  at  nest  29  does  not 
approach  one,  the  appearance  of  the  sediments  erupting  from 
upward  pressures  at  this  point  -  like  a  bulge  in  a  tire  - 
is  unmistakable.  The  wetland  sediments  that  are  void  of 
roots  around  this  area  will  not  support  any  foot  traffic 
without  the  aid  of  wide  boards  to  distribute  the  load. 
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Additionally,  water  can  actually  be  seen  flowing  out  of  the 


ground  and  down  the  mound  of  eruption  at  this  location. 

The  fact  that  none  of  the  piezometers  register  a  gradient 
of  one  in  this  particularly  wet  area  of  wetland  merely 
indicates  that  the  sampling  grid  has  failed  to  pinpoint  the 
exact  location  where  the  breach  has  occurred. 

Table  4-2  provides  a  measure  of  the  variability  within 
and  between  layers  for  all  dates  measured.  Although  the 
sample  sizes  (displayed  in  parentheses  next  to  the  column 
headings)  are  too  small  for  convincing  statistical  analysis, 
they  results  may  imply  that  a  reduced  flow  rate  reduces  the 
variation  of  hydraulic  head  within  layers. 


Table  4-2.  Mean  Standard  Deviation  in  Hydraulic 
Head  Distribution  for  All  Dates  Measured. 


Fall  2001  (7) 

Winter  2002  (1) 

Mean  (8) 

Top  Soil  Layer  (ft) 

0.376 

0.18 

0.352 

Middle  Soil  Layer  (ft) 

0.572 

0.351 

0.544 

Bottom  Layer  (ft) 

0.107 

0.068 

0.102 

Note:  sample  size  displayed  in  parentheses  next  to  the  column  heading. 


Parameter  Analysis 

Hydraulic  Conductivity .  Appendix  E  displays  the 
results  of  the  slug  tests  performed.  The  results  range 
over  values  indicated  by  various  groundwater  texts  for 
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silts,  sandy  silts,  clayey  sands,  and  tills  to  silty  sands 


and  fine  sands  (Fetter,  1994;  Charbeneau,  2000)  .  Overall, 
the  results  indicate  heterogeneities  and  variations  within 
and  between  layers  as  seen  in  Table  4-3. 


Table  4-3.  Hydraulic  Conductivity  Estimates  from  In-situ  Tests  (ft/sec) 


Top  Soil  Layer 

Middle  Soil  Layer 

Bottom  Soil  Layer 

Arithmetic 

Mean 

Standard 

Deviation 

Arithmetic 

Mean 

Standard 

Deviation 

Arithmetic 

Mean 

Standard 

Deviation 

Before  Developing 

6.4x10® 

4.9x10® 

1.5x10® 

3.1x10® 

6.8x1 0"^ 

3.4x1 0"^ 

After  Developing 

8.1x10® 

1.1x10® 

5.6x10® 

8.0x10® 

6.8x1 0"^ 

3.4x1 0"^ 

In  Table  4-3,  the  magnitude  of  standard  deviations  relative 
to  the  means  indicate  variation  within  layers  while  the 
means  demonstrate  variation  between  layers. 

The  impact  of  developing  the  piezometers  can  be  seen 
in  Table  4-4  by  comparing  the  geometric  means  of  the  layers 
before  and  after  development.  (Many  hydrologists  believe 
the  geometric  mean  provides  a  better  comparison  assessment 
of  mean  hydraulic  conductivities  because  an  arithmetic  mean 
tends  to  give  more  weight  to  larger  values  (Fetter,  1994) .) 


Table  4-4.  Geometric  Means  of  Hydraulic  Conductivity  (ft/sec) 


Top  Soil  Layer 

Middle  Soil  Layer 

Bottom  Soil  Layer 

Before  Developing 

4.2x10® 

2.5x10'^ 

5.6x1 0"^ 

After  Developing 

3.4x10® 

2.2x10® 

5.6x1 0"^ 
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Although  hydraulic  conductivity  appears  to  increase  with 
depth,  the  trend  is  not  conclusive.  Prior  to  development, 
hydraulic  conductivity  clearly  improved  with  depth,  while 
after  development,  the  (geometric)  mean  conductivities  of 
the  top  and  middle  layers  were  approximately  equal. 

The  concern  over  the  choice  of  the  appropriate  slug 
test  for  this  study  was  alleviated  somewhat  by  comparing 
the  results  obtained  from  the  Bouwer  and  Rice  (1976)  and 
Hvorslev  (1951)  methods;  results  obtained  with  the  Hvorslev 
method  using  a  shape  factor  that  may  approximate  that  of 
the  piezometers  used  in  this  study  (Brand  and  Permichitt, 
1980)  over-performed  the  results  obtained  from  the  Bouwer 
and  Rice  test  by  only  30%  on  average.  (See  Appendix  F  for 
an  example  of  the  MATHCAD®  template  used  to  compute 
hydraulic  conductivity  and  compare  the  results  of  both  test 
methods.)  The  close  results  may  not  be  surprising  given 
the  similar  form  of  the  key  equations  as  well  as  the 
similar  piezometer  intake  geometries,  although  it  did 
provide  some  degree  of  reassurance.  Also,  a  30%  disparity 
is  well  within  the  order  of  magnitude  precision  that  is 
possible  for  hydraulic  conductivity  tests. 

The  hydraulic  conductivity  estimates  shown  in  Tables 
4-3  and  4-4  for  the  bottom  soil  layer  reflect  the  results 
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of  the  pump  tests  conducted  on  four  piezometers  in  the 
layer.  The  fit  for  the  data  from  two  of  the  tests  implied 
some  leakage  from  the  overlying  aquitard  while  the  other 
two  did  not.  The  effects  of  partial  penetration  were 
negligible  for  all  four  tests.  Appendix  G  contains 
calculations  and  plots  for  one  of  the  four  tests. 

Porosity .  The  results  of  the  porosity  tests  range 
from  0.46  to  0.63  (sample  mean:  0.53) .  The  low  end  was 
typical  of  samples  taken  from  the  compacted  portion  of  the 
wetland,  while  the  higher  values  were  common  in  samples 
taken  from  wet  areas  in  the  wetland.  The  values  in  this 
range  are  typical  of  silts  and  clays  (Charbeneau,  2000) . 

Specific  Yield  and  Effective  Porosity .  The  results  of 
the  specific  yield  tests  range  from  2.2x10”^  to  0.069 
(sample  mean:  0.022.)  Because  these  results  do  not  account 
for  the  water  retained  by  capillary  forces  in  the 
interconnected  void  spaces,  a  value  of  0.05  -  which  is 
higher  than  the  mean  for  specific  yield  but  representative 
of  soils  containing  significant  amounts  of  fine  particles  - 
was  assumed  for  effective  porosity. 

Particle  Density .  The  results  of  the  particle  density 
analysis  range  from  2.2  to  2.4  g/cm^.  These  values  indicate 
a  presence  of  organic  material  in  the  wetland  sediments. 
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Nvimerical  Model 


The  output  from  the  model  implies  a  tendency  for  water 
to  preferentially  flow  towards  the  sides  of  the  wetland  in 
the  more  permeable  layers.  Flow  in  the  bottom  soil  and 
underlying  gravel  and  sand  layers  is  predominately 
horizontal,  while  water  in  the  top  two  soil  layers  flows 
predominately  straight  up.  This  is  illustrated  in  Figure 
4-3.  The  wider  spacing  of  contour  intervals  along  the 
sides  of  the  wetland  indicates  that  the  sediments  have  a 
higher  transmissivity  and  facilitate  greater  flows  relative 
to  the  more  central  portions  of  the  wetland.  In  general, 
the  output  from  the  model  indicates  that  a  greater  amount 
of  water  moves  towards  the  north  side  of  the  wetland. 

Visual  MODFLOW' s  particle  tracking  feature  computed  an 
average  hydraulic  residence  time  of  three  days  for  water 
particles  released  from  the  distribution  pipes.  The 
minimum  residence  time  was  approximately  16.5  hours.  This 
generally  occurred  in  the  vicinity  of  piezometers  10,  24, 
29,  30,  and  36  based  on  observations  of  the  model  output. 
The  maximum  residence  time  for  particles  released  was  15 
days  while  approximately  64%  of  the  particles  released  had 
a  residence  time  less  than  the  average.  The  large 
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(Near  Piezometer  Nests  25-30) 


SCIS)  C(19)  NC19> 


(Near  Piezometer  Nests  55-60) 


Figure  4-3.  Cross  Section  Views  of  the  Wetland  (Displays  from 
Visual  MODFLOW®) 

variation  in  computed  residence  time  implies  preferential 
flow  in  the  wetland. 

The  program' s  water  budget  calculation  feature 
indicated  that  water  losses  from  evapotranspiration  were 
minor,  amounting  to  less  than  five  percent  of  the  wetland' s 
overall  water  budget  for  the  baseline  conditions.  This  was 
based  on  a  potential  evapotranspiration  rate  of  10  mm/day 
for  November  1,  2001  (Renwick,  2002).  Evapotraspiration 
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was  assumed  to  be  zero  for  the  validation  measurements 
taken  during  February  2002. 

Approximately  50  calibration  runs  were  required  to 
achieve  a  close  fit  between  computed  and  observed  hydraulic 
head  values  for  the  bottom  two  layers  of  soil.  (As  noted 
in  Chapter  3,  hydraulic  conductivity  was  not  calibrated  for 
the  top  layer  of  soil.)  The  fit  achieved  for  these  layers 
is  displayed  in  Appendix  H.  The  general  method  of 
adjusting  hydraulic  conductivity  for  cells  at  piezometer 
locations  in  the  model,  and  then  kriging  these  values  to 
interpolate  conductivities  for  the  other  cells  in  the  layer, 
was  not  100%  successful  in  bringing  all  points  close  to  the 
ideal  fit  line.  Most  notably,  piezometers  4C,  25C,  IB,  40B, 
41B  and  61B  were  particularly  problematic  and  continually 
overestimated  by  the  model.  This  helps  to  illustrate  the 
point  that  although  groundwater  models  can  be  helpful  in 
characterizing  the  general  behavior  of  a  subsurface  system, 
they  are  generally  not  capable  of  capturing  all  geological 
features  hidden  within  a  natural  system.  Numerical  models 
generate  approximate  solutions,  and  the  weight  placed  on 
such  solutions  must  be  viewed  with  this  truth  in  mind. 

Figure  4-4  displays  the  overall  fit  achieved  during 
model  calibration  for  the  conditions  and  heads  observed  on 
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Calculated  vs.  Observed  Heads  ;  Steady  state 


•• 


5.8  6.8  7.8 

Obs.  Heads  (ft) 

Num.  Fbints  :  1 83 


Mean  Error  :  -0.09226004  (ft) 
Mean  Absolute  :  0.1446608  (ft) 
Standard  Error  of  the  Btimate  :  0.018  85738  (ft) 
Root  mean  squared  :  0.27  06128  (ft) 
Normalized  RMS:  12.00855  (%) 


Figure  4-4.  Calibrated  Model.  Grouping  of  points  in  the  upper  right  corner 
represent  piezometers  in  the  bottom  soil  layer.  Points  forming  a  horizontal 
line  along  the  bottom  represent  piezometers  in  the  top  soil  layer  (not 
calibrated).  Points  between  these  two  groups  represent  piezometers  in  the 
middle  soil  layer. 


November  1,  2001.  A  comparison  of  Figure  4-4  with  Figure 
4-5  provides  an  assessment  of  how  well  the  calibrated  model 
was  able  to  predict  the  validation  conditions  and  heads 
observed  on  February  19,  2002.  The  deviation  from  an  exact 
fit  may  indicate  calibration  error,  errors  in  parameter 
estimation,  or  transient  subsurface  conditions.  The 
majority  of  data  points  in  the  center  of  Figure  4-5  (which 
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Calculated  vs.  Observed  Heads  :  Steady  state 


fi 


•  •  •  • 


•  *  •  • 

V  ^9 

^  *  Validation 

0^  Conditions 


(02-19-02) 


5.7 


6.7 

Obs.  Heads  (ft) 


7.7 


Num.  Fbints  :  1 97 
Mean  Bror  :  0.03247213  (ft) 
Mean  Absolute  :  0.1597922  (ft) 
Standard  Error  of  the  Btimate  :  0  .01858658  (ft) 
Root  mean  squared  :  0.2622304  (ft) 
Namalized  RMS  :  12.79797  (%) 


Figure  4-5.  Model  Validation.  The  middle  soil  layer  piezometers 
appear  to  bracket  a  line  slightly  above  the  best-fit  line. 


represent  the  middle  soil  layer)  do  appear  to  bracket  a 
straight  line,  and  this  provides  for  some  degree  of 
confidence  in  the  model.  The  position  of  these  points 
above  the  best  fit  line  may  indicate  that  the  mean 
hydraulic  conductivity  for  this  layer  (or  throughout  the 
wetland)  has  shifted,  or  it  could  possibly  reflect  the 
influence  of  the  colder  temperatures. 
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V.  Conclusions  and  Recommendations  for  Further  Study 


The  Wright-Patterson  Air  Force  Base  wetland  treatment 
cells  are  designed  to  achieve  uniform  vertical  flows  by 
distributing  water  into  a  highly  permeable  layer  near  the 
bottom  of  the  wetland  and  forcing  it  up  through  the  largely 
anaerobic  subsurface  sediments.  Conceptually,  uniform 
flows  will  achieve  the  most  efficient  degree  of  contaminant 
removal  possible  by  evenly  dispersing  the  groundwater 
contaminants  throughout  the  full  volume  of  the  subsurface 
media.  The  existence  of  uniform,  predictable  flows  would 
also  aid  in  the  development  of  design  criteria  for  such  a 
treatment  technology. 

This  study  shows  that  the  flow  of  groundwater  is  not 
entirely  uniform  in  wetland  cell  number  one.  Admittedly, 
what  constitutes  uniform  flows  in  such  a  system  is  a  matter 
of  perspective.  Aside  from  being  a  subjective  label,  even 
the  seemingly  variable  distribution  of  hydraulic  heads  in 
the  middle  soil  layer  could  be  scaled  in  such  a  manner  to 
appear  uniform.  Still,  the  variability  in  hydraulic  heads 
for  the  top  and  middle  soil  layers  was  three  to  five  times 
greater  than  that  for  the  bottom  soil  layer.  The  fact  that 
differences  in  head  of  feet  are  common  between  groups 
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of  piezometers  located  in  close  proximity  to  one  another  - 
rather  than  being  gradually  dispersed  over  the  full  range 
of  the  wetland  -  supports  the  case  for  preferential  flow. 

The  high  values  of  hydraulic  conductivity  measured  in  these 
areas  are  not  surprising  from  an  analytical  perspective  and 
add  weight  to  this  possibility.  The  output  from  the 
numerical  model  further  supports  the  notion  that  regions  of 
higher  flow  are  inter-dispersed  among  regions  of  nominal 
flows.  The  residence  time  of  water  particles  observed  in 
the  model  is  significantly  less  in  regions  of  the  wetland 
where  the  observed  heads  and  hydraulic  conductivities  are 
high.  Additionally,  the  flow  vectors  in  the  model  indicate 
that  groundwater  in  the  more  permeable  bottom  layers 
prefers  to  flow  horizontally  towards  the  sides  of  the 
wetland,  and  higher  hydraulic  heads  are  observed  along  the 
north  side  and  weir  end  of  the  potentiometric  contour  maps. 
Finally,  the  observations  noted  on  the  surface  of  the 
wetland  also  support  the  case  for  preferential  flows. 

The  above  observations  do  not  imply  that  only  a  small 
portion  of  the  wetland  facilitates  flow.  The  existence  of 
water  levels  above  the  surface  of  the  water  table  in 
virtually  all  of  the  functioning  piezometers  indicates  that 
upward  vertical  flows  occur  throughout  the  subsurface  media. 
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The  numerical  model  supports  this  evidence  by  displaying 
vertical  flow  vectors  and  water  particle  movement  through 
the  top  and  middle  soil  layers.  A  large  proportion  of  the 
wetland  media  may,  therefore,  be  performing  as  intended 
albeit  at  a  somewhat  reduced  level  of  efficiency. 

Study  Strengths 

The  study  used  well-established  methods  to  investigate 
groundwater  flow  in  wetland  cell  number  one.  The  evidence 
from  these  different  methods  is  consistent  and  supports  the 
finding  that  some  degree  of  preferential  flow  is  occurring. 
At  the  same  time,  evidence  from  the  piezometers  and  model 
simulations  confirm  that  vertical  flows  also  occur  through 
at  least  36"  of  saturated  wetland  soil. 

Study  Weaknesses 

The  effect  of  the  suspected  preferential  flows  on  the 
contaminant  degradation  processes  at  work  in  the  wetland 
was  not  determined  in  this  study.  The  eight  measurements 
of  hydraulic  head  taken  during  this  study  do  not  permit 
meaningful  statistical  analysis,  and  the  short  time  span 
over  which  the  measurements  were  taken  does  little  to 
indicate  whether  the  observed  flow  patterns  reflect  a 
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steady  state  condition  or  are  merely  snapshots  of  a  flow 
system  in  transition.  The  question  of  how  varying  loading 
rates  affect  patterns  in  groundwater  flow  was  not  studied 
as  fully  as  intended.  The  method  used  to  develop  the 
piezometers  was  not  well  established  in  the  literature,  and 
the  reason  behind  the  order-of-magnitude  improvements  that 
were  typical  after  development,  as  well  as  the  deviation 
from  straight-line  appearance  in  the  data  plots,  was  never 
fully  resolved.  The  reliance  on  specific  yield  estimates 
to  assume  an  effective  porosity  (which  provides  the  basis 
of  hydraulic  residence  time  in  the  numerical  model)  lacks 
precision,  although  the  effective  porosity  parameter  is 
inherently  difficult  to  measure.  Possible  strategies  to 
mitigate  or  discourage  preferential  flows  are  not  discussed 
in  this  study,  and  the  flow  behavior  in  the  wetland  cell 
number  two  was  not  studied.  Last,  the  measured  hydraulic 
conductivity  values  required  a  significant  amount  of 
adjustments  to  calibrate  the  numerical  model,  and  the  fit 
obtained  with  the  validation  conditions,  while  close,  was 
less  than  exact. 
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Recommendations  for  Further  Study 

(1)  Study  the  flow  in  wetland  cell  number  two  to  determine 
how  it  compares  to  that  observed  in  this  study. 

(2)  Conduct  a  tracer  study  in  wetland  cell  number  one  to 
substantiate  or  refute  the  results  obtained  in  this  study. 

(3)  Continue  to  monitor  the  distribution  of  hydraulic 
heads  in  wetland  cell  number  one  to  determine  whether 
changes  are  occurring  over  time. 

(4)  Vary  the  flow  into  the  wetland  to  determine  whether 
preferential  flows  occur  under  all  practical  loading 
conditions.  Attempt  to  relate  soil  parameters  to  the 
findings  to  develop  design  guidelines  appropriate  for  like 
systems . 

(5)  Improve  on  the  numerical  model  by  increasing  the  grid 
resolution  and  by  conducting  additional  parameter  tests  to 
verify  the  model's  results.  Use  the  model  as  a  predictive 
tool  to  investigate  the  flow  behavior  of  different  design 
approaches . 
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Appendix  A:  Piezometer  Grid  and  Construction  Data 


Local  Coordinates  (ft) 

Elevations 

(ft) 

In-gr 

Piez  # 

X 

y 

coupling 

screen 

tube 

Lengt 

OlA 

15.902 

10.028 

6.094 

4 . 646 

6.995 

1  ' 

OIB 

14 . 854 

10.984 

6.610 

3.694 

8 . 616 

30" 

OIC 

14 . 928 

10.039 

6.638 

2.253 

9.550 

2  'x2 

02A 

16.072 

18.113 

6.028 

4 . 627 

6.945 

1  ' 

02B 

14 . 995 

19.190 

6.020 

3.619 

7 . 822 

2  ' 

02C 

14 . 973 

18.195 

6.560 

2 . 170 

9.425 

4  ' 

03A 

15.863 

26.835 

5.948 

4.532 

6.782 

1  ' 

03B 

14 .895 

27.794 

5.982 

3.566 

7.795 

2  ' 

03C 

14.890 

26.873 

6.481 

2 .091 

9.403 

4  ' 

04A 

16.053 

35.007 

6.046 

4 . 650 

6.895 

1  ' 

04B 

14 . 987 

36.085 

6.029 

3.644 

7.389 

2  ' 

04C 

14 . 948 

35.041 

6.615 

2 . 121 

9.344 

3  'xl 

05A 

15.722 

43.392 

5.995 

4.594 

6.860 

1 ' 

05B 

14 .808 

44.247 

6.065 

3.669 

7 . 977 

2  ' 

05C 

14 . 633 

43.319 

6.571 

2 . 077 

9.358 

3  'xl 

0  6A 

15.786 

51.436 

6.111 

4 . 695 

6.992 

1 ' 

06B 

14 .897 

52.229 

6.091 

3.685 

7 . 925 

2  ' 

06C 

14 . 870 

51.358 

6.599 

2 .198 

9.584 

2  'x2 

07A 

26.509 

10.162 

6.030 

4 . 634 

6.890 

1  ' 

07B 

25.465 

11.280 

6.070 

3.674 

7 .893 

2  ' 

07C 

25.670 

10.153 

6.558 

2 . 152 

9.454 

4  ' 

08A 

26.559 

18.357 

6.098 

4.702 

6.952 

1  ' 

08B 

25.690 

19.058 

6.171 

3.739 

7 . 953 

2  ' 

08C 

25.693 

18.228 

6.620 

2.224 

9.454 

4  ' 

0  9A 

26.412 

26.823 

6.113 

4.702 

6.926 

1  ' 

09B 

25.430 

27 .898 

6.130 

3.719 

8 . 026 

2  ' 

09C 

25.552 

26.901 

6.614 

2.218 

9.541 

4  ' 

lOA 

26.291 

34 .882 

6.205 

4.789 

7 . 039 

1  ' 

lOB 

25.411 

35.918 

6.167 

3.761 

8.204 

2  ' 

IOC 

25.461 

35.019 

6.770 

2.286 

9.583 

3  'xl 

llA 

26.225 

43.166 

6.118 

4.702 

6.957 

1 ' 

IIB 

25.295 

44.224 

6.117 

3.716 

8 . 076 

2  ' 

lie 

25.375 

43.207 

6.678 

2.215 

9.480 

3  'xl 

12A 

26.168 

51 . 604 

6.064 

4 . 663 

6.918 

1 ' 

12B 

25.230 

52.563 

6.090 

3.689 

7 .991 

2  ' 

12C 

25.211 

51 . 681 

6.564 

2 . 174 

9.465 

4  ' 

Riser 

Type 

S 

S 


Local  Coordinates  (ft) 

Elevations 

(ft) 

In-gr 

Piez  # 

X 

y 

coupling 

screen 

tube 

Lengt 

13A 

37.019 

10.086 

6.113 

4.712 

7 . 020 

1  ' 

13B 

35.924 

11 . 092 

6.572 

3.682 

8.593 

30" 

13C 

35.904 

10.083 

6.671 

2 . 171 

9.661 

2  'x2 

14A 

36.869 

18 . 623 

6.125 

4.703 

6.969 

1  ' 

14B 

35.949 

19.545 

6.531 

3.646 

8.552 

30" 

14C 

35.853 

18.766 

6.672 

2 . 178 

9.646 

2  'x2 

15A 

36.942 

26.953 

6.135 

4.719 

6.948 

1  ' 

15B 

35.953 

27 . 851 

6.704 

3.793 

8.715 

30" 

15C 

35.857 

26.876 

6.730 

2.246 

9.694 

2  'x2 

16A 

36.929 

35.092 

6.126 

4.715 

6.991 

1  ' 

16B 

35.833 

36.183 

6.047 

3.651 

7 .886 

2  ' 

16C 

35.941 

35.228 

6.666 

2 . 172 

9.375 

3  'xl 

17A 

36.905 

43.101 

6.062 

4 . 666 

6.979 

1 ' 

17B 

35.890 

44 . 100 

6.152 

3.730 

8 . 038 

2  ' 

17C 

35.895 

43.225 

6.692 

2.208 

9.458 

3  'xl 

18A 

36.701 

51.522 

6.128 

4.717 

6.998 

1 ' 

18B 

35.833 

52.417 

6.130 

3.714 

8 . 047 

2  ' 

18C 

35.707 

51.500 

6.642 

2 . 158 

9.371 

3  'xl 

19A 

47.412 

10.328 

6.092 

4 . 691 

6.952 

1 ' 

19B 

46.562 

11.378 

6.104 

3.698 

8 . 078 

2  ' 

19C 

46.569 

10.241 

6.521 

2 . 125 

9.495 

4  ' 

2  0A 

47.297 

18.387 

6.201 

4 .800 

7 . 087 

1  ' 

20B 

46.471 

19.379 

6.203 

3.792 

8 . 047 

2  ' 

20C 

46.314 

18.344 

6.622 

2.221 

9.523 

4  ' 

21A 

47 .197 

26.743 

6.011 

4.584 

6.886 

1  ' 

21B 

46.331 

27.752 

6.080 

3.669 

7 . 955 

2  ' 

21C 

46.368 

26.757 

6.424 

2 . 034 

9.325 

4  ' 

22A 

47  .111 

35.156 

6.176 

4.770 

7 . 020 

1  ' 

22B 

46.214 

36.240 

6.132 

3.721 

7 . 950 

2  ' 

22C 

46.182 

35.193 

6.676 

2.202 

9.473 

3  'xl 

2  3A 

47.311 

43.143 

6.224 

4 . 802 

7 . 078 

1 ' 

23B 

46.407 

44.362 

6.102 

3.686 

7 . 930 

2  ' 

23C 

46.323 

43.347 

6.728 

2.228 

9.510 

3  'xl 

24A 

47.258 

51.362 

6.048 

4 . 652 

6.913 

1 ' 

24B 

46.315 

52.453 

5.978 

3.572 

7.780 

2  ' 

24C 

46.268 

51.443 

6.584 

2 .090 

9.303 

3  'xl 

2  5A 

57 .867 

10.239 

6.084 

4 . 673 

6.928 

1 ' 

25B 

56.944 

11.205 

6.591 

3.690 

8 . 628 

30" 

25C 

56.951 

10.284 

6.622 

2.226 

9.565 

2  'x2 

Riser 

Type 

S 

S 


Local  Coordinates  (ft) 

Elevations 

(ft) 

In-gr 

Piez  # 

X 

y 

coupling 

screen 

tube 

Lengt 

2  6A 

57 .817 

18.458 

6.169 

4.758 

7 . 013 

1  ' 

2  6B 

56.791 

19.418 

6.652 

3.772 

8 . 668 

30" 

2  6C 

56.832 

18.388 

6.704 

2.230 

9.631 

2  'x2 

27A 

57 . 829 

26.856 

6.117 

4.701 

6.909 

1  ' 

27B 

56.730 

27 . 684 

6.094 

3.693 

7 . 922 

2  ' 

27C 

56.816 

26.827 

6.546 

2 . 150 

9.515 

4  ' 

2  8A 

57 . 932 

35.430 

6.024 

4 . 618 

6.868 

1  ' 

28B 

56.876 

36.484 

5.976 

3.560 

7.768 

2  ' 

28C 

56.863 

35.553 

6.584 

2 .110 

9.366 

3  'xl 

2  9A 

57.796 

43.487 

6.322 

4 . 900 

7 . 176 

1 ' 

29B 

56.660 

44.382 

6.145 

3.744 

8 . 020 

2  ' 

29C 

56.827 

43.576 

7.704 

3.225 

9.460 

3  'xl 

30A 

57.428 

51 . 677 

6.169 

4.737 

6.966 

1 ' 

30B 

56.580 

52.530 

6.128 

3.722 

7 . 946 

2  ' 

30C 

56.574 

51 . 683 

6.700 

2.195 

9.466 

3  'xl 

31A 

68.427 

9.960 

6.016 

4 . 610 

6.865 

1 ' 

31B 

67.506 

11.055 

6.504 

3.614 

8.530 

30" 

31C 

67.354 

10.115 

6.510 

2 . 052 

9.500 

2  'x2 

32A 

68 .161 

18.386 

6.149 

4.722 

6.993 

1  ' 

32B 

67.280 

19.499 

6.582 

3.697 

8 . 619 

30" 

32C 

67.270 

18.503 

6.677 

2.203 

9.615 

2  'x2 

33A 

68 .166 

26.824 

6.152 

4.762 

7 . 048 

1  ' 

33B 

67.346 

27 .865 

6.684 

3.788 

8.549 

30" 

33C 

67.242 

26.981 

6.745 

2.277 

9.745 

2  'x2 

34A 

68.214 

35.694 

6.004 

4.598 

6.838 

1  ' 

34B 

67.257 

36.813 

5.982 

3.566 

7 .889 

2  ' 

34C 

67 . 159 

35.799 

6.446 

1 . 988 

9.280 

3  'xl 

35A 

68.456 

43.356 

6.071 

4 . 634 

6.931 

1 ' 

35B 

67.307 

44.434 

6.050 

3.639 

7 . 821 

2  ' 

35C 

67.403 

43.409 

6.616 

2 . 142 

9.356 

3  'xl 

36A 

68.058 

51.702 

6.141 

4.709 

6.954 

1 ' 

36B 

66.995 

52.534 

6.107 

3.691 

7 . 956 

2  ' 

36C 

67.051 

51 . 689 

6.669 

2.206 

9.497 

3  'xl 

37A 

78 . 677 

10.092 

6.103 

4 . 676 

6.921 

1 ' 

37B 

77 . 622 

11.099 

6.634 

3.744 

8.718 

30" 

37C 

77 . 669 

10.015 

6.683 

2 .189 

9.569 

2  'x2 

38A 

78.715 

18.647 

6.145 

4.734 

6.963 

1  ' 

38B 

77.497 

19.554 

6.702 

3.780 

8.702 

30" 

38C 

77 . 612 

18.739 

6.716 

2.258 

9.706 

2  'x2 

Riser 

Type 

S 

S 


Local  Coordinates  (ft) 

Elevations 

(ft) 

In-gr 

Piez  # 

X 

y 

coupling 

screen 

tube 

Lengt 

39A 

78.744 

26. 80S 

6. 14S 

4.749 

6.968 

1  ' 

39B 

77.792 

27 . 828 

6.702 

3.801 

8.739 

30" 

39C 

77 . 840 

26.789 

6.700 

2.216 

9.622 

2  'x2 

4  0A 

78 . 678 

3S.129 

6.074 

4 . 663 

6.892 

1  ' 

40B 

77.S61 

36.118 

6.0S4 

3.669 

7.4SS 

2  ' 

40C 

77.702 

3S.039 

6.S99 

2  .  IIS 

9.36S 

3  'xl 

41A 

78.S00 

43.090 

6.  lOS 

4 . 678 

6.933 

1 ' 

41B 

77.727 

44.080 

6. 08S 

3.663 

7 .867 

2  ' 

41C 

77.S42 

43.141 

6.717 

2.233 

9.446 

3  'xl 

42A 

78.S90 

SI . 624 

6.064 

4 . 6S8 

6. 6S3 

1 ' 

42B 

77 . 647 

S2 . 600 

6.0S4 

3.643 

7 . 9SS 

2  ' 

42C 

11 .518 

SI . 69S 

6.667 

2.199 

9.407 

3  'xl 

4  3A 

89.17S 

10.163 

6.048 

4 . 632 

6.840 

1 ' 

43B 

88.211 

11 . 134 

6.SS0 

3.649 

8.S82 

30" 

43C 

88.278 

10.318 

6.S90 

2 . 106 

9.SS9 

2  'x2 

44A 

89.181 

18.737 

6.162 

4.740 

6.824 

1  ' 

44B 

88  .  ISS 

19.840 

6. 6S1 

3.740 

8 . 667 

30" 

44C 

88 . 132 

18 . 827 

6.719 

2.230 

9.6S2 

2  'x2 

4  5A 

89.029 

26. 666 

6. OSS 

4 . 649 

6.873 

1  ' 

45B 

88.076 

27.717 

6.624 

3.702 

8 . 63S 

30" 

45C 

88.024 

26.834 

6.662 

2 .183 

9.S74 

2  'x2 

4  6A 

89.107 

34 . 983 

6.136 

4.720 

6.918 

1  ' 

4  6B 

88 . 148 

3S.984 

6.104 

3.688 

7 . 990 

2  ' 

4  6C 

88 . 173 

3S.044 

6.698 

2.230 

9.448 

3  'xl 

47A 

89.0S7 

43.2S1 

6. 1S4 

4.732 

7  .  OSS 

1 ' 

47B 

87 . 964 

44.413 

6.049 

3.648 

7 . 924 

2  ' 

47C 

87 . 978 

43.319 

6.736 

2.226 

9.S96 

3  'Xl 

4  8A 

89.176 

SI  .800 

6.112 

4.701 

6.784 

1  ' 

48B 

88.089 

S2 . 64S 

6. 12S 

3.724 

7 . 98S 

2  ' 

48C 

88.033 

SI. 760 

6.630 

2 . 1S6 

9.417 

3  'Xl 

4  9A 

99.776 

10.3S9 

S.942 

4.S26 

6.7SS 

1  ' 

49B 

98 . 6S2 

11.331 

6.S80 

3.664 

8 . 601 

30" 

49C 

98 .891 

10.317 

6.636 

2 . 1S2 

9.647 

2  'x2 

50A 

99.672 

18 . 6S7 

6. 1S2 

4.741 

6.928 

1  ' 

SOB 

98.6S0 

19.803 

6.106 

3.690 

7 . 9S0 

2  ' 

50C 

98.660 

18.791 

6.679 

2 .20S 

9.69S 

2  'x2 

51A 

99.634 

26.672 

6. 1S6 

4.7S0 

6.984 

1  ' 

51B 

98.739 

27 . 612 

6.192 

3.791 

8 . 088 

2  ' 

SIC 

98.711 

26.733 

6.731 

2 .2S2 

9.6S8 

2  'x2 

Riser 

Type 

S 

S 


Local  Coordinates  (ft) 

Elevations 

(ft) 

Piez  # 

X 

y 

coupling 

screen 

tube 

52A 

99.G13 

35.029 

G.228 

4.79G 

G.  890 

52B 

98.578 

3G. 144 

G.210 

3.778 

8 . 09G 

52C 

98 . GG3 

35.119 

G.717 

2.259 

9.551 

53A 

99.G57 

43.225 

G.  172 

4.77G 

G.  933 

53B 

98.543 

44.323 

G.043 

3.G42 

7 . 929 

53C 

98.575 

43.354 

G.7G3 

2.248 

9.451 

54A 

99.59G 

51.771 

G.043 

4  .  G42 

G.  G89 

54B 

98.550 

52 . 821 

G.  148 

3.732 

7 . 9GG 

54C 

98.572 

51 . 815 

G.  G71 

2.192 

9.499 

55A 

110.178 

10.128 

G.030 

4  .  G24 

G.  G71 

55B 

109.229 

11.250 

G.  153 

3.72G 

8 . 122 

55C 

109. 2G8 

10.101 

G.  G30 

2 . 14G 

9.5G8 

56A 

110. 23G 

18.435 

G.0G9 

4  .  G73 

G.710 

56B 

109.282 

19.GG7 

G.  117 

3.721 

8 . 018 

56C 

109. 32G 

18.5GG 

G.  G73 

2.210 

9.G21 

57A 

110.057 

2G.71G 

G.  109 

4.708 

G.7G0 

57B 

109.154 

27 .804 

G.  192 

3.770 

7.989 

57C 

109.052 

2G.  823 

G.725 

2.241 

9.757 

58A 

110. 1G3 

35.170 

G.  108 

4.707 

G.733 

58B 

108 . 920 

3G. 192 

G.  1G8 

3.757 

8 . 012 

58C 

109. OGO 

35.214 

G.709 

2.220 

9.449 

59A 

109.990 

42 . 997 

G.  13G 

4.730 

G.782 

59B 

108 . 991 

44.258 

G.  175 

3.748 

8 . 035 

59C 

109.000 

43.153 

G.735 

2.251 

9.5G9 

GOA 

110.135 

51 . 133 

G.  172 

4.771 

G.  8G0 

GOB 

109.137 

52.304 

G.  149 

3.727 

7 . 972 

GOC 

109.093 

51.249 

G.  G72 

2 .198 

9.459 

GIA 

120.532 

9.941 

5.922 

4.521 

G.579 

GIB 

119.570 

10.94G 

G.0G5 

3.G43 

7 . 951 

GIC 

119. 50G 

9.959 

G.597 

2.212 

9.57G 

G2A 

120. G9G 

18.438 

G.OGl 

4  .  GG5 

G.  994 

G2B 

119.707 

19.402 

G.  080 

3.G38 

7 .997 

G2C 

119. GG9 

18.490 

G.  GIO 

2 . 12G 

9.548 

G3A 

120.714 

2G.71G 

G.055 

4  .  G49 

G.  878 

G3B 

119.790 

27.749 

G.089 

3.G78 

7 . 975 

G3C 

119.792 

2G.793 

G.713 

2.208 

9.713 

G4A 

120. G17 

35.092 

G.  IGG 

4.7G5 

7 . 047 

G4B 

119.475 

3G. 183 

G.22G 

3.794 

8  .  OGO 

G4C 

119. G79 

35.152 

G.  G99 

2.215 

9.49G 

In-ground  Riser 

Length  Type 

S 

2  '  S 

3'Xl'  S 

S 

2  '  S 

3'Xl'  S 

1 '  G 

2  '  S 

3'Xl'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

3'Xl'  S 

1 '  G 

2  '  S 

3'Xl'  S 

1 '  G 

2  '  S 

3'Xl'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

2'x2'  S 

1 '  G 

2  '  S 

3'Xl'  S 


Piez  # 

Local  Coordinates  (ft) 

Elevations  (ft) 

In-ground  Riser 

X 

y 

coupling 

screen 

tube 

Length 

Type 

65A 

120.688 

43.250 

6.127 

4.726 

6.992 

1  ' 

G 

65B 

119.620 

44.427 

6.104 

3.708 

7 . 990 

2  ' 

S 

65C 

119.632 

43.212 

6.779 

2.279 

9.581 

3  'XI  ' 

S 

66A 

120.530 

51.529 

6.049 

4 . 659 

6.867 

1  ' 

G 

66B 

119.612 

52.427 

6.113 

3.717 

8 . 046 

2  ' 

S 

66C 

119.656 

51.456 

6.682 

2.286 

9.641 

2  'x2  ' 

S 

SI 

0.000 

0.000 

6.572 

S2 

134 . 648 

0.000 

7.589 

S3 

0.003 

67.409 

8.357 

S4 

134 . 616 

67.551 

9.256 

WEIR 

127 . 682 

29.790 

7.190 

Notes  : 

(1)  The  conversion  to  local  coordinates  is  based  on  a  survey  conducted 
12-05-2001  using  Total  Station  Survey  equipment.  Construction  stake  SI 
forms  the  origin  for  the  coordinate  system.  Elevations  are  relative  to 
the  approximate  elevation  of  the  impermeable  wetland  liner,  assumed  to 
be  at  817.6097'  mean  sea  level. 

(2)  "Coupling"  elevation  is  the  elevation  of  the  above  ground  coupling; 
"screen  elevation"  is  the  elevation  at  the  midpoint  of  the  piezometer 
screen;  "tube  elevation"  is  the  elevation  at  the  top  of  the  exposed 
1/2"  Teflon  tube. 

(3)  "S"  indicates  stainless  steel  is  used  for  the  riser;  "G"  indicates 
galvanized  steel  is  used  for  the  riser. 

(4)  Piezometer  61B  is  severely  canted.  Piezometer  6C  has  never 
registered  a  water  level.  Piezometer  29C  hit  a  rock  during 
installation  and  is  positioned  in  the  middle  soil  layer. 
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Appendix  B:  Contours  of  Hydraulic  Head  (11-01-2001) 


(Vertical  Exaggeration: 


Profile  View  of  Equipotentials :  11-01-2001 


Profile  View  of  Equipotentials :  02-19-2002 
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Profile  View  of  Equipotentials :  02-19-2002 


Piezometers:  4,  10,  16,  22,  28,  34,  40,  46,  52,  58,  64 
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Appendix  D:  Water  Level  Measurements 


Piez  # 
OlA 
OIB 
QIC 
02A 
02B 
02C 
03A 
03B 
03C 
04A 
04B 
04C 
05A 
05B 
05C 
0  6A 
06B 
06C 
07A 
07B 
07C 
08A 
08B 
08C 
09A 
09B 
09C 
lOA 
lOB 
IOC 
llA 
IIB 


3/17/0l|l0/18/0l|l0/19/0l|l0/22/0l|l0/2  9/0l| 

11/1/01 

11/8/01 

2/19/02 

0.850 

0.125 

0.219 

0.292 

0.810 

0.980 

1.050 

1.040 

2.490 

2.200 

2.270 

2.380 

2.590 

2 . 610 

2.700 

2.780 

1.520 

1.530 

1.500 

1.510 

1 . 600 

1.580 

1.515 

1 . 920 

1.410 

1.300 

1.190 

0.990 

1 . 170 

1.080 

1 . 620 

0.950 

0.580 

0.640 

0.610 

0.530 

0.855 

0.580 

1.235 

1.380 

1.370 

1.360 

1.330 

1.294 

1.390 

1.370 

1.320 

1.770 

0.510 

0.520 

0.500 

0.443 

0.645 

0.550 

1.278 

0.945 

1.307 

1.220 

1 . 075 

0.980 

1.360 

1.090 

1 . 150 

1.455 

1.350 

1.350 

1.310 

1.273 

1.370 

1.360 

1.298 

1.755 

0.600 

0.570 

0.530 

0.469 

0.680 

0.575 

0.570 

1 . 120 

0.960 

0.880 

0.860 

0.820 

1.010 

0.970 

0.930 

1 .185 

1.720 

1 . 680 

1 . 630 

1 . 620 

1.770 

1.710 

1 . 640 

1.720 

0.990 

0.960 

0.890 

0.820 

1 . 050 

0.910 

0.920 

0.915 

1 .  930 

1.580 

1 .800 

1.770 

1 .  985 

1.890 

1 .  960 

2.050 

1.370 

1.350 

1.315 

1.268 

1.370 

1.360 

1.304 

1 .800 

1.260 

1.240 

1 . 170 

1 . 060 

1.320 

1 . 170 

1.304 

1 .190 

2.200 

2 . 150 

2.070 

2.020 

2.210 

2.095 

2.075 

2 . 140 

0.700 

0.590 

0.500 

0.440 

0.745 

0.670 

0.770 

0.910 

1.580 

1.510 

1.440 

1.380 

1 . 635 

1.510 

1.590 

1 . 855 

1.400 

1.380 

1.370 

1.350 

1.420 

1.400 

1.360 

1 .813 

1.025 

0.940 

0.795 

0.700 

1.010 

0.810 

0.875 

0.960 

1 .  910 

1.780 

1 . 650 

1.565 

1 .885 

1 . 680 

1.735 

1 .885 

1.390 

1.370 

1.360 

1.350 

1.415 

1.400 

1.350 

1 .800 

0.430 

0.370 

0.360 

0.470 

0.710 

0.690 

0.635 

0.810 

1.460 

1.410 

1.385 

1.465 

1.705 

1.595 

2 . 145 

1 .  980 

1.490 

1.480 

1.450 

1.435 

1.510 

1.490 

1.445 

1 .  920 

0.525 

0.490 

0.480 

0.500 

0.710 

0.700 

0.725 

0.840 

1.570 

1.520 

1.520 

1.510 

1.720 

1.730 

1 . 820 

1.780 

1.525 

1.500 

1.490 

1.470 

1.540 

1.530 

1.480 

1 .  940 

1.030 

1 . 020 

0.975 

0.855 

1 . 130 

0.930 

1 . 615 

0.935 

2.220 

2.220 

2 .185 

2.080 

2.330 

2 . 140 

2 . 130 

Piez  # 
lie 
12A 
12B 
12C 
13A 
13B 
13C 
14A 
14B 
14C 
15A 
15B 
15C 
16A 
16B 
16C 
17A 
17B 
17C 
18A 
18B 
18C 
19A 
19B 
19C 
2  0A 
20B 
20C 
21A 
21B 
21C 
22A 
22B 
22C 
2  3A 


3/17/0l|l0/18/0l|l0/19/0l|l0/22/0l|l0/2  9/0l| 

11/1/01 

11/8/01 

1.420 

1.400 

1.385 

1.370 

1.430 

1.420 

1.370 

0.795 

0.910 

0.910 

0.770 

1 .180 

0.790 

1 . 170 

1.700 

1.790 

1.700 

1.520 

2.030 

1.490 

1.580 

1.410 

1.400 

1.380 

1.365 

1.440 

1.430 

1.380 

1 .110 

1 . 030 

0.900 

0.755 

1 . 135 

0.870 

1.000 

1.400 

3.210 

2 . 150 

2.030 

2.490 

2.090 

2.220 

1 . 610 

1.570 

1.560 

1.550 

1 . 620 

1 . 600 

1.560 

0.840 

0.900 

0.840 

0.670 

1 . 120 

0.715 

0.820 

2.290 

2.350 

2.200 

2.020 

2.580 

2.050 

2 . 140 

1.580 

1.550 

1.550 

1.530 

1 . 600 

1.580 

1.535 

0.600 

0.450 

0.380 

0.370 

0.710 

0.535 

0.479 

2.290 

2.190 

2.090 

2.080 

2.440 

2.210 

2 .180 

1 . 630 

1 . 600 

1 . 600 

1.590 

1 . 650 

1 . 640 

1 . 600 

0.445 

0.385 

0.375 

0.417 

0.620 

0.540 

0.417 

1.310 

1.290 

1.325 

1 .190 

1.495 

1.289 

1.220 

1.318 

1.310 

1.283 

1.236 

1.340 

1.312 

1.268 

0.980 

0.980 

0.810 

0.565 

1.215 

0.610 

0.730 

2.210 

2 . 130 

1 . 840 

1.595 

2.240 

1 . 630 

1.745 

1.390 

1.350 

1.360 

1.350 

1.420 

1.400 

1.350 

1.770 

1.560 

1.340 

0.995 

1 .180 

1 . 035 

2 . 980 

1 . 910 

1.710 

1 . 660 

2 . 050 

1.780 

1 . 855 

1.318 

1.270 

1.273 

1.242 

1.320 

1.315 

1.256 

0.990 

0.930 

0.920 

0.900 

1 . 080 

1.045 

1.060 

1 .  930 

1 . 840 

1.725 

1 . 680 

2 . 000 

1 . 850 

1 . 875 

1.500 

1.470 

1.490 

1.460 

1.540 

1.500 

1.470 

1.580 

1.580 

1.420 

1 .160 

1 . 845 

1.225 

1.345 

1.460 

1.410 

1.430 

1.410 

1.480 

1.450 

1.410 

0.420 

0.430 

0.400 

0.203 

0.620 

0.443 

1.400 

1.470 

1.435 

1.290 

1 .185 

1 . 650 

1.278 

1 . 120 

1.245 

1.220 

1.225 

1.210 

1.772 

1.269 

1.210 

0.188 

0.172 

0.177 

0.198 

0.420 

0.417 

0.182 

0.770 

0.760 

0.750 

0.795 

1 . 040 

1.050 

0.510 

1.420 

1.370 

1.380 

1.360 

1.430 

1.420 

1.370 

0.560 

0.520 

0.530 

0.545 

0.820 

0.780 

0.720 

2/19/02 
1 . 870 
0.780 
1 . 670 
1 . 850 
1 . 170 
2.265 
2.005 
0.960 
2.485 
2 . 005 
0.870 
2.475 
2.068 
0.820 
1.730 
1.740 
0.945 
1 . 960 
1.849 
1.200 
2.040 
1.780 
1.060 
1.780 
1 . 920 
0.935 
1.520 
1 . 875 
0.850 
1.720 
1 . 680 
1.200 
1 .  950 
1.755 
1.030 


Piez  # 
23B 
23C 
24A 
24B 
24C 
2  5A 
25B 
25C 
2  6A 
2  6B 
2  6C 
27A 
27B 
27C 
2  8A 
28B 
28C 
2  9A 
29B 
29C 
30A 
30B 
30C 
31A 
31B 
31C 
32A 
32B 
32C 
33A 
33B 
33C 
34A 
34B 
34C 


10/17/0110/18/0110/19/0110/22/0110/29/01  11/1/01 
1.276  1.240  1.170  1.130  1.435  1.200 

1.450  1.410  1.420  1.410  1.480  1.470 

1.050  0.960  0.810  0.650  1.000  0.700 

1.150  1.140  1.150  1.145  1.729  1.180 

1.250  1.210  1.220  1.200  1.623  1.250 

0.520  0.490  0.430  0.375  0.750  0.890 

2.330  2.320  2.240  2.020  2.555  2.100 

1.800  1.760  1.700  1.695  1.810  1.730 

2.060  2.190  2.080  1.810  2.470  1.960 

1.570  1.530  1.545  1.520  1.590  1.580 


11/8/01 

1.220 

1.410 

0.875 

2.445 

1.200 

0.840 

2.345 

1.730 

2.020 

1.585 


3.470 

3.220 

2 . 960 

2.290 

1.500 

1.380 

1.440 

1.410 

1.415 

1.400 

1.460 

1.440 

1.400 

0.266 

0.271 

0.203 

0.047 

0.390 

0.167 

0.375 

1.281 

1.210 

1 . 020 

0.950 

1 . 835 

0.970 

0.052 

1.297 

1.260 

1.268 

1.216 

1.320 

1.309 

1.256 

0.208 

0.193 

0.255 

0.240 

0.375 

0.417 

0.156 

0.830 

0.850 

0.860 

0.840 

0.975 

0.940 

0.930 

2 . 135 

2.090 

2 . 000 

1 . 970 

2 .110 

2 . 000 

1 .890 

0.182 

0.151 

0.115 

0.130 

0.328 

0.333 

0.313 

0.610 

0.570 

0.550 

0.545 

0.700 

0.630 

0.650 

1.410 

1.380 

1.380 

1.370 

1.450 

1.420 

1.390 

0.560 

0.420 

0.156 

0.042 

0.510 

0.156 

0.229 

2.210 

2 . 140 

1 . 930 

1.770 

2.230 

1 . 840 

1 .  910 

1.440 

1.400 

1.400 

1.385 

1.460 

1.445 

1.410 

0.890 

0.850 

0.800 

0.905 

0.980 

0.830 

0.820 

2 . 150 

2 . 125 

2.020 

1 . 875 

2.290 

1 . 975 

2.040 

1.540 

1.510 

1.515 

1.500 

1.570 

1.550 

1.520 

0.750 

0.385 

0.198 

0.830 

0.840 

0.835 

0.980 

1 . 000 

0.745 

1.200 

1 . 670 

1 . 630 

1 . 650 

1 . 640 

1.700 

1 . 680 

1 . 645 

0.670 

0.550 

0.490 

0.453 

0.680 

0.610 

0.555 

0.318 

0.292 

0.240 

0.188 

0.339 

0.229 

0.161 

1.210 

1 .185 

1.195 

1.570 

1.245 

1.230 

1 .195 

2/19/02 
1.560 
1 . 870 
0.950 
1.230 
1 . 660 
1.045 
2.330 
1 . 900 
1.250 
2 .800 
2.036 
1 .100 
1.210 
1 .  910 
0.520 
2.000 
1.730 
0.950 
1.570 
2.510 
0.745 
1.283 
1 . 823 
1.040 
2 . 650 
1 . 830 
1 . 170 
2.580 
2.000 
0.594 
2 . 620 
2  .  Ill 
0.910 
0.698 
1 . 640 


Piez  # 
35A 
35B 
35C 
36A 
36B 
36C 
37A 
37B 
37C 
38A 
38B 
38C 
39A 
39B 
39C 
4  0A 
40B 
40C 
41A 
41B 
41C 
42A 
42B 
42C 
4  3A 
43B 
43C 
44A 
44B 
44C 
4  5A 
45B 
45C 
4  6A 
46B 


3/17/0l|l0/18/0l|l0/19/0l|l0/22/0l|l0/2  9/0l| 

11/1/01 

11/8/01 

0.610 

0.560 

0.580 

0.560 

0.770 

0.760 

0.730 

0.370 

0.348 

0.292 

0.281 

0.318 

0.188 

0.057 

1.271 

1.260 

1.252 

1.250 

1.320 

1.304 

1.266 

0.151 

0.073 

0.026 

0.000 

0.188 

0.000 

-0.135 

0.835 

0.810 

0.755 

0.705 

0.790 

0.700 

0.665 

1.445 

1.410 

1.420 

1.405 

1.480 

1.460 

1.420 

0.980 

0.900 

0.905 

0.880 

1 . 070 

1.050 

1.040 

2.370 

2.300 

2 . 170 

2 . 120 

2.565 

2.270 

2.370 

1.500 

1.475 

1.480 

1.465 

1.535 

1.510 

1.480 

0.385 

0.346 

0.276 

0.109 

0.302 

0.219 

0.980 

1 .  960 

1 . 940 

1 . 870 

1.785 

2.200 

2.050 

2.500 

1 . 640 

1 . 610 

1 . 630 

1 . 620 

1 . 680 

1 . 660 

1 . 630 

1.330 

1 . 140 

1.580 

2 . 040 

2.070 

1.750 

2 . 105 

1 . 640 

1.200 

1.550 

1.520 

1.530 

1.560 

1.585 

1.565 

1.530 

0.980 

1.055 

0.980 

0.720 

1 . 135 

0.880 

0.890 

1.770 

1.560 

1.520 

1.250 

1 . 600 

1.360 

1.390 

1.305 

1.290 

1.290 

1.320 

1.350 

1.320 

1.283 

1.785 

1.320 

1 .  915 

1.550 

1.530 

1.520 

1.420 

1.590 

1.410 

1.595 

1.380 

1.355 

1.365 

1.380 

1.420 

1.400 

1.365 

0.760 

0.620 

0.470 

0.286 

0.710 

0.490 

0.515 

1 . 120 

1 .110 

1.010 

0.895 

1.220 

0.950 

1.010 

1.330 

1.315 

1.304 

1.320 

1.380 

1.350 

1.320 

1.495 

1 . 940 

2.080 

2.080 

1 . 830 

2.360 

1 . 985 

2.775 

1.490 

1.460 

1.480 

1.480 

1.530 

1.510 

1.470 

2 . 135 

2 .180 

2.200 

2 .165 

2 .160 

2.205 

1.570 

1.540 

1.550 

1.555 

1 . 610 

1.590 

1.555 

1 . 000 

0.950 

0.890 

0.780 

1 .110 

0.910 

1 .100 

2 .810 

2.700 

2.570 

2.420 

2 .810 

2.515 

2 . 630 

1.500 

1.460 

1.475 

1.470 

1.530 

1.520 

1.480 

1.275 

1.263 

1.065 

0.810 

1.370 

0.755 

1.215 

2.500 

2.440 

2 . 150 

1 . 615 

2.450 

1.730 

2.080 

2/19/02 
0.970 
0.860 
1.700 
1.070 
1.400 
1.755 
0.850 
1.740 
1 .990 
0.635 
2 . 600 
2.075 
1 . 150 
2 .880 
2.010 
0.930 
1.450 
1.750 
1.050 
2 . 010 
1 .802 
0.670 
1 . 610 
1.755 
1.070 
2 .160 
1 .995 
0.542 
2.240 
2.040 
1 .180 
2.720 
1 .  970 
1.025 
1 .  960 


Piez  # 
4  6C 
47A 
47B 
47C 
4  8A 
48B 
48C 
4  9A 
49B 
49C 
50A 
SOB 
50C 
51A 
51B 
SIC 
S2A 
S2B 
S2C 
S3A 
S3B 
S3C 
S4A 
S4B 
S4C 
SSA 
SSB 
SSC 
S6A 
S6B 
S6C 
S7A 
S7B 
S7C 
S8A 


10/17/01 10/18/01 10/19/01 10/22 /Oil  0/2 9/0111/1/01  11/8/01 


1.390 

1.3S0 

1.370 

1.360 

1.430 

1.39S 

1.36S 

1.400 

1.4S0 

1.370 

0.900 

1.S90 

1.030 

2.260 

2.370 

2 .26S 

1.780 

2.SS0 

1 . 840 

2.390 

1.S20 

1.480 

1.490 

1.490 

l.SSO 

1.S30 

1.490 

1.330 

1.220 

1 . 100 

0.910 

1 . 080 

0.940 

1.294 

3.170 

3.070 

2.780 

2.030 

1.430 

1.320 

1.360 

1.320 

1.330 

1.283 

1.380 

1.370 

1.32S 

0.89S 

0.920 

0.820 

0.690 

1 . 08S 

0.82S 

1 .180 

2.370 

2.S20 

2.S40 

2.410 

2 . 87S 

2.490 

3.0S0 

1.S80 

l.SSO 

l.SSO 

1.S40 

1 . 610 

1.S90 

1.S60 

0.970 

0.960 

0.9SS 

0. 9S0 

1 .100 

1 .100 

1 . 100 

3.910 


1 . 630 

1 . 600 

1 . 60S 

1.S90 

1 . 66S 

1 . 63S 

1 . 610 

1.430 

1.294 

1 . 170 

1.S60 

1.070 

2.210 

2 . 040 

1 . 94S 

1 . 820 

2 . 140 

2.000 

2.030 

1.S80 

l.SSS 

1.S70 

1.S40 

1 . 62S 

1 . 610 

1.S70 

0.810 

0.770 

0.680 

1.273 

0.800 

0.600 

1.340 

2.710 

2.190 

2.0S0 

1.490 

1.470 

1.470 

1.4S0 

1.S30 

l.SOO 

1.470 

0.870 

0.800 

0.770 

0.690 

0.9S0 

0.890 

0.880 

1 . 940 

1 . 87S 

1 . 840 

1.7S0 

2 . 030 

1 . 9SS 

1 .  9SS 

1.380 

1.3S0 

1.360 

1.33S 

1.420 

1.39S 

1.360 

0.860 

0.670 

0.S20 

0.286 

0.76S 

0.S90 

0.640 

1.010 

1.080 

l.OSO 

0.870 

1.708 

0.990 

2.200 

1.43S 

1.410 

1.41S 

1.400 

1.470 

1.4S0 

1.41S 

0.710 

0.660 

0.680 

0.660 

0.830 

0.83S 

0.810 

1.470 

1.S90 

1 . 640 

1 . 600 

1 .  920 

1.72S 

3.S70 

1.490 

1.460 

1.480 

1.46S 

1.S20 

I.SIO 

1.470 

1.S40 

l.SOO 

1.S30 

l.SOO 

1.S80 

1.S60 

1.S20 

0.990 

0.890 

0.830 

0.760 

0.930 

0. 90S 

2 .110 

2.020 

1 . 98S 

1 . 930 

2 .180 

2 .110 

2  .  lOS 

1 . 690 

1 . 660 

1 . 670 

1 . 6S0 

1.720 

1.700 

1 . 66S 

0.60S 

O.SIS 

0.440 

0.344 

0.S60 

O.SOO 

O.SOO 

2/19/02 
1 .880 
1.060 
1.70S 
1 . 99S 
0.960 
1 . 63S 
1.781 
0.910 
2.710 
2.021 
1.07S 

1 .  910 
2.060 
0.980 
2.000 

2.  QOS 
0.800 
2.000 
2.000 
0.920 
1 . 8S0 
1 .  900 
0.830 
1 . 130 
1 . 917 
0.780 
1 .  920 
1 . 9S0 
0.920 
2.230 
1 . 960 
0.940 
2  .  ISO 
2  .  ISO 
0.740 


Piez  #10/17/0110/18/0110/19/0110/22/0110/29/0111/1/01  11/8/01  2/19/02 


58B 

1 .880 

1.780 

1.720 

1 .  G50 

1 . 830 

1.7G5 

1.7G0 

1 .815 

58C 

1.390 

1.3G0 

1.370 

1.350 

1.420 

1.400 

1.3G0 

1 . 870 

5  9A 

1.000 

0.910 

0.730 

0.4G9 

0.870 

O.GIO 

1.410 

0.800 

59B 

2 . 050 

2 . 000 

1 .  900 

1.740 

2 . 080 

1 . 875 

2.385 

2.0G0 

59C 

1.500 

1.470 

1.490 

1.4G5 

1.530 

1.520 

1.480 

2.020 

GOA 

0.302 

1.570 

0.400 

0.000 

0.28G 

0.052 

0.115 

0.730 

GOB 

1 . 170 

1 .100 

1.380 

1.330 

1.400 

0.945 

1.210 

2.030 

GOC 

1.400 

1.390 

1.380 

1.370 

1.440 

1.420 

1.385 

1 .885 

GIA 

0.417 

0.400 

0.335 

0.094 

0.570 

0.292 

O.GIO 

GIB 

1.790 

1.700 

1 . 830 

1.785 

2 . 005 

1.710 

2 . 870 

1.700 

GIC 

1.570 

1.550 

1.570 

1.5G0 

1 .  G40 

1 .  G20 

1.590 

2.070 

G2A 

1 .100 

G2B 

1.080 

1.200 

1.190 

0.920 

1.580 

l.OGO 

1.5G0 

1.435 

G2C 

1.470 

1.430 

1.445 

1.435 

1.500 

1.480 

1.445 

1 .  953 

G3A 

0.G70 

0.G70 

0.595 

0.450 

0.830 

0.570 

0.830 

1.050 

G3B 

0.970 

1.050 

1.015 

0.870 

1.335 

1.000 

1.500 

1.415 

G3C 

1 .  G40 

1 .  G15 

1 .  G30 

1 .  GIO 

1 .  G80 

1 .  GGO 

1 .  G15 

2.099 

G4A 

0.920 

0.8G5 

0.800 

0.G70 

0.990 

0.820 

0.880 

0.990 

G4B 

0.830 

0.780 

0.G80 

0.595 

0.930 

0.700 

0.735 

1.020 

G4C 

1.440 

1.405 

1.415 

1.400 

1.4G0 

1.440 

1.415 

1 . 938 

G5A 

0.1G7 

0.1G2 

0.104 

0.135 

0.328 

0.281 

0.339 

0.G5G 

G5B 

0.354 

0.345 

0.297 

0.28G 

0.427 

0.417 

0.380 

0.830 

G5C 

1.520 

1.485 

1.500 

1.480 

1.550 

1.530 

1.510 

1 . 990 

GGA 

0.078 

0.078 

0.021 

O.OIG 

0.120 

0.042 

0.3G5 

0.531 

GGB 

0.G05 

0.890 

0.595 

O.GIO 

0.705 

0.G70 

0.720 

1 . 170 

GGC 

1.580 

1.530 

1.550 

1.540 

1 .  GIO 

1.595 

1.5G0 

2.020 

water 

table 

1.208 

1.214 

1.214 

1.214 

1.370 

1.375 

1.370 

1.391 

Notes  : 

(1)  All  measurements  in  feet. 

(2)  Distance  to  the  water  level  in  the  piezometers  is  measured  from  the 
top  of  the  W  plastic  piezometer  tube. 

(3)  Distance  to  the  water  table  is  measured  from  the  south  edge  of  the 
horizontal  metal  bar  on  the  weir. 
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Appendix  E:  Bouwer  and  Rice  Slug  Test  Results 


Top  Soil  Layer 


Piez 

Test  Date 

Early  K 

Middle  K 

Overall  K 

OlA 

12/07/01 

1.54E-05 

7 . 9  6E- 0  6 

7 .22E-06 

0.92 

OlA 

12/07/01 

1 . 65E-05 

1 . 16E-05 

1.04E-05 

0.97 

02A 

12/07/01 

1.31E-05 

1.00 

03A 

12/07/01 

too  fast 

04A 

12/07/01 

3.66E-05 

0.99 

05A 

12/07/01 

8.37E-06 

1.50E-06 

1 . llE-06 

0.94 

0  6A 

12/07/01 

1.53E-05 

1.00 

07A 

12/07/01 

1.47E-06 

0.98 

08A 

0  9A 

11/08/01 

1 . 68E-08 

1.00 

09A 

12/04/01 

3.45E-07 

1.00 

lOA 

llA 

11/16/01 

2 .28E-06 

0.99 

12A 

13A 

12/07/01 

1 .22E-05 

1 . 90E-06 

1.57E-06 

0.82 

14A 

12/11/01 

3.32E-06 

0.99 

15A 

11/08/01 

1.59E-07 

0.99 

15A 

11/16/01 

1 .22E-07 

0.99 

15A 

12/15/01 

2.43E-06 

1.09E-06 

7.72E-07 

0.90 

16A 

12/07/01 

3.22E-05 

1.06E-05 

1 .27E-05 

0.90 

17A 

18A 

12/07/01 

1 .20E-05 

6 . 92E- 0 6 

4.77E-06 

0.97 

19A 

12/07/01 

4 . 60E-05 

1.77E-05 

1 . 91E-05 

0.94 

2  0A 

21A 

12/07/01 

1 . 16E-05 

1 . 60E-06 

1.58E-06 

0.94 

22A 

2  3A 

11/14/01 

1 .27E-07 

1.00 

23  A 

12/11/01 

1.41E-05 

0.99 

23  A 

12/11/01 

3.06E-05 

1.00 

24A 

2  5A 

12/04/01 

6.03E-06 

3.59E-06 

2 . 9 6E- 0 6 

0.97 

2  6A 

12/07/01 

1 . 84E-06 

7 .28E-07 

2 . 61E-07 

0.72 

27A 

12/04/01 

5.13E-06 

8.56E-07 

6. 16E-07 

0.82 

2  8A 

12/04/01 

4 . 98E-06 

1 . 15E-06 

0.97 

29  A 

12/04/01 

too  fast 

Type  Recov 


74.4% 
68 . 9% 


Test  Date  Early  K  Middle  K  Overall  K 
12/11/01  6.71E-05  3.27E-05 
11/02/01  8.56E-08 
12/11/01  1.33E-06 
12/11/01  1.59E-05 
12/11/01  1.50E-05 
11/02/01  3.92E-09 
11/21/01  1.02E-08 
12/15/01  2.30E-06 
11/02/01  5.82E-08 
11/16/01  9.15E-08 
12/07/01  9.91E-06  2.74E-06 
12/11/01  1.72E-05  l.OOE-05 
11/02/01  1.04E-07 
12/11/01  1.28E-05  7.08E-06  8.04E-06 


Recov 


10/19/01 

11/16/01 

12/15/01 

11/02/01 

11/21/01 

12/04/01 


5.56E-08 
5.14E-08 
too  fast 
9. 92E-09 
3.87E-06 
2.07E-06 


12/04/01  6.68E-06  4.07E-06  2.97E-06 


11/21/01  4.77E-06  2.62E-06  1.85E-06  0.96 


12/05/01  6.50E-06  2.28E-06  1.51E-06  0.90 


12/05/01 

12/05/01 


5.68E-08  0.98 

4.92E-05  1.00 


12/05/01  5.23E-06  1.24E-06  9.42E-07  0.98 

12/05/01  4.52E-06  1.16E-06  3.21E-07  0.60 

12/05/01  1.05E-05  2.03E-07  5.25E-07  0.88 


12/11/01 


11/16/01 


2.49E-05  1.00 


2.42E-06  0.98 


12/11/01  1.15E-05  5.58E-06  1.96E-06  0.79 


12/11/01  4.63E-05 


12/11/01  3.70E-05 


2.84E-05  0.97 


1.18E-05  0.89 


Piez 

Test  Date 

Early  K 

Middle  K 

Overall  K 

R^ 

Dev 

Type 

Recov 

58A 

59A 

12/11/01 

7 . 94E-06 

2 . 63E-06 

2 . 06E-06 

0.93 

2 

+ 

82 . 6% 

GOA 

11/20/01 

2 . 66E-08 

0.99 

- 

- 

22.7% 

61A 

11/16/01 

3.75E-07 

1 . 88E-07 

1.71E-07 

0.94 

1 

+ 

76.4% 

62A 

12/11/01 

6.34E-06 

3.14E-06 

2 . 06E-06 

0.86 

1 

+ 

58.7% 

63A 

64A 

12/11/01 

3.34E-06 

1 . 80E-06 

0.98 

2 

+ 

87 . 6% 

65A 

66A 

12/11/01 

2.43E-05 

4 . 90E-06 

6.71E-06 

0.87 

2 

+ 

92 . 1% 

Notes  : 

(1)  Units  for  hydraulic  conductivity  are  feet/second. 

(2)  Hydraulic  conductivity  for  early  and  late  data  computed  for  tests 
with  <  0.98. 

(3)  "Dev"  indicates  whether  the  test  was  conducted  on  an  undeveloped 

piezometer  (-)  or  on  a  piezometer  developed  on  11/08/2001  (1)  or 

11/28/2001  (2) . 

(4)  "Type"  indicates  whether  water  was  added  (+) ,  extracted  (-)  or 
displaced  (d) . 

(5)  "%Recov"  indicates  percent  recovery  relative  to  the  pre-test 
equilibrium  water  levels. 
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Middle  Soil  Layer 


Test  Date  Early  K  Middle  K  Overall  K 
12/15/01  too  fast 


11/15/01 


too  fast 
1.25E-08  1.000 


12/03/01  6.23E-06  1.48E-06  8.70E-07  0.934 


11/02/01 


1.33E-07  0.994 


12/07/01  6.39E-06  1.08E-06  l.lOE-06  0.788 
12/07/01  1.15E-05  7.89E-06  5.71E-06  0.960 
12/15/01  9.56E-07  4.83E-07  4.63E-07  0.973 


11/15/01 


1.26E-07  0.982 


12/15/01  3.76E-06  1.51E-06  8.24E-07  0.806 


12/15/01  4.22E-05  1.39E-05  4.64E-06  0.807 


11/15/01 

12/15/01 

12/15/01  5.14E-05 


5.67E-07  0.993 
4.40E-05  0.985 
3.44E-05  0.941 


12/15/01  1.76E-06  5.83E-07  4.96E-07  0.916 
12/15/01  1.85E-06  4.30E-07  1.93E-07  0.714 


11/15/01 


7.19E-08  0.978 


12/15/01  3.85E-06  1.31E-06  7.66E-07  0.804 


12/15/01  1.99E-05  7.64E-06  4.57E-06  0.784 


11/15/01 


4.50E-07  0.995 


12/16/01  8.15E-06  4.75E-06  3.58E-06  0.931 


10/18/01 

10/31/01 

11/05/01 


3.83E-06  0.999 
1.38E-06  0.993 
1.52E-06  0.987 


12/15/01  3.80E-05  2.12E-05  1.91E-05  0.945 


12/16/01  2 . 66E-06 


1.16E-07  0.357 


12/16/01  1.31E-06  2.63E-07  1.43E-07  0.938 
12/16/01  9.78E-06  5.02E-06  3.17E-06  0.891 


11/15/01 

12/16/01 

11/05/01 


1.44E-06  0.918 
1.28E-06  0.981 
7.79E-06  0.991 


Type  Recov 


Test  Date  Early  K  Middle  K  Overall  K  R 
12/04/01  2.37E-05  0.981 

12/04/01  1.16E-06  5.47E-07  6.36E-07  0.977 


Type  Recov 


11/15/01 


2.03E-08  0.999 


12/04/01  7.87E-07  5.73E-07  4.23E-07  0.885 


12/16/01  1.66E-05  5.94E-06  5.54E-06  0.959 


11/01/01 


3.53E-07  0.997 


12/16/01  7.54E-07  1.70E-07  1.93E-07  0.978 
11/05/01  2.49E-05  1.13E-05  1.23E-05  0.853 
12/16/01  7.38E-06  2.04E-06  1.96E-06  0.890 
12/05/01  1.39E-05  3.71E-06  2.82E-06  0.893 

12/16/01  2.44E-06  6.15E-07  5.72E-07  0.898 

12/16/01  1.63E-05  8.02E-06  3.93E-06  0.847 

12/05/01  2.59E-06  8.13E-07  3.85E-07  0.858 
11/13/01  4.50E-06  3.57E-07  1.03E-07  0.824 
12/05/01  6.01E-06  2.16E-06  1.54E-06  0.807 


11/13/01  3.15E-08 


1.45E-08  0.965 


49B 

12/11/01 

1.42E-05 

7.07E-06 

3.85E-06 

0.922 

2 

+ 

91 . 9% 

50B 

51B 

11/13/01 

6.73E-08 

0.991 

- 

+ 

49.0% 

51B 

12/11/01 

1 . 13E-05 

4.57E-06 

3.44E-06 

0.845 

2 

+ 

91 . 1% 

52B 

53B 

11/13/01 

8.37E-07 

5.23E-07 

6. 14E-07 

0.964 

- 

+ 

74 . 9% 

53B 

12/11/01 

2.44E-05 

1 .22E-05 

1.39E-05 

0.895 

2 

+ 

97 . 9% 

56B 

11/13/01 

9.17E-08 

3. 92E-08 

3. 99E-08 

0.890 

1 

+ 

24.3% 

56B 

12/05/01 

2 .20E-07 

3.09E-08 

0.826 

1 

+ 

33.7% 

57B 

58B 

11/13/01 

1.48E-08 

0.980 

- 

+ 

61 . 1% 

58B 

12/05/01 

1 . 96E-05 

l.OOE-05 

6. 85E-06 

0.901 

2 

+ 

99.2% 

58B 

12/05/01 

2 .21E-05 

1 . lOE-05 

8.31E-06 

0.893 

2 

+ 

97.7% 

59B 

12/05/01 

9.12E-06 

4 . 18E-06 

1 . 87E-06 

0.847 

1 

+ 

97 . 1% 

Piez 

Test  Date 

Early  K 

Middle  K 

Overall  K 

R^ 

Dev 

Type 

Recov 

61B 

12/11/01 

3.30E-06 

2 . lOE-06 

1 .20E-06 

0.952 

2 

+ 

96.3% 

62B 

63B 

12/11/01 

1 . 16E-06 

4.55E-07 

2 . 94E-07 

0.912 

2 

+ 

72 . 6% 

66B 

10/18/01 

6. 02E-07 

0.984 

- 

d 

31 . 9% 

66B 

11/05/01 

3.74E-06 

2 . 16E-06 

1 . 83E-06 

0.902 

- 

d 

82 . 6% 

66B 

12/11/01 

1.45E-05 

5.42E-06 

3.99E-06 

0.958 

2 

+ 

85.3% 

Notes  : 

(1)  Units  for  hydraulic  conductivity  are  feet/second. 

(2)  Hydraulic  conductivity  for  early  and  late  data  computed  for  tests 
with  <  0.98. 

(3)  "Dev"  indicates  whether  the  test  was  conducted  on  an  undeveloped 

piezometer  (-)  or  on  a  piezometer  developed  on  11/08/2001  (1)  or 

11/28/2001  (2) . 

(4)  "Type"  indicates  whether  water  was  added  (+) ,  extracted  (-)  or 
displaced  (d) . 

(5)  "%Recov"  indicates  percent  recovery  relative  to  the  pre-test 
equilibrium  water  levels. 
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Appendix  F:  Mathcad®  Template  for  Slug  Tests 


Piezometer:  05A 

Date:  12-07-2001 

Type :  Add  Water 


1.  Display  and  plot  the  test  data: 


^0.17^ 
0 . 33 
0 . 50 
0 . 67 
0 . 83 
1 

1 . 5 
2 

2 . 5 

3 

4 

5 

6 

7 

8 

10 

15 

25 

35 

V  50  , 


^24 .79  ^ 
23.19 
22 .09 
20 . 08 
18 . 50 
17 . 37 

16.28 

15.29 
14 .56 
13.90 
12 . 72 
11 . 68 
10 . 77 
9.92 

9.29 
8  . 18 
5.55 
3.35 
2 . 13 

V  1-37  , 


Equilibrium 

Water 

Level 


30 


20 


10 


0 


t 


F-1 


2.  Define  relevant  piezometer  data: 


Casing  radius:  rc:=. 25-in 

Intake  radius:  :=  .  5  •  in 

Effective  screen  length:  Lg:=  2.625-in 

Distance  from  piezometer  tip  to  center  of  screen:  tip  :=  2 . 1875 • in 

Piezometer  length  (tip  to  above-ground  coupling):  PL  :=  30 . 625 • in 

Midpoint  elevation  of  piezometer  screen:  Pgj_g^  :=  4 . 5943  •  f t 

Water  level  in  piezometer  down  from  top  of  tube:  wl  :=  .915 -ft 

Distance  from  equilibrium  water  level  to  bottom  of  screen  intake 
(D)  and  distance  of  equilibrium  water  level  above  the  impermeable 
wetland  liner  (H) : 

Le 

D:=PL-tip-wl  h:=Pgigv-  — 

H  :=  D  +  h 


3.  Determine  best  fit  line  and  compute  R2  for  semi-log  plot  of  y(t) 
versus  time.  This  is  done  by  using  built-in  Mathcad  functions  and 
adapting  the  principle  of  least  squares  as  presented  by  Devore  for 
an  intrinsically  linear  (in  this  case  exponential)  function  (see 
pages  497-506  &  552  of  Devore  (2000)  for  a  theoretical  explanation) . 

ORIGIN  =1  n  :=  rows  (t)  i  :=  1  ..  n 


Define  the  best-fit  line: 

Po  :=  exp  (  intercept  ( t  ,  In  (  y )  )  ) 

Pi't 

Pi  :=  slope  (t,ln(y))  line(t):=Po-e 


Compute  the  error  sum  of  squares: 


s:=ln(y)  -ln(line(t)) 


SSE:=  ^  (s,)" 

i  =  1 


F-2 


Compute  the  total  sum  of  squares: 


l^_Ybar  ■“ 


n 

SST;=  ^  (in(yi)  -  ln_y]3ar)^ 
i  =  1 


Compute  the  coefficient  of  determination  and  plot  the  best  fit 
line  against  the  observed  data: 

SSE 

'^squared  ^  ^squared  =  0.9405 


0  10  20  30  40  50 


t 


4.  Use  the  Bouwer  and  Rice  Method  to  Compute  K: 
Y0:=P0  Yn:=line(tn)  tgcaie-min 


Define  constants  "A"  and  "B" : ,  for 


=  5.25 


A:=1.75  B:=.25 


(Source:  figure  3.3.6,  Charbeneau,  2000) 


Determine  the  value  for  In [ (H-D) /r ^] : 

In  I  -  1  =  4.679  (Note :  if  ln[(H-D)/r  ]>6,  use  6.) 
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Determine  ln(Rg/r^)  and  define  as  "X"  : 


X  := 


f 


1 . 1 


A  +  B  •  in 


V 


<’  H  -  D 


)) 


in 


f  \ 


-  1 


(Equation  3.6.14,  p.l48, 
Charbeneau ,  2000) 


Compute  hydraulic  conductivity  K  : 


Kjjr  :— 


2-L, 


■  tscale 


■  in 


^  YO^ 


yYny 


(Equation  3.6.13,  p.l47, 
Charbeneau ,  2000) 


=  1.113x  10" 


ft 


sec 


5.  Use  Hvorslev's  Method  (Source:  Charbeneau) : 
Define  the  Shape  Factor  F: 


2 . 4 -m-L, 


F  := 


in 


1 . 2Lp 


2  -r.. 


+  ,  1  + 


(  1.2Le^ 

V  J 


(Equation  9,  Brand  and 
Permchitt,  1980) 


Compute  K: 


Khv  ;= 


■  tscale 


•in 


^  YO  ^ 


y  Yny 


(Equation  3.6.5,  p.l43, 
Charbeneau ,  2000) 


Khv  =  1.497  X  10  - 

sec 

Compare  Methods:  -  =  1.345 

Kjjr 
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6.  Determine  K  for  Early  Data: 


24 . 7892 


Compute  coefficient  of  determination  (R^) 


ORIGIN  =  I  n  :=  rows  ( t)  i  :=  I  ..  n 


Po  :=  exp  (  intercept  ( t ,  in  (  y)  )  ) 


Pi  :=  slope  (t,ln(y))  line(t)  :=Po-e 


s:=ln(y)  -  in  (line  (t))  SSE  :=  (^i)" 


23. 1892 
22 .0892 


20 . 0792 


18 .4992 
17 . 3692 


ln_ybar:=  ^  In(yi)' 


SST  := 


^  (in(yi)  -  ln_ybar)" 


'"squared  • 


^squared  “  ^  .  99286 


Compute  and  Compare  K  Values 


Yo:=Po  Yn:=line(tn) 


1  f  YO  1  - fi  ft 

Khr;= - in  -  Khr  =  8. 373  x  10  - 

2-Le  tn-tgcale  I  Yn  1  sec 


7t-(rc)‘ 


ft 


Khv- - in  -  Khv=  1-126  x  10  -  -  =  1.345 

1’  ta'tgQgi0  I  Yn  j  sec 
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7.  Determine  K  for  Middle  Data: 

^  1.50  ^ 

^  16.2792  ^ 

2 . 00 

15 . 2892 

2 . 50 

14 . 5592 

3.00 

13.8992 

4 . 00 

12.7192 

t  := 

5.00 

Y  ;= 

11 . 6792 

Compute  coefficient  of  determination  (R^)  : 

6.00 

10 .7692 

ORIGIN  =  I  n  :=  rows  ( t)  i  :=  I  ..  n 

7 . 00 

9.9192 

8 . 00 

9.2892 

Po  :=  exp  (  intercept  ( t ,  In  (  y)  )  ) 

10 .00 

8.1792 

^15.00  j 

^  5.5492  ^ 

Pi' t 

Pi  ;=  slope  (t,ln(y))  line(t)  :=Po'e 


s:=ln(y)  -ln(line(t)) 


ln_ybar:=  ^  In(yi)'  — 
i  =  1 


SSE:=  ^  (s,)" 

i  =  1 

n 

SST:=  ^  (in(yi)  -  ln_ybar)^ 
i  =  1 


SSE 

^^squared  ■=  1  “  1 

^squared  “  0  .  9 957 9 

YO 


Compute  and  Compare  K  Values  :  ^  ^  I _ ^ ^ _ I 

0  5  10  15 

YO:=Po  Yn;=line(tn) 


^br  ■ 


2-L, 


■  tscale 


■  In 


^  YO^ 


Yn 


=  1 . 497  X  10 


ft 

sec 


Khv  := 


•N‘ 


■  tscale 


■  In 


^  YO^ 


y  Yny 


Khv  =  2.015x  10 


-  6 


ft 


sec 


Khv 

-  = 1 . 345 

Kbr 
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8.  Determine  K  for  Late  Data: 


^25  ^ 

^3.3528  ^ 

t  := 

35 

Y  ;= 

2 . 1336 

,50; 

,1 . 3716  7 

Compute  coefficient  of  determination  (R2) : 

ORIGIN  =  1  n  :=  rows  ( t)  i  :=  1  ..  n 

Po  :=  exp  (  intercept  ( t ,  in  (  y)  )  ) 

Pi't 

Pi  :=  slope  (t,ln(y))  line(t)  :=Po-e 

n 

s:=ln(y)  -  in  (line  (t))  SSE  :=  (^i)^ 

i  =  1 

n  ^  n 

ln_ybar:=  ^  In(yi)-—  SST  :=  ^  (in(yi)  -  ln_ybar)^ 
i  =  1  i  =  1 


SSE 

I^squared  :=  1  “ 

^squared  =  0.98531 

Compute  and  Compare  K  Values  : 
Yo:=Po  Yn:=line(tn) 


Yi 


Yo 


1 


0.1 


0.01 


20  30  40  50 


ti 


K]3r 


rc^-X 

2  -Lf. 


■  in 


-n ■ Rscale 


^  YO^ 


yYny 


-  7  ft 

Kbr  =  6.733  X  10  - 

sec 


Khv  ;= 


F 


tn ■ tscale 


■  in 


^  YO^ 


yYny 


-  7  ft 

Khv=  9.059  X  10  - 

sec 


Khv 

-  = 1 . 345 

Kpr 
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Appendix  G:  Mathcad®  Template  for  Pump  Tests 


Pump  Test:  11-20-2001 

Pumping  Well:  21C 
Observation  Piezometer:  15C 


Hantush  Match-point  Method  for  a  Leaky  Confined  Aquifer 
with  Partially  Screened,  Partially  Penetrating  Wells  : 


1 .  Define  the  distance  between  the  pumping  and  observation 
wells  (r)  ,  the  aquifer  thickness  (b) ,  the  aquitard  thickness 
(b')/  and  the  pump  rate  (Q) . 


r  ;=  10 . 51139026  -ft 


ORIGIN  =  1 


b  :=  (1.5  +  1.25)  -ft 


tprime  ■—  3  '  ft 

1  •  gal 

Q:= - 

7 1 • sec 


2 .  Define  elevation  of  layer  surface  and  the  approximate 
piezometer  screen  length  (ASL) : 

Layer_C_elev-t-op  :=  2 . 75  •  f t 
ASL  ;=  2 . 625  •  in 

3.  Define  screen  midpoint  elevations  (SME)  for  the  pumping  and 
observation  well  and  compute  1^,  d^,  1^,  d^: 

;=  2 . 0337  -  ft 

ASL 

ip  :=  Layer_C_elev|-Qp  —  SMEp^jfip^-j^gPP  H  ^ 

ASL 

dp  ;=  Layer_C_elev|-Qp  —  SMEp^jfip^-j^gPP  ^ 


ip  =  0 . 826  ft 
dp  =  0 . 607  ft 
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SMEgbs  well  ■“  2 . 2  4  5  9  •  f  t 


ASL 

Ig  :=  Layer_C_elevtgp  -  SMEgbs^well  + 


ASL 

dg  ;=  Layer_C_elev|-gp  —  SMEgjgg^well  1 


4.  Estimate  K^/K^  and  compute  x: 


Kz  over  Kr  :=  1 


r 

X  :=  —  •  Kz  over  Kr 
b  ■  - 


5.  Define  the  range  of  values  for  "u^",  "1/u^",  and  "r/B^": 


(Note:  "Uj,"  spans  values  from 
.00001  to  10  but  is  not  defined 
in  the  margins . ) 

i  :=  1  ,  2  ..  rows 

1 

one_over_Ur  := - 

i  U  j;- 

i 

j:=l,2..  10 


6.  Define  the  Well  Function,  W(Uj,,r/B^),  for  partial  penetration: 

Note:  The  sumation  of  the  partial  penetration  correction  factor, 

n,  ranges  from  1  to  100  (instead  of  1  to  infinity  as  directed  in 
the  reference)  to  allow  solution  with  this  software  application. 


Ig  =  0 . 613  ft 

dg  =  0 . 395  f t 
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1 

Csl 

CM 

1 _ 

y  fyl 

2  /  \ 

2^  2 

!-■ 

1 

d 

!-■ 

0 

1 

d 

0 

1 

c 

\ — 1 

II 

(  ( 
s  in 
V  V 


n  •  7t  •  Ip 
b 


y 


( 

sin 

V 


n  •  71  •  dp 
b 


id 

)) 


(  ( 
sin 
V  V 


n  -TC  •  Iq 
b 


) 


n  -TC  •  dr 


id 

yy 


rco 


(r_over_Bj)^+  (n-7i-x) 


(equation  12.14,  Dawson 
and  Istok,  1991) 


7.  Define  the  Well  Function,  W(u^,r/B^),  for  partial  penetration: 


W 


1  ,  D 


1 

—  •  e 
Y 


(r_over_Bj^ 

4y 


dy  +  Fi , j 


(equation  12.12,  Dawson 
and  Istok,  1991) 


Display  the  time /drawdown  data  from  the  aquifer  test: 


Fobs  ■“ 


^  0.5  d 
0 . 75 
1 

1.5 
2 

2.5 

3 

3.5 

4 

4.5 
6 

13 

16 

V  20  y 


•min  Sobs  ■“ 


^0 . 31  ^ 
0 . 52 
0 . 75 
0 . 92 
1 .  14 
1 . 4 
1 . 53 
1 . 75 
1 . 95 
2 . 07 
2 . 31 
3.08 
3.29 

^3.45; 


■  cm 
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11.  Calculate  aquifer  hydraulic  conductivity: 


a.  List  match-points  for  W  (u, r/B)  and  drawdown,  s: 
'^star  •—  1  •  15  :=  .  108  •  ft 


b.  Compute  transmissivity,  T,  and  hydraulic  conductivity,  K: 


Q''^star  _  o  ft^ 

T  := -  T  =  2. 428x  10  - 

s star  sec 


(equation  9.8,  Dawson 
and  Istok,  1991) 


T  ft 

K:=—  K  =  8 . 828  X  10“ - 

b  sec 


12.  Calculate  aquifer  storativity,  S: 

a.  List  match-point  for  1/u  and  define  u* : 

1 

one_over_Ustar  :=  1 9  Ugtar  := - 

one_over_U3t3r 

b.  List  the  match-point  for  time,  t: 
tstar  :=  960  -sec 


c.  Compute  storativity,  S: 


(equation  9.10,  Dawson  and 
Istok,  1991) 


■  ^star  ■  ^star 

S  := -  S  =  4.441x  10 

2 

r 


G-6 


13.  Calculate  hydraulic  conductivity,  K',  for  the  overlying 
aquitard: 


a.  List  the  match-point  value  for  r/B  and  determine  the 
leakage  factor,  B: 

r 

r_over_Bg|-32-  :=  .  15  B  := - 

r_over_Bgtar 


b .  Compute  K ' : 


(equation  9.2,  Dawson  and  Istok,  1991) 


^prime  •“ 


'bprime  _  c  f  t 

7  ^prime  =  l-  483x  10 


B 


sec 


14.  Determine  hydraulic  conductivity  for  the  bottom  soil  layer 
of  the  wetland: 

a.  Estimate  gravel  and  sand  layer  conductivities  and  define 
layer  thicknesses  (source  of  estimates:  Charbeneau,  2000) : 

0 

bgravel  ■“  •  15  •  ft  1  x  10 

_  2  cm 

^sand  ■=  •  5  •  ft  '^sand  ■=  1  ^  10 


b.  Define  soil  layer  thickness  and  compute 


^soil  •—  1  •  5  •  ft 


(Equation  2.2.10, 
Charbeneau,  2000.) 


_  4  f  t 

Kgoii  =  9.566  X  10  - 

sec 


bsoil 


^soil 


b  bggnfp  bgravel 
^  ^sand  '^gravel 
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Hantush  Match-point  Method  for  Leaky  Confined  Aquifers 
with  Aquitard  Storage  : 


Note:  this  method  assumes  the  effects  of  partial  penetration 
introduce  little  error  when  r  >  1 . 5m  (Kj^/K^)  1/2  (p.  199,  Dawson 

and  Istok,  1991)  and  can  be  ignored. 

1.  Test  Range  of  versus  that  assumption  is  valid  over: 

b:=  (1.5+1.25)  -ft  r:=  10. 51139026- ft 

assume:  :=  6  Kz  :=  1 

Kj. 

^test  •“  1 • 5  -b •  ^test  “  10.104  ft 

^7. 


2.  Define  values  for  "u"  and  "p": 


P:= 


^0.01'' 

0 .05 
0 . 1 
0.2 
0 . 5 
1 
2 
5 

10 

V  20  , 


Note:  "u"  spans  values  from  .000001  to  3  but 
is  not  defined  in  the  margins. 
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3.  Define  the  range  of  variables  for  the  Well  Function,  H(u,  p),  for 
Leaky  Aquifers  with  Aquitard  Storage: 


j  :=  1  ,  2  ..  10  i  :=  1  ,  2  ..  rows  (u) 


4.  Define  1/u: 


1 

one_over_uj_  := - 

Ui 


5.  Define  the  Well  Function,  H(u,p)  : 


•  erf  c 


y 


(y  -  Ui) 


dy 


(equation  10.18,  Dawson 
and  Istok,  1991) 


6.  Define  relevant  test  and  aquifer  parameters  to  support  below 
calculations : 

a.  Distance  between  pumping  well  and  observation  well: 
r  :=  10 . 51139026  -ft 


b.  Aquifer  thickness: 

b  :=  (1.5  +  1.25)  -ft 


c.  Aquitard  thickness: 

bprime  :=  3  •  f t 

d.  Pump  Rate : 

1  •  gal 

Q  := - 

7 1 • sec 
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9.  Calculate  hydraulic  conductivity  for  the  bottom  soil  layer: 


a.  List  the  match-points  H(u,(3)*,  1/u*,  s* ,  and  t*  and  define  u* 


1 

^star  •=  1  •  6  one_over_u3t;g2-  :=  29  ^31^32-  := 

one_over_U3|-32- 


Sstar:=3.08-cm  t3tar  :=  13 -min 


b.  Compute  transmissivity  (T),  hydraulic  conductivity  (K) ,  and 
storativity  (S) : 


Q ■ ^star 

T  := - 

4  ■  ■  s  star 


T  =  2 .372  X 


sec 


(equation  10.15,  Dawson 
and  Istok,  1991) 


T  ft 

K:=—  K  =  8 . 627  X  10“ - 

b  sec 


“1  ■  1  ■  ^star  ■  Ustar 

S  := -  S  =  2.31x  10 

2 

r 


(equation  10.16,  Dawson 
and  Istok,  1991) 


c.  Estimate  gravel  and  sand  layer  conductivities  and  define 
layer  thicknesses  (source  of  estimates:  Charbeneau,  2000): 

0 

bgravel  •“  •  7 5  •  f  t  ^Igravel  •—  1  ^  10 

_  2  cm 

l^sand  •=  •  5  •  ft  ^^sand  •=  1  ^  10 

sec 

bsoil  •=  1  •  5  •  ft 


d.  Determine  the  hydralic  conductivity  of  the  bottom  soil 

layer  of  the  wetland  (K  .,) : 

'  soil 


^soil 


b 

K 


bsoil 


^sand  bgravel 
^sand  ^Igravel 


_  4  f  t 

Ksoii  =  9.142  X  10  - 

sec 

(Equation  2.2.10, 
Charbeneau,  2000.) 
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Appendix  H:  Calibration  Plots 


Calculated  vs.  Observed  Heads  :  Steady  state 


Obs.  Heads  (ft) 


Num.RDints  :  63 
Mean  Error  :  0.04819473  (ft) 
Mean  Absolute  :  0.06677023  (ft) 
Standard  Error  of  the  Estimate  :  0.01456262  (ft) 
Root  mean  squared  :  0.1243828  (ft) 
Normalized  RMS  :  6.300094  (  %  ) 


Calculated  vs.  Observed  Heads  :  Steady  state 


7.63  7.73  7.83  7.93  8.03 


Obs.  Heads  (ft) 

Num.RDints  :  64 


Mean  Error  :  0.007224859  (ft) 
Mean  Absolute  :  0.02538688  (ft) 
Standard  Error  of  the  Estimate  :  0.007304927  (ft) 
Root  mean  squared  :  0.05842946  (ft) 
Normalized  RMS  :  13.26435  (  %  ) 
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